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I.  BACKGROUND  AND  PURPOSE 

This  thesis  details  the  development  of  user  friendly 
mathematical  applications  capable  of  computing  the  radiation 
pattern  and  other  pertinent  antenna  parameters  of  an  antenna 
or  antenna  system  based  on  available  information,  primarily 
dimensional  information  obtained  from  photographs.  The 
applications  developed  are  compatible  with  any  IBM  personnel 
computer  using  MS-DOS  version  3.2  or  higher  and  a  math 
coprocessor.  The  applications  are  written  in  Mathcad  version 
3.1  engineering  software  and  can  be  run  using  either  Mathcad 
3  . 1  or  Mathcad  version  4.0.  These  applications  are  intended 
for  use  by  the  Naval  Maritime  Intelligence  Center 
(NAVMARINTCEN) . 

As  already  mentioned,  user  inputs  are  limited  to  antenna 
dimensions  based  on  physical  measurements  primarily  obtained 
from  photographs,  although  other  source  data  can  be  used.  The 
applications  are  flexible  to  the  extent  that  other  parameters 
may  be  estimated  and  used  as  inputs  to  increase  the  accuracy 
of  the  computations  to  gain  better  insight  into  the  antenna's 
performance.  The  Mathcad  applications  provide  various 
performance  predictions  as  well  as  a  graphical  representation 
of  the  antenna's  far-field  radiation  pattern.  The  necessary 
background  information  needed  to  interpret  the  application's 


formulas  and  displays  are  provided  in  the  corresponding  thesis 
chapters. 

Dietrich  [Ref.  1]  and  Gerry  [Ref.  2]  have  completed  the 
first  and  second  reports  of  this  project,  respectively.  This 
thesis  is  the  last  in  a  series  of  three  reports  intended  to 
fulfill  the  NAVMARINTCEN  statement  of  work. 


II.  INTRODUCTION 

The  increasingly  sophisticated  design  by  foreign 
countries  of  antennas  warrants  the  development  of  improved 
analytical  methods  to  provide  timely  and  accurate  technical 
assessment  of  these  antennas.  Through  the  process  of  reverse 
engineering,  foreign  communication,  navigation, 
Identification,  Friend  or  Foe  (IFF) ,  and  radar  antennas  are 
evaluated  to  determine  their  capabilities,  limitations,  and 
vulnerabilities.  Having  limited  information  available 
concerning  the  parameters  of  an  antenna,  usually  restricted  to 
photographs  that  provide  only  the  type  and  physical  dimensions 
of  the  antenna,  current  methods  of  antenna  evaluation  are  slow 
and  tedious. 

With  the  availability  of  high-speed  personal  computers 
matched  with  sophisticated  off-the-shelf  engineering  software, 
the  assessment  of  antennas  using  only  photographic  information 
may  now  be  obtained  rapidly.  The  objective  of  this  report  is 
to  provide  NAVMARINTCEN  with  user  friendly  Mathcad  software  to 
achieve  this  type  of  performance  in  their  evaluation  of 
antenna  systems. 

As  with  any  computer  analysis  of  antennas,  the  accuracy 
of  the  results  depend  upon  the  sophistication  of  software 
used,  the  complexity  of  calculations,  the  reliability  of  input 
data  obtained  from  photographs,  and  the  accuracy  of  any 


required  input  data  that  must  be  estimated.  This  report  and 
the  Mathcad  applications  are  written  with  existing  engineering 
equations  for  the  type  of  antenna  analyzed  and  inputs  are 
limited  to  physical  dimensions.  Additional  information  such 
as  ground  effects  and  characteristic  impedance  may  be  included 
to  increase  the  accuracy  of  the  analysis.  To  this  end,  this 
report  should  provide  an  initial  understanding  as  to  the 
capabilities  of  the  antenna  system  in  question. 

The  chapters  of  this  report  document  the  equations  and 
assumptions  used  for  each  specific  antenna  type  and  is  written 
as  a  comprehensive  reference  for  the  software.  Copies  of  each 
application  are  included  as  appendices  to  provide  the  user 
with  a  printed  illustration  of  the  software. 


III.  THE  UNIFORM  LINEAR  ARRAY  ANTENNA 

Arrays  are  highly  directional  antenna  system  with  narrow, 
steerable  beams  and  low  side  lobes  used  for  long  distance 
communication  and  radar  systems.  The  linear  array  is  one 
geometrical  configuration  that  consists  of  radiating  elements 
lying  along  a  straight  line.  As  a  result  of  the  dependance  of 
inter-element  spacing  (d)  to  the  wavelength  (A)  ,  linear  arrays 
are  inherently  narrowband  antenna  systems. 

In  this  Mathcad  application  uniformly  spaced,  equally 
excited,  isotropic  point  sources  aligned  on  the  z-axis  are 
used  as  the  radiating  elements.  In  using  isotropic  point 
sources,  polarization  of  the  array  is  not  calculated  and  an 
assumption  is  made  that  the  antenna  input  resistance  (Rin)  is 
equal  to  the  radiation  resistance  (Rr)  .  This  assumption 
results  in  calculations  for  gain  (G)  that  are  idealized 
values.  The  following  inputs  are  needed  to  implement  the 
software: 

N  =  number  of  isotropic  radiating  elements 

d  =  inter-element  spacing  between  adjacent  elements 

6  =  direction  of  main  lobe 

0 

f  =  frequency  of  interest 
I  =  antenna  feed  current 

o 

Zo  =  characteristic  feed  impedance 
The  first  two  inputs  are  obtained  from  visual  data  such  as 


photographs,  and  the  indicates  input  parameters  that  are 
either  known  or  are  estimated.  The  linear  array  geometry  is 
illustrated  in  Figure  3.1. 


FIGURE  3.1  Geometry  of  N-element  linear  array  of 
isotropic  point  sources 


The  total  electric  field  of  a  linear  array  in  the  far- 
field,  neglecting  mutual  coupling  between  adjacent  elements, 
is  determined  by  the  vector  addition  of  the  fields  of 
individual  radiating  elements.  It  can  be  shown  that  the  total 
field  is  equivalent  to  the  product  of  the  field  of  a  single 
element,  selected  at  a  reference  point,  and  a  factor  which  is 
referred  to  as  the  array  factor  (AF) : 

E(  total)    =  [E  (element)  ]   x  [AF]  (3.1) 

In  (3.1),  E  is  the  vector  electric  field  intensity,  while  AF 
is  a  scalar  quantity.  This  result  is  referred  to  as  pattern 
multiplication  for  arrays  of  identical  elements.  The  array 
factor  is  a  function  of  the  geometry  of  the  array  and  the 
element  excitation  amplitude  and  phase.  For  uniform 
excitation,  by  varying  the  inter-element  spacing  (d)  and/or 
the  phase  (/?)  between  the  elements,  the  characteristics  of  the 
array  factor,  and  consequently  the  total  electric  field,  of 
the  array  can  be  controlled.   [Ref.  3:  pp.  204-207] 

This  Mathcad  application  is  written  for  an  N-element 
linear  array  with  equally  spaced,  uniformly  excited  isotropic 
radiating  elements.  In  addition,  element  phase  is  assumed  to 
vary  linearly  along  the  length  of  the  array.  The  radiation 
pattern  of  a  specific  element  type  is  neglected  in  (3.1)  since 
in  normal  usage  it  will  have  little  effect  when  the  array 
consists  of  a  large  number  of  elements.   However,  the  user  of 


this  software  should  be  able  to  use  information  regarding  the 
type  of  element  and  the  orientation  of  the  elements  to  narrow 
the  range  of  possible  6's.  The  normalized  array  factor  (AFp) 
is  [Ref.  3:  pp.  212-214]: 


^<e)  =i 


■l-f  *  (9) 


sin(-|i|r  (6) 


(dimensionless) 


(3.2) 


where 


i|/(6)  =  kdcos(Q)    +p    (radians) 


(3.3) 


P  =  -icdcos(60)    (radians) 


(3.4) 


In  (3.4),  6  is  the  direction  of  the  maximum  array  factor 
value  (i.e.,  for  a  broadside  array  0o  =  tt/2,  for  an  end-fire 
array  with  the  main  beam  directed  at  180°,  0o  =  w ,  and  for  a 
phased  array  with  the  main  beam  directed  to  ±60°,  6o  =  ir/3)  , 
(3  is  the  phase  excitation  difference  between  the  elements,  k 
is  the  wavenumber  (2ir/X)  ,  and  d  is  the  inter-element  spacing. 
With  the  linear  array  aligned  along  the  z-axis,  the  array 
factor  is  a  function  of  6  but  not  of  <p.  The  visible  region  of 
the  array  in  terms  of  8,  ij»(6),  and  the  inter-element  spacing 
(d/X)    is: 

0  <  6  <  7i    (radians)  (3.5) 


P-kd    <  t|/(6)  <  p+£d    (radians)  (3.6) 

d/A  determines  how  much  of  the  array  factor  appears  in  the 
visible  region  as  defined  by  (3.5)  and  (3.6).  Although  d/X 
provides  guidance  in  determining  the  frequency  to  select  for 
this  Mathcad  application,  it  should  be  noted  that  most  arrays 
are  designed  with  the  inter-element  spacing  less  than  one 
wavelength,  usually  close  to  a  half -wavelength.  Thus,  linear 
arrays  are  characteristic  narrowband  antenna  systems.  [Ref. 
4:  p.  123] 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

i  z.   1.6  A.    {meters)  (3.7) 

i  >  5  (Nd)         (meters)  (3.8) 


r  £ 


2  {"d)  2    (meters)  (3.9) 


The  array  length,  Nd,  in  (3.8)  and  (3.9)  is  the  maximum 
dimension  of  the  antenna  (the  array  length  is  assumed  to 
include  a  distance  of  d/2  beyond  each  end  element  [Ref.  5: 
p. 55]).  The  minimum  distance  to  the  far-field  is  found  by 
selecting  the  largest  value  of  (3.7),  (3.8),  and  (3.9). 

Radiation  intensity   (U(6))   is  related  to  the  power 
radiated  in  a  given  direction  and  is  independent  of  the 


distance  to  the  observation  point  [Ref.  3:  p.  27;  Ref.  4:  p. 
33].   For  linear  arrays  U(0)  is  given  by  [Ref.  3:  pp.  229- 
233]  : 

U(B)    =  (AFn)2         (W/  solid  angle)  (3.10) 

Radiated  power  (Prad)  is  the  total  power  radiated  by  the 
antenna  and  is  obtained  by  integrating  U(0)  over  a  surface 
surrounding  the  antenna  [Ref.  4:  p.  33].  Since  the  geometry 
of  the  linear  array  is  aligned  on  the  z-axis  and  the  array 
factor  is  independent  of  <p ,    Ppad  is  thus  [Ref.  3:  pp.  229-233]: 

Prad  =  2tz    fnU(G)  sin(Q)  dQ        (W)  (3.11) 

Jo 

Directivity  (Do)  is  the  maximum  value  of  directive  gain 
where  directive  gain  is  the  ratio  of  U(6)  in  a  specific 
direction  to  the  average  radiation  intensity  [Ref.  4:  pp.  34- 
36],   Hence,  DQ  is  given  by  [Ref.  3:  pp.  229-233]: 

D    =     max    (dimensionless)  (3.12) 

^rad 

In  (3.12),  U|nax  equals  unity  and  occurs  at  6o  for  the  uniform 
linear  array. 

Effective  isotropic  radiated  power  (EIRP)  is  a  commonly 
used  communications  term  that  is  defined  as  the  product  of 
antenna  gain  and  total  power  accepted  by  the  antenna  from  the 
transmitter.   EIRP  is  determined  as  [Ref.  4:  p.  62]: 
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EIRP   =  PraHDn  (W)  (3.13) 


Assuming  that  the  power  delivered  to  the  array  consisting 
of  N-isotropic  point  sources  is  equivalent  to  the  radiated 
power,  we  can  estimate  the  radiation  resistance  (Rp)  of  the 
antenna  as  [Ref.  3:  p.  55]: 


Rr   =        ■  Ia°  (Q)  (3.14) 


o 


Using  the  same  logic,  we  obtain  the  input  resistance  (Rin)  of 
the  array  made  up  of  isotropic  point  sources  as  equivalent  to 
the  radiation  resistance: 

Rin  =  Rr         (Q)  (3.15) 

The  exact  gain  (G)  of  the  array  cannot  be  calculated  with 
the  data  assumed,  but  an  idealized  value  can  be  computed  based 
on  the  assumption  made  to  determine  Rjn  and  the  assumptions 
made  to  calculate  total  antenna  efficiency  (et)  .  The  total 
antenna  efficiency  is  given  by  [Ref.  3:  pp.  44-45]: 

et  =  titcd  =   l  -  i r I2    (diwensionless)  (3.16) 


where 


z^  =  1         (dimensionless)  (3.17) 


er  =  l-|r|2        (dimensionless)  (3.18) 
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r  =  -^ ^         (dimensionless)  (3.19) 

Rin    +    Zo 


In  (3.17),  ecd  is  the  conduction  and  dielectric  efficiency 
which  is  unity  for  an  ideal  lossless  antenna.  In  (3.19),  Zo 
is  the  characteristic  impedance  of  the  transmission  line  and 
r  is  voltage  reflection  coefficient.  Therefore,  the  total 
antenna  efficiency  is  equal  to  the  mismatch  efficiency  of 
(3.18)  and  the  ideal  gain  for  a  linear  array  is  [Ref.  3:  pp. 
43-44] : 

G  =  ttD0        (dimensionless)  (3.20) 

G(dB)    =   10log10  (ttD0)  (dB)  (3.21) 

As  a  result  of  constructing  the  linear  array  with 
isotropic  point  sources,  the  polarization  loss  factor  (PLF) 
and  antenna  polarization  are  not  calculated.  Therefore,  the 
maximum  effective  aperture  (Aem)  of  the  array  when  used  as  a 
receiver  is  calculated  assuming  negligible  polarization 
mismatches.   The  Aom  is  then  [Ref.  3:  p.  631: 

em  L  c  j 

Km   =  G[X)2  PLF         (m2)  (3.22) 


where 


PLF  =  1         (dimensionless)  (3.23) 
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The  assumptions  used  in  finding  Rr  and  Aem  result  in  an 
idealized  value  in  calculating  the  effective  height  (h^)  .  The 
effective  height  is  [Ref.  6:  p.  42]:] 


Km  =  2 


\ 


■^ss    (meters)  (3.24) 

1o 


where  rjo   is  the  characteristic  impedance  of  free  space  (12  07T). 

Two  displays  are  produced  in  this  Mathcad  application. 
The  first  display  is  the  rectangular  representation  of  the 
magnitude  of  the  array  factor  from  -n  to  n .  The  second 
display  is  the  polar  representation  of  the  magnitude  of  the 
array  factor. 

The  Hansen-Woodyard  end-fire  array  was  not  addressed  in 
this  application  nor  were  arrays  of  uniform  spacing  but 
nonuniform  amplitude  (i.e.,  Dolph-Tschebyschef f  array).  In 
addition,  nonuniformly  spaced,  uniformly  excited  arrays  were 
not  addressed  by  this  application. 

The  results  obtained  with  this  application  are  identical 
to  those  found  in  [Ref.  3:  pp.  216-240]  with  respect  to  array 
factor  patterns  and  directivity.  Table  3.1  is  a  comparison  of 
data  from  [Ref.  3]  for  a  10-element  equally  excited,  uniformly 
spaced,  end-fire  linear  array  while  Figure  3.2  is  a  comparison 
of  the  array  factor  pattern  of  the  same  array.  Of  note,  the 
array  factor  patterns  from  this  application  and  the  patterns 
from  [Ref.  3]  are  identical  as  expected. 
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Table  3.1   End-fire  Array  Data  Comparison 
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FIGURE  3.2  Array  factor  pattern  of  10-element  equally 
excited,  uniformly  spaced  end-fire  linear  array  with  d  = 
A./4  m,  X0   =  0°  ,  and  f3   =   -kd 
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Table  3.2  is  a  comparison  of  data  from  [Ref.  3]  for  a  10- 
element  equally  excited,  uniformly  spaced,  broadside  linear 
array,  and  Figure  3.3  is  a  comparison  of  the  array  factor 
pattern  of  the  same  array.  The  array  factor  patterns  from 
this  application  and  the  patterns  from  [Ref.  3]  are  identical. 

Table  3.2   Broadside  Array  Data  Comparison 


ANTENNA 

REFERENCE 

CALCULATED 

DIMENSIONS 

DIRECTIVITY  (DQ) 

DIRECTIVITY  (DQ) 

N  =  10 

d  =  A./ 4  m 

6.99  dB 

7.14  dB 

6o  =  7T/2 

(3    =    0 
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FIGURE  3.3  Array  factor  pattern  of  10-element  equally 
excited,  uniformly  spaced  broadside  linear  array  with  d  = 
X/4  m,  60  =  it/2,    and  (3   =   0 
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IV.  THE  PLANAR  ARRAY  ANTENNA 

Planar  arrays,  frequently  used  long  distance 
communications  and  radar  systems,  exhibit  characteristics 
analogous  to  those  of  linear  arrays  but  have  additional 
capability  to  control  and  shape  the  radiation  pattern.  Planar 
arrays  are  more  versatile  than  linear  arrays,  providing 
symmetrical  patterns  with  lower  side  lobes,  and  can  be  used  to 
scan  the  main  beam  of  the  antenna  array  toward  any  point  in 
space  [Ref.  3:  p.  261].  As  with  linear  arrays,  planar  arrays 
are  narrowband  antenna  systems  due  to  the  dependance  of  the 
inter-element  spacing  (dy  and  d)    on  the  wavelength  (A) .   The 

x         y 

following  inputs  are  required  analyze  a  planar  array: 
M  =  number  of  isotropic  radiating  elements  in  the 

x-direction 
N  =  number  of  isotropic  radiating  elements  in  the 

y-direction 
dx  =  inter-element  spacing  between  adjacent  elements  in 

the  x-direction 
dy  =  inter-element  spacing  between  adjacent  elements  in 

the  y-direction 
f  =  frequency  of  interest 
6n  =  direction  of  main  lobe  at  6  =  0 

o  o 

0O  =  direction  of  main  lobe  at  <p  =   0o 
Io  =  antenna  feed  current 
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ZQ  =  characteristic  feed  impedance 
The  first  four  inputs  are  data  obtained  primarily  from  visual 
data  such  as  photographs,  and  the  *  indicates  input  parameters 
that  are  either  known  or  estimated.  The  planar  array  geometry 
is  illustrated  in  Figure  4.1. 


z 

i 
X 

e    / 

4fN   >y 

FIGURE  4.1  Geometry  of  M  x  N  element  planar  array 
of  isotropic  point  sources 
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This  Mathcad  application  is  written  for  an  M  x  N  planar 
array  with  uniformly  spaced,  equally  excited,  identical 
radiating  elements  positioned  on  the  x-y  plane.  This 
application  may  also  be  used  to  approximate  the  array  factor 
and  field  patterns  for  arrays  that  have  elliptical  or  circular 
geometries  by  assuming  an  ellipse  is  equivalent  to  a  rectangle 
or  a  circle  is  equivalent  to  a  square,  respectively  [Ref.  6: 
p. 187],  Since  isotropic  point  sources  are  assumed  as  array 
elements,  the  polarization  of  the  array  is  not  calculated,  and 
an  assumption  is  made  that  the  antenna  input  resistance  (Rjn) 
is  equivalent  to  the  radiation  resistance  (Rr)  .  This 
assumption  results  in  values  for  gain  (G)  that  are  idealized. 
As  with  the  linear  array,  the  radiation  pattern  of  a  specific 
element  type  is  neglected  in  the  analysis  of  the  planar  array 
since  in  normal  usage  it  will  have  little  effect  for  a  planar 
array  with  a  large  number  of  elements.  However,  as  with  the 
linear  array,  the  user  of  this  software  should  be  able  to  use 
information  regarding  the  type  of  element  and  the  orientation 
of  the  elements  to  narrow  the  range  of  the  possible  0o,s  and 

The  normalized  array  factor  (AFn)  is  calculated  for  far- 
field  observations  neglecting  mutual  coupling  between  adjacent 
elements  with  equal  amplitude  excitation  and  is  [Ref.  3:  pp. 
260-263] : 
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AFn(0,4> 


sin/^tx(ef*) 


M  sin 


f  *x(8f<J>) 


Hi! 


8in(-S*v<8,*J 


N  sin  -I£ 


icje,*) 


(4.1) 


where 


\|rx(6,4>)    =  Jcdxsin(0)  cos  (<|>)    +  Px        (radians) 


(4.2) 


11/(0,(1))  =  Jed  sin(  0)  sin(4>)  +  P„    (radians) 


(4.3) 


Px  =  -Jcdxsin(0o)  cos  ( 4>0)    (radians) 


(4.4) 


p  =  -icd sin(80)  sin(4>0)    (radians) 


(4.5) 


Equation  (4.1)  is  the  normalized  array  factor  that  is  a 
function  of  both  6  and  0  where  (4.2)  and  (4.3)  are  the  array 
factor  phase  shifts  in  the  x  and  y  directions,  respectively. 
(4.4)  and  (4.5)  are  the  progressive  phase  shift  between 
adjacent  elements  in  the  x  and  y  directions,  respectively, 
with  the  main  beam  directed  along  0  =  0o  and  <p  =  <po.  In  (4.2) 
through  (4.5),  k  is  the  wavenumber  (2ir/X)  . 

When  the  inter-element  spacing  between  the  elements  (dx 
and  d  )  equals  or  is  greater  than  the  wavelength  (k)  ,  multiple 
maxima  of  equal  magnitude  are  formed.  To  avoid  creating 
multiple  maxima,  the  inter-element  spacing  must  be  less  than 
a  wavelength  at  the  frequency  of  operation  [Ref.  3:  p.  263]. 
With   the   dependance   of   inter-element   spacing   on   the 
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wavelength,   planar  arrays   are   characteristic   narrowband 
antenna  systems . 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

r  z   1.6  A.    (meters)  (4.6) 

r  ^  5 y  ( Mdx ) 2  +  ( Ndy ) 2    (meters)  (4.7) 


2((Mdx)2   +  (iVdy)2  )  ,   B1 

r  ^  — i - r ^ 1    (meters)  (4.8) 


The  length  of  an  array  is  assumed  to  include  a  distance  d/2 
beyond  each  end  element  [Ref.  5:  p.  55].  Therefore,  the 
maximum  dimension  of  the  planar  array  is: 


J(Mdx)2   +  (Ndy)2         (meters)  (4.9) 

The  minimum  distance  to  the  far-field  is  found  selecting  the 
largest  value  of  (4.6),  (4.7),  and  (4.8). 

Directivity   (DQ)   for   planar   arrays   is   determined 
approximately  by  [Ref.  7:  p.  78]: 

n      -   A  tt  A-FnAFn  Imax  /  tat\ 

Do    -    47t     r2nrn/2/ " ~ {W)  (4.10) 

Jo  J       (AFnAF^)sin(d)   d6d<t) 

The  maximum  value  of  the  numerator  in  (4.10)  is  unity  at  0  = 
6o  and  0  =  0o  since  the  array  factor  is  normalized.  The 
denominator  is  the  radiated  power  (Prad)  of  the  antenna  where 
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Prad  ^S  obtained  by  integrating  over  all  angles  around  the 
antenna  above  the  x-y  plane  (0  <  0  <  tt/2)  [Ref.  3:  p.  28;  Ref. 
7:  p.  78].  This  assumes  the  array  is  always  pointed  toward 
one  side  or  the  other  of  the  x-y  plane. 

Effective  isotropic  radiated  power  (EIRP)  for  planar 
arrays  is  defined  as  the  product  of  antenna  gain  and  the  total 
power  radiated.   EIRP  is  determined  as  [Ref.  4:  p.  62]: 

EIRP  =  PradD0         (W)  (4.11) 

Assuming  that  the  power  delivered  to  planar  arrays 
constructed  of  M   x  N   isotropic   radiating   elements   is 
eguivalent  to  the  radiated  power,  we  estimate  the  radiation 
resistance  (R  )  of  the  array  as  [Ref.  3:  p.  55]: 


?  (  P      ) 
Rr  =  Z  (Q)  (4.12) 

\I/ 


Since  the  radiating  elements  are  assumed  to  be  ideal  point 
sources,  the  input  resistance  of  a  planar  array  (R,n)  is 
assumed  to  be  eguivalent  to  the  radiation  resistance: 

Rin  =  Rr         (Q)  (4.13) 

The  gain  (G)  of  a  planar  array  cannot  be  calculated  given 
the  information  assumed  to  be  available,  but  an  idealized 
value  can  be  derived.  Another  factor  affecting  the  gain  is 
the  assumption  made  when  determining  the  total  antenna 
efficiency  (et) .   The  total  antenna  efficiency  for  an  antenna 
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is    [Ref.    3:    pp.    44-45]: 

ec  =  erecd  =  l-|r|2         (dimensionless)  (4.14) 


where 


zrd  =  1    (dimensionless)  (4.15) 


er  =  l-|r|2    (dimensionless)  (4.16) 


T  =  Rin    Z°         (dimensionless)  (4.17) 


In  (4.14),  ecd  is  the  conduction  and  dielectric  efficiency, 
which  is  unity  for  an  ideal  lossless  antenna.  In  (4.17),  Zo 
is  the  characteristic  impedance  of  the  transmission  line  and 
r  is  the  voltage  reflection  coefficient.  Therefore,  the  total 
antenna  efficiency  is  equal  to  the  mismatch  efficiency  of 
(4.15)  when  ideal  lossless  antenna  is  assumed.  The  ideal  gain 
of  a  planar  array  in  dimensionless  and  decibel  quantities  is 
then  [Ref.  3:  pp.  43-44]: 

G  =  ttD0        (dimensionless)  (4.18) 

G(dB)    =  101og10  (ttD0)         (dB)  (4.19) 

The  maximum  effective  aperture  (Aem)  for  a  planar  array 
is  derived  by  assuming  polarization  mismatches  are  negligible 
(this  application  is  constructed  using  isotropic  radiating 
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elements  for  which  polarization  loss  factor  (PLF)  and  antenna 
polarization  are  not  calculated)  [Ref.  3:  p.  63]: 

»   _  CU)2 


4tt 


PLF         ( m2  ) 


(4.20) 


PLF  =  1    {dimensionless) 


(4.21) 


The  effective  height  (hem)  of  the  planar  array  is  now 
estimated  by  [Ref.  6:  p.  42]: 


Km  =   2 


\ 


RrKm 


(m) 


(4.22) 


In  (4.22),  r)0  is  the  characteristic  impedance  of  free  space 

(1207T)  . 

Four  two-dimensional  array  pattern  displays  are  produced 
in  this  Mathcad  application:  the  first  is  in  the  x-z  plane; 
the  second  is  in  the  y-z  plane;  the  third  display  generated  is 
a  plane  perpendicular  to  the  x-y  plane  when  <p  -  7r/4;  the  last 
display  is  a  plane  perpendicular  to  the  x-y  plane  when  <p  =  <p0. 
The  last  display  is  redundant  when  0  is  either  0,  7r/4,  or  ir/2. 

Results  obtained  with  this  application  are  the  same  as 
those  found  in  [Ref.  3:  pp.  260-274].  Directivity  was 
calculated  using  the  approximate  eguation  (4.10).  An 
alternative  approximation  for  directivity  is  given  by  either 
[Ref.  3:  p.  272-273]: 
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DQ  =  n  cos  (0O)  DxDy        (dimensionless) 


(4.23) 


or 


D, 


JC' 


QA (rads2) 


(dimensionless) 


(4.24) 


In  (4.23),  Dx  and  Dy  are  the  directivities  of  broadside  linear 
arrays  of  length  and  number  of  elements  Mdx,  M,  and  Ndy,  N  in 
the  x  and  y  directions,  respectively.  In  (4.24),  nA  is  the 
beam  solid  angle. 

Table  4.1  is  a  comparison  of  directivity  computed  for  a 
10x10  planar  array.  The  reference  directivity  is  derived 
using  (4.24)  while  the  calculated  directivity  is  from  (4.10). 

Table  4.1  Planar  Array  Data  Comparison 


ANTENNA 

REFERENCE 

CALCULATED 

DIMENSIONS 

DIRECTIVITY  (D0) 

DIRECTIVITY  (D0) 

M  =  N  =  10 

dx  =  dy  =  X/2   m 

0O  =  30° 

23.67  dB 

24.07  dB 

0o  =  45° 

/3X  =  (3y    =    -tt/(2J2) 

Figure  4.2  and  Figure  4.3  are  comparisons  of  array  factor 
patterns  for  a  5x5  planar  array  with  patterns  from  [Ref.  3] 
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and  the  calculated  patterns.  Figure  4.2  is  a  representation 
of  the  array  factor  patterns  in  the  x-z  and  y-z  planes,  and 
the  scale  of  this  figure  has  been  increased  in  order  to 
increase  the  clarity  of  the  pattern.  Figure  4.3  is  the  array 
pattern  in  the  plane  at  <p  =  45°  ,  where  the  maximum  beam  is 
directed  at  <p0  =  45°.  The  array  patterns  produced  by  this 
application  are  identical  to  the  patterns  given  in  [Ref.  3]. 
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FIGURE  4.2  Array  factor  pattern  of  a  5x5  element,  equally 
excited,  planar  array  with  dx  =  dy  =  A/2  m,  0O  =  30°, <p0  = 
45°,  in  the  x-z  and  y-z  planes 
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FIGURE  4.3  Array  factor  pattern  of  a  5x5  element,  equally 
excited,  planar  array  with  dx  =  dy  =  A/2  ™,  60  =  30°  ,  0O  = 
45°  ,  in  the  plane  <p   =  45° 
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V.  THE  CIRCULAR  ARRAY  ANTENNA 

Circular  arrays,  or  ring  arrays,  are  antennas  in  which 
the  radiating  elements  are  placed  on  a  circle  with  no  elements 
positioned  inside  the  circle.  Circular  antenna  arrays  with 
elements  positioned  inside  the  circle  are  planar  arrays,  and 
these  antennas  are  addressed  in  Chapter  IV.  Circular  array 
antenna  are  used  for  radio  direction  finding,  air  and  space 
navigation,  underground  propagation,  radar,  and  sonar  systems 
[Ref.  3:  p.  274].  As  with  linear  and  planar  arrays,  circular 
arrays  with  uniformly  spaced,  equally  excited  elements  are 
narrowband  antenna  systems  due  to  the  dependance  of  the  number 
of  elements  (N)  and  the  circumference  of  the  circle  to  the 
wavelength  (A)  [Ref.  3:  p.  277]. 

The  following  inputs  are  required  to  analyze  a  circular 
array: 

N  =  number  of  isotropic  radiating  elements 

a  =  radius  of  the  circle 

f  =  frequency  of  interest 

6„  =  direction  of  main  lobe  at  6  =  6 

o  o 

0O  =  direction  of  main  lobe  at  <p   =  (pQ 

IQ  =  antenna  feed  current 

Zo  =  characteristic  feed  impedance 
The  first  two  inputs  are  data  primarily  obtained  from  visual 
data  such  as  photographs,  and  the  *  indicates  input  parameters 
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that  are  either  known  or  estimated.  The  planar  array  geometry 
is  illustrated  in  Figure  5.1. 


FIGURE  5.1  Geometry  of  N-element  circular  array  of 
isotropic  point  sources 
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This  Mathcad  application  is  written  for  N-element 
circular  arrays  with  uniformly  spaced,  equally  excited, 
isotropic  radiating  elements  on  the  x-y  plane.  The 
polarization  of  the  array  is  not  calculated  since  isotropic 
point  sources  are  used,  and  the  antenna  input  resistance  (Rin) 
is  assumed  to  be  equivalent  to  the  radiation  resistance  (Rr)  . 
This  assumption  results  in  the  computation  of  ideal  gain  (G) . 
As  with  the  linear  and  planar  arrays,  the  radiation  pattern  of 
a  specific  element  type  is  neglected  in  the  analysis  of  the 
circular  array  since  in  normal  usage  it  will  have  little 
effect  for  a  circular  array  with  a  large  number  of  elements. 
However,  as  with  both  the  linear  and  planar  arrays,  the  user 
of  this  application  should  be  able  to  use  information 
regarding  the  type  of  element  and  the  orientation  of  the 
elements  to  narrow  the  range  of  possible  60's  and  <P0's. 

The  normalized  array  factor  (AF)  is  calculated  for  far- 
field  observations  neglecting  mutual  coupling  between  adjacent 
elements  with  equal  amplitude  excitation  and  is  [Ref.  3:  pp. 
274-278;  Ref.  4:  350-354]: 

AF  (6,4>)  =  ±Y   eJl*asineco8  «,-♦„)  *  «J  (dimensionless)    (5.1) 


where 

Qn   =  2ic/-j2\,  n   =  1,  2,  .  .  .  N        (radians)  (5.2) 


<x.n   =  -icasin0ocos(<|)o-$;i)    (radians)  (5.3) 
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Equation  (5.1)  is  the  normalized  array  factor.  Equation  (5.2) 
is  the  anqular  position  of  the  nth  element  on  the  x-y  plane. 
Equation  (5.3)  is  the  phase  excitation  of  the  nth  element  with 
the  position  of  the  main  beam  directed  at  0  =  0Q  and  0  =  <po, 
and  k  is  the  wavenumber  (2n/X)  and  a  is  the  radius  of  the 
array. 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

r   ;>  1 .6 1         (meters)  (5.4) 

r  z  5  (2a)         (meters)  (5.5) 


r  z 


2  (2,a)2    (meters)  (5.6) 


The  minimum  distance  to  the  far-field  is  found  by  taking  the 
maximum  of  (5.4),  (5.4),  and  (5.6),  where  the  diameter  of  the 
circle  is  taken  to  be  the  maximum  dimension  of  the  antenna. 

To  determine  the  directivity  (DQ)  for  a  circular  array, 
the  approximate  equations  of  radiation  intensity  (U(8,0))  for 
arrays  and  radiated  power  (Prad)  are  used.  Radiation  intensity 
is  [Ref.  3:  pp.  229-233]: 

[7(0,4))  =  (AF)2         (W/  solid  angle)  (5.7) 

The  radiated  power  is  [Ref.  3:  pp.  28]: 

Prad    =  J2*/*72^6'^  Sin(6)  <*<**     {W)  (5>8) 
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In  (5.8)  ,  Prad  is  computed  by  integrating  over  all  angles  above 
the  x-y  plane  (0  <  0  <  n/2)  [Ref.  7:  p.  78].  This  assumes 
that  the  array  has  a  pattern  maxima  above  the  x-y  plane  with 
no  radiation  below  the  x-y  plane.  The  directivity  for  the 
circular  array  is  [Ref.  3:  pp.  229-233]: 

DQ  =   —    (dimensionless)  (5.9) 

^rad 

In  (5.9),  Umax  is  unity  and  occurs  at  0  =  0O  and  0  =  0O  since 
the  array  factor  is  normalized. 

Effective  isotropic  radiated  power  (EIRP)  is  determined 
by  [Ref.  4:  p.  62]: 

EIRP  =   PradD0  (W)  (5.10) 

Assuming  that  the  power  delivered  to  a  circular  array  of 
N  isotropic,  uniformly  spaced  radiating  elements  is  eguivalent 
to  the  radiated  power,  we  can  estimate  then  the  radiation 
resistance  (Rr)  of  the  array  as  [Ref.  3:  p.  55]: 

?  (  P      ) 

Rr    =  p       (Q)  (5'11) 

Since  the  radiating  elements  are  assumed  to  be  ideal  point 
sources,  the  input  resistance  of  a  circular  array  (Rin)  is 
estimated  to  be  equivalent  to  the  radiation  resistance: 

Rin  =  Rr         (Q)  (5.12) 
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Since  the  input  resistance  is  estimated,  the  gain  (G)  of 
a  circular  array  cannot  be  calculated,  but  an  ideal  value  can 
be  derived.  As  with  linear  and  planar  arrays,  addressed  in 
Chapter  III  and  Chapter  IV,  respectively,  the  total  antenna 
efficiency  (et)  is  [Ref.  3:  pp.  44-45]: 

ec  =  tztcd   =  1-|r|2    (dimensionless)  (5.13) 


where 


e  d   =  1    (dimensionless)  (5.14) 


er  =  l-|r|2    (dimensionless)  (5.15) 


T  =  -^ ?    (dimensionless)  (5.16) 


In  (5.14),  ecd  is  the  conduction  and  dielectric  efficiency, 
which  is  unity  for  an  ideal  lossless  antenna.  In  (5.16),  Z0 
is  the  characteristic  impedance  of  a  transmission  line,  and  r 
is  the  voltage  reflection  coefficient.  Therefore,  the  total 
antenna  efficiency  is  equal  to  the  mismatch  efficiency  of 
(5.15)  when  an  ideal  lossless  antenna  is  assumed,  and  the  gain 
of  a  circular  array  is  [Ref.  3:  pp.  43-44]: 

G  =  ttD0        (dimensionless)  (5.17) 

G(dB)    =  10log10  (ttDQ)  (dB)  (5.18) 
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The  maximum  effective  aperture  (Aem)  of  a  circular  array 
is  derived  assuming  polarization  mismatches  are  negligible 
(this  Mathcad  application  is  constructed  using  isotropic 
radiating  elements  for  which  polarization  loss  factor  (PLF) 
and  antenna  polarization  are  not  calculated)  [Ref.  3:  p.  63]: 

Km   =  Gi})2  PLF         im2)  (5.19) 

PLF  =   1    (dimensionless)  (5.20) 

The  effective  height  (hem)  of  a  circular  array  is 
estimated  by  [Ref.  6:  p.  42]: 


Km   =  2 


^Z  {m)  (5.21) 

N    ^ 


In  (5.21),  rj0    is  the  characteristic  impedance  of  free  space 

(1207T)  . 

Four  two-dimensional  array  pattern  displays  are  produced 
in  this  Mathcad  application:  the  first  is  the  x-z  plane;  the 
second  is  the  y-z  plane;  the  third  is  a  plane  perpendicular  to 
the  x-y  plane  when  0  =  7r/4;  the  last  is  a  plane  perpendicular 
to  the  x-y  plane  when  0  =  <p0.  The  last  display  is  redundant 
when  0  is  either  0,  7r/4,  or  tt/2. 

The  directivity  of  a  uniform  circular  array  approaches 
the  number  of  elements  (N)  as  the  radius  (a)  of  the  array 
becomes  very  large  [Ref.  3:  p.  277].   Table  5.1  provides  the 
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directivity  of  a  uniform  circular  array  as  a  function  of  the 
radius  of  the  array. 

Table  5.1  Circular  Array  Directivity  Data 


RADIUS  (m) 

NUMBER  OF 
ELEMENTS 

DIRECTIVITY  (DJ 

1.59 

10 

13.70  dB 

10 

10 

12.83  dB 

20 

10 

12.89  dB 

50 

10 

11.14  dB 

100 

10 

10.33  dB 

Figures  5.2  and  5.3  provide  comparisons  of  the  array 
radiation  patterns  computed  for  a  10-element,  uniformly 
spaced,  equally  excited  circular  array  with  a  radius  of  1.59 
meters  with  those  in  [Ref.  3].  Figure  5.2  is  the  radiation 
pattern  in  the  x-z  plane,  and  the  radiation  pattern  in  the  y-z 
plane  is  shown  in  Figure  5.3.  For  these  patterns,  0O  =  0O  = 
0°.  The  array  radiation  patterns  generated  by  this 
application  are  identical  to  the  patterns  given  in  [Ref.  3]. 
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FIGURE  5.2  Array  factor  pattern  of  a  10-element,  uniformly 
spaced,  equally  excited  circular  array  with  a  =  1.59  m, 
and  60=00=0°,  in  the  x-z  plane 


38 


-1   -0,8  -0.6  -0.4  -0.2   0   0.2   0.4   0.6   0.8   1 


FIGURE  5.3  Array  factor  pattern  of  a  10-element,  uniform 
spaced,  equally  excited  circular  array  with  a  =  1.59  m, 
and  60  =  <p0   =  0°,  in  the  y-z  plane 
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VI.  THE  FOLDED  PI POLE  ANTENNA 

The  folded  dipole  is  a  center-fed  wire  antenna  that 
performs  in  a  manner  similar  to  that  of  a  linear,  center-fed 
thin-wire  antenna  with  the  added  advantage  of  a  more  easily 
matched  input  impedance  [Ref.  8:  p.  7-36].  The  folded  dipole 
antenna  consists  of  two  parallel  dipoles  connected  at  the  ends 
to  form  a  narrow  wire  loop.  The  folded  dipole  antenna  is 
illustrated  in  Figure  6.1. 


FIGURE   6.1  Geometry  of   folded  dipole 
antenna 
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The  folded  dipole  can  be  assumed  to  behave  as  a  thin-wire 
dipole  when  the  following  are  satisfied:  the  length  of  the 
antenna  is  much  greater  than  the  diameter  of  the  feed 
conductor  (L  >  2a)  ;  the  distance  between  the  conductors  is 
greater  than  the  radius  of  the  feed  conductor  (d  >  a)  ; 
(2nd/k)2  <  1;  all  conductors  are  assumed  to  be  perfect  [Ref. 
9:   p.  172]  ; 

This  Mathcad  application  calculates  the  parameters  of  a 
folded  dipole  by  assuming  the  antenna  behaves  as  a  thin-wire 
dipole  with  an  equivalent  radius  (ae)  .  The  current  maximum  is 
normalized  to  one  amp.  The  equivalent  radius  is  [Ref.  8:  p. 
7-38] : 

ra2ln(a)  +  b2ln(b)   +  2abln(d) 


a„  *  exp 


(a+b)2 


(meters)       (6.1) 


This  application  provides  the  characteristics  of  the  antenna 
both  in  free  space  and  positioned  horizontally  over  the  earth. 
Although  the  folded  dipole  positioned  vertically  over  the 
earth  occurs  frequently  in  VHF/UHF  communications,  it  is  not 
addressed  in  this  research. 

The  following  inputs  are  required  to  analysis  a  folded 
dipole: 

a  =  radius  of  feed  conductor 

b  =  radius  of  second  conductor 

d  =  spacing  between  conductors 

L  =  length  of  antenna 
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h  =  height  of  antenna  above  ground  plane 

f  =  frequency  of  interest 

ecd  =  conduction/dielectric  efficiency  of  conductors 

er  =  relative  dielectric  constant  of  ground  plane 

a   =   conductivity  of  ground  plane 

rff  =  distance  for  far-field  parameter  calculations 

aw  =  incoming  wave  electric  field  unit  vector  for  antenna 
in  free  space 

aa  =  unit  polarization  vector  for  antenna  in  free  space* 

ahw  =  incoming  wave  electric  field  unit  vector  for 
horizontally  positioned  antenna 

6hp  =  antenna  polarization  direction  for  horizontally 
positioned  antenna 
The  first  five  inputs  are  physical  dimensions  obtained  through 
photographs,  the  '  indicates  input  parameters  that  are  either 
known  or  estimated. 

A.  FOLDED  DIPOLE  IN  FREE  SPACE 

Assuming  the  folded  dipole  behaves  as  a  thin-wire  dipole 
with  ae  <  d  <  L,  we  have  the  far-field  electric  field 
intensity  of  the  antenna  aligned  with  the  z-axis  as  [Ref.  3: 
p.  120]: 


£"(9)  =JT)o^f 


-Jkrff 


2nrff 


|  —  cos(8))  -  cos/  — 
sin(8) 


/     /    v       (6.2) 

( v/m)      y        ' 


In  (6.2),  T70  is  the  intrinsic  impedance  of  free  space,  equal 
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to  1207T,  and  k  =  2ir/\  is  the  wavenumber.  The  normalized 
electric  field  is  found  by  dividing  (6.2)  by  the  maximum  value 
of  the  electric  field  which  occurs  at  0  =  ir/2  for  L  less  than 
about  1.25A..  Clearly,  the  radiation  pattern  of  the  folded 
dipole  is  identical  to  that  of  a  conventional  dipole  of  length 
L. 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

i  z   1.6  A.    (meters)  (6.3) 

r  ^  51    (meters)  (6.4) 


x   ^  -2Lf    (meters)  (6.5) 


The  minimum  distance  to  the  far-field  is  found  by  comparing 
the  values  of  (6.3),  (6.4),  and  (6.5)  and  selecting  the 
maximum  value.  The  distance  between  the  conductors  is  not 
taken  into  account  in  determining  the  minimum  distance  to  the 
far-field  since  it  is  assumed  that  the  folded  dipole  has  an 
over  all  effective  radius  (ae)  . 

Directivity  (D0)  for  the  folded  dipole  is  found  by 
determining  the  radiation  intensity  (U(0))  and  radiated  power 
(Prad)  .  Since  the  radiation  pattern  is  symmetric  in  the  x-y 
plane,  U(6)  and  Prad  are  functions  of  0  only.  Radiation 
intensity  is  [Ref.  3:  p.  28]: 
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r  2 
(7(0)  =_I^|e(6)|2    (W/ solid  angle)  (6.6) 

2r\o 


The  radiated  power  is  [Ref.  3:  pp.  28]: 

Prad  =  2nfnU(d)  siniB)  dQ         (W)  (6.7) 

Jo 

Directivity  for  a  folded  dipole  is  [Ref.  3:  pp.  29]: 

4  7t  t/ 

£>o  =  —         idimensionless)  (6.8) 

where   Umax   =   U(tt/2)  . 

Effective  isotropic  radiated  power  (EIRP)  is  determined 
by  [Ref.  4:  p.  62] : 

EIRP  =   PradDQ  (W)  (6.9) 

The  folded  dipole  is  unique  in  its  capability  to  act  as 
a  step  up  impedance  transformer  when  the  antenna  is  resonant. 
The  antenna  effectively  operates  as  an  unbalanced  transmission 
line,  and  it  can  be  analyzed  by  decomposing  the  current  into 
two  distinct  modes:  a  transmission  line  mode  [Figure  6.2]  and 
an  antenna  mode  [Figure  6.3].  [Ref.  3:  p.  341;  Ref.  4:  pp. 
205-207] 
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FIGURE   6.2   Transmission 
line  mode 


FIGURE  6 . 3  Antenna  mode 


The  input  impedance  (Zjn)  at  the  feed  for  the  folded 

dipole  is  given  in  (6.10),  and  the  radiation  resistance  (Rin) 

for  the  antenna  is  the  real  component  of  (Zin)  )   [Ref.  3: 
p. 342]  : 


«£c£d 

2Zd   +  Zc 


(Q) 


(6.10) 
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In  (6.10),  Zt  is  the  input  impedance  of  the  transmission  line 
mode  and  is  obtained  from  the  impedance  transfer  equation 

kL 


Zt  -   jz0tan|^) 


(CI) 


(6.11) 


where  Zo  in  (6.11)  is  the  characteristic  impedance  of  a  two- 
wire  transmission  line  and  is  given  in  (6.12)  [Ref.  3:  pp. 
341-342;  Ref.  8:  p.  7-39]: 


Z0   =  ^cosh"1 

71 


\ 


v  2/al5, 


(0) 


(6.12) 


The  input  impedance  of  a  linear  dipole  of  length  L  and 
equivalent  radius  ae  is  calculated  by  [Ref.  3:  pp.  124,  127, 
294,    and    342] : 

Zd  =  Rd  +  jXd         (CI)  (6.13) 


*d  = 


»1#) 


SI 


(CI 


(6.14) 


Xd  = 


X, 


sin' 


kL 


(CI) 


(6.15) 


Rr  = 


*\o 


j^[y  +  ln(icL)    -  Cd(kD] 


2n 


nc 


2tc 


-|sin(icL)  [Si(2kD    -  2Si(kL)} 


±cos(kL)ty  +  ln|— )  +  Ci(2kL)-2Ci(kL)\ 


(6.16) 
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xr   = 


^■[2Si(icL)  +cos(kL)  [2Si(kL)-Si(2kD]] 


(6.17) 


n0 


4  TC 


/ 


sin ( kL) 


2Ci(kL)  -Ci(2kL)  -Ci 


2ka( 


In  (6.16)  and  (6.17),  Ct(x)  and  S^x)  are  the  cosine  and  sine 
integrals,  respectively,  [Ref.  3:  pp.  124  and  743-746]: 


x cos (t) 


:.{x)   .Jj*  cob  it  j  dx      Y+ln(x)_£ 


n  =  l 


(  -l)n-xx2n 
2n (2n) ! 


(6.18) 


sd(x)  =  r 

Jo 


x  sin(x) 


*-  e 


(  -1  )  n_1  x  *2n_1' 
^  (2/2-1)  (2/2-1)  ! 


(6.19) 


and  y  -  0.57721  (Euler's  constant). 

Mathcad  computes  summations  more  quickly  than  integrals. 
In  this  application,  the  summation  forms  of  (6.18)  and  (6.19) 
are  used  with  the  index  n  ranging  from  1  to  50  to  ensure 
proper  convergence.  The  maximum  value  of  the  argument  of 
either  C^x)  or  S,(x)  is  x  =  47TL/A.,  and  for  L  =  A  the 
succeeding  terms  in  the  summation  rapidly  decrease  for  n  >  20. 

The  gain  (G)  of  a  folded  dipole  is  obtained  as  the 
product  of  the  antenna  reflection  efficiency  (et)  and  the 
directivity  of  (6.8)  where  the  antenna  reflection  efficiency 
is  [Ref.  3:  pp.  43-45] : 

zt  =  trecd  =   ecd(l-|r|2)    idimensionless)  (6.20) 
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where 

er  =  l-|T|2         (dimensionless)  (6.21) 


R     —  Z 

T  =  —is £    (dimensionless)  (6.22) 


In  (6.21)  er  is  the  mismatch  efficiency  and  ecd  is  the 
conduction  and  dielectric  efficiency  in  (6.22).  In  (6.22),  ZQ 
is  the  characteristic  impedance  of  the  transmission  line  given 
in  (6.12).   Therefore,  the  gain  of  the  antenna  is: 

G  =  ecD0        {dimensionless)  (6.23) 

G(dB)    =   10log10  (ttD0)  (dB)  (6.24) 

The  folded  dipole  in  free  space  aligned  on  the  z-axis  is 
horizontally  polarized  with  respect  to  the  x-y  plane. 
Therefore,  given  the  direction  of  the  incoming  wave's  unit 
polarization  vector  (a  )  ,  the  polarization  loss  factor  (PLF) 
is  [Ref.  3:  pp.  48-53]: 

PLF  =  \aw-  o*|2    (dimensionless)  (6.25) 

The  maximum  effective  aperture  (A^)  of  a  folded  dipole 
is  [Ref.  3:  p. 63] : 

z   =  GiX)2  PLF         (m2)  (6.26) 

ew  All 
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The  effective  height  (hem)  of  the  folded  dipole  is  [Ref. 
6:  p.  42]  : 


hem  =  2\ 


R'A™  (m)  (6.27) 


The  bandwidth  of  the  folded  dipole  antenna  is  dependant 
on  the  length  of  the  wires  and  has  better  bandwidth 
characteristics  than  a  conventional  dipole  of  the  same  length. 
It  can  be  assumed  that  the  bandwidth  of  the  folded  dipole  is 
essentially  the  same  as  that  of  a  conventional  dipole  of  an 
radius  ae  (a  <  ae  <  d/2)  [Ref.  3:  p.  346].  The  bandwidth  is 
increased  when  the  distance  between  the  conductors  (d)  is 
increased  and/or  when  the  conductor  radii  a  and  b  are 
increased  [Ref.  10:  p.  190]. 

Tables  6.1  and  6.2  are  comparisons  of  measured  and 
calculated  input  impedance  data  for  the  folded  dipole  as  a 
function  of  its  electrical  length.  The  radii  of  the 
conductors  (a,  b)  are  the  same  (0.0005  m)  while  the  spacing 
between  the  conductors  (d)  is  held  constant  at  0.00625  m  and 
X    =    1  m  [Ref.  8:  pp.  7-40  -  7-41]. 
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TABLE  6.1  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 
LENGTH  ( L) 

MEASURED 
Re(Zin)  n 

CALCULATED 
Re(Zin)  n 

k/ir 

2000 

6533  ** 

2A/7T 

2400 

2592 

3A/7T 

0.0 

0.1 

AX/ti 

300 

324 

5A/7T 

1700 

1685 

6A/7T 

0.5 

1.0 

7A/7T 

100 

100 
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TABLE  6.2  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 
LENGTH  (L) 

MEASURED 
Im(Zin)  n 

CALCULATED 

lm(Zin)  n 

k/TT 

5000 

5437 

2X/7T 

-650 

-1204  ** 

3X/TT 

-90 

-85 

4A/7T 

1000 

1061 

5A/7T 

700 

793 

6X/7T 

-190 

-174 

7A/7T 

600 

633 

To  further  test  the  validity  of  the  Mathcad  application, 
a  numerical  antenna  analysis  program  was  also  used.  The 
program  used  was  ELNEC,  a  powerful,  easy-to-use  program  for 
modeling  and  analyzing  a  wide  variety  of  antenna  types 
including  ground  effects  and  parasitic  structures.  The 
fundamental  computation  portion  of  ELNEC  is  the  same  as 
MININEC  (Version  3)  ,  that  was  developed  by  the  Naval  Ocean 
Systems  Center  [Ref.  11:  pp.  5-7]. 

Table  6.3  is  a  comparison  of  the  gain  for  a  folded  dipole 
in  free  space  with  different  dimensions  [Ref.  8:  p.  3-29]. 
For  each  case,  the  radii  of  the  conductors  (0.0005  m)  and  the 
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distance  between  the  conductors  (12.5a  m)  are  the  same.  As 
can  be  seen,  the  comparison  between  the  gains  obtained  from 
various  sources  is  very  good. 

TABLE  6.3  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 
DIMENSION 

MEASURED 
Gain  (G) 

ELNEC 
Gain  (G) 

Mathcad 
Gain  (G) 

L  =  A/4  m 
d  =  L/13  m 

1.71  dB 

1.70  dB 

1.64  dB 

L  =  X/2    m 
d  =  L/2  5.5  m 

2.14  dB 

2.20  dB 

2.15  dB 

The  computed  parameters  in  Tables  6.1,  6.2,  and  6.3  are 
essentially  the  same  as  those  measured  or  computed  with  ELNEC 
with  exception  of  the  ones  marked  with  **.  These  discrepancies 
result  in  using  the  transmission  line  model  and  equivalent 
radius  to  analysis  the  folded  dipole  in  this  application. 

As  previously  mentioned,  the  folded  dipole  is  a  center- 
fed  wire  antenna  that  performs  in  a  manner  similar  to  a 
linear,  center-fed,  thin-wire  antenna.  Therefore,  the 
radiation  pattern  of  a  folded  dipole  is  also  similar  to  that 
of  a  thin-wire  antenna.  The  folded  dipole  radiation  pattern 
in  free  space  is  illustrated  in  Figure  6.4. 
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FIGURE  6.4  Radiation  Pattern  folded  dipole  in  free  space 
with  L  =  A/2  m,  a  =  b  =  0.0005  m,  and  ae  =  0.004  m 


Figure  6.5  is  a  comparison  of  the  radiation  pattern  of  a 
folded  dipole  in  free  space  obtained  using  ELNEC  and  the 
Mathcad  application.  The  patterns  are  not  identical,  but  are 
similar  and  oriented  in  the  same  direction  as  expected.  The 
difference  can  be  attributed  to  the  fact  that  ELNEC  takes  into 
account  mutual  coupling  between  the  conductors  while  the 
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Mathcad  application  models  the  folded  dipole  as  a  thin-wire, 
center-fed  dipole. 


-1       -0.8     -0.6     -0,4     -0.2       0        0.2       0.4       0.6       0.8        1 


-  Mathcad  Pattern,   ...  ELNEC  Pattern 


Figure  6.5  Comparison  of  radiation  patterns  using  ELNEC 
and  Mathcad  of  a  Folded  dipole  in  free  space  with  L  =  k/2 
m,  a  =  b  =  0.0005  m,  and  ae  =  0.004  m 
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B.  FOLDED  DIPOLE  POSITIONED  HORIZONTALLY  OVER  THE  EARTH 

Through  the  use  of  image  theory,  the  far-field  electric 
field  intensity  of  a  folded  dipole  oriented  horizontally  above 
a  flat  earth  can  be  obtained.  The  0  component  of  the  electric 
field  of  a  horizontal  dipole  (parallel  to  the  y-axis)  in  free 
space  is  [Ref.  3:  pp.  143-144]: 

kL 


2?e<e,4>)  =jti 


-jkr. 


2tzi 


ff 


cos 


sin  (6)  sin(<|>)\  -cos/— ) 
v/l-sin2(6)  sin2(<}>) 


iV/m) 


(6.28) 


From  image  theory,  the  6  component  of  the  electric  field 
over  ground  is  [Ref.  4:  pp.  229-235]: 

£6^(6)  =  £0(0,(t>)(e^lcos(e>  -  rve-^cos(e>)    (V/m)  (6.29) 


In  (6.28)  and  (6.29),  Ee(0,0)  is  valid  only  above  the  ground 
(0  <  0  <  7T/2)  [Ref.  3:  pp.  135-142].  In  (6.29),  h  is  the 
distance  from  the  earth  to  the  antenna  and  the  vertical 
reflection  coefficient  (Tv)  of  the  ground  is  [Ref.  4:  pp.  231- 
233]  : 

e'rcos(0)  -  yjz' z-   sin2  (0) 


r„  = 


e/rcos(0)  +  yjt' z-   sin2  (0) 


(dimensionless)     (6.30) 


where 


e      =  e    -   7 

r         br       J 


2  71  ft. 


(dimensionless) 


(6.31) 
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In  (6.30)  and  (6.31),  (e'r)  is  the  relative  complex  effective 
dielectric  constant  of  the  ground,  (er)  is  the  relative 
dielectric  constant  of  the  ground,  and  (a)  is  the  conductivity 
of  the  ground.  For  a  perfect  ground  plane,  a  -*■  oo  and  rv  =  1. 
The  $  component  of  the  electric  field  of  a  horizontal 
dipole  (parallel  to  the  y-axis)  in  free  space  is  [Ref.  3:  pp. 
135-142] : 


I  e'jkZff 
F<J>(8,4>)  =jx\0 


2nrff 


cos/  —  sin(8)sin(4>))-cos(— ) 
yjl  -sin2  (6)  sin2(<t>) 


..„    .   (6.32) 
( V/m) 


Now  from  image  theory,  the  §  component  of  the  electric  field 
over  earth  is  [Ref.  4:  pp.  229-235]: 

£4)^(6)  =  E<$(Q,<b)(ejkhcos{e)  +  rhe-^cos(8))    (V/m)         (6.33) 

In  (6.33),  (Th)  is  the  horizontal  reflection  coefficient  of 
the  ground  and  is  [Ref.  4:  pp.  231-232]: 

cos(6)  -  Jt'   -   sin2  (6) 
rh  =  v  (dimensionless)         (6.34) 

cos(0)  +  yjt'z-   sin2  (6) 

When  a  -+  oo,  rh  =  -1. 

The  radiation  intensity  (Uh(6,0))  of  a  folded  dipole  over 
a  earth  must  take  into  account  both  the  6  and  §  components  of 
the  electric  field.  Therefore,  the  radiation  intensity  is 
[Ref.  3:  p.  28]: 
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r    2 

1  ff 


uh^  =^££-[l£6pP(0)  l2  +  |£*sp(6)  I2]     (W/ solid  angle)      (6.35) 


2T1 


The  radiated  power  (Phrad)  follows  (6.7)  except  that  the  field 
is  present  only  above  the  earth.   Thus, 

Phrad  =  2n[  2  £7(8)  sin(0)  dQ         (W)  (6.36) 

J  o 

Directivity  (Dh0)  and  effective  isotropic  radiated  power 
(EIRPh)  for  a  horizontal  folded  dipole  over  the  earth  are 
identical  to  (6.8)  and  (6.9),  respectively,  with  the  results 
of  (6.36)  inserted. 

Input  resistance  (Rhin)  for  the  folded  dipole  over  the 

earth  plane  is  found  by  (where  I0  in  (6.32)  is  normalized  to 

one  amp)  [Ref.  3:  p.  124]: 

Rh 
**»•        .     /kL\  <Q)  (6-37) 

where  the  radiation  resistance  (Rhr)  is, 

2  Ph 
Rhr   =   IS*  (CI)  (6.38) 

The  voltage  reflection  coefficient  (TH) ,  gain  (Gh)  , 
maximum  effective  aperture  (Ahem)  ,  and  maximum  effective  height 
(hhem)  for  a  horizontal  folded  dipole  over  the  earth  are 
determined  using  the  equations  for  those  parameters  for  a 
folded  dipole  in  free  space,  as  listed  above. 
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The  Polarization  Loss  Factor  (PLFh)  for  the  antenna  over 
the  earth  can  be  determined  when  the  direction  of  the  main 
beam  of  the  polarized  antenna  (8hp)  is  estimated.  The  antenna 
polarization  vector  (aha)  is  then  calculated,  and  estimating 
the  unit  polarization  vector  of  the  incoming  wave  (o\x„)  ,  the 
PLF  is  found  from  (6.25). 

Table  6.4  is  a  comparison  of  measured  [Ref.  8:  p.  3-29], 
ELNEC,  and  Mathcad  software  gains  for  a  horizontal  folded 
dipole  over  a  reflecting  sheet.  The  radii  of  the  conductors 
(0.0067  m)  are  the  same.  As  can  be  seen,  the  Mathcad  software 
gain  is  similar  to  the  measured  and  ELNEC  gains. 

TABLE  6.4  Folded  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

ELNEC 

Mathcad 

DIMENSION 

Gain  (G) 

Gain  (G) 

Gain  (G) 

L  =  A/2  m 

d  =  L/2  5.5  m 

7.14  dB 

6.98  dB 

7.19  dB 

h  =  A/8  m 

No  input  impedance  data  has  been  found  to  compare  a 
folded  dipole  oriented  horizontally  over  the  earth  with  that 
computed  by  this  Mathcad  application. 

Figures  6.6  and  6.7  are  the  H-plane  and  E-plane  radiation 
patterns  of  a  folded  dipole  oriented  horizontally  over  the 
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earth,  respectively.  Since  the  folded  dipole  is  similar  to  a 
thin-wire  center-fed  dipole  of  length  L,  the  illustrations 
below  are  similar  to  those  of  a  dipole  of  length  L. 


O.B  - 


0.6  - 


0.4-  - 


0.2 


•1     -0.8    -0.6    -0.4    -0.2 


0.2     0.4     0.6     0.8 


FIGURE  6.6  H-plane  radiation  pattern  of  a  folded  dipole 
with  L  =  X/2  m,  a  =  b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5 
m,  a   =  0.01  S/m,  er  =  15 
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FIGURE  6.7  E-plane  radiation  pattern  of  a  folded  dipole 
with  L  =  A/2  m,  a  =  b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5 
m,  a   =  0.01  S/m,  er  =  15 


Figures  6.8  and  6.9  are  comparisons  of  the  H-plane  and  E- 
plane  radiation  patterns  between  those  obtained  with  ELNEC  and 
Mathcad  software  of  a  folded  dipole  oriented  horizontally  over 
the  earth,  respectively.  As  expected,  the  patterns  are 
similar.  As  previosly  mentioned,  the  difference  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
coupling  between  the  conductors  while  the  Mathcad  application 
is  modeled  after  a  thin-wire,  center-fed  dipole. 
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Figure  6.8  Comparison  between  ELNEC  and  Mathcad  of  H-plane 
radiation  patterns  of  a  folded  dipole  with  L  =  A/2  m,  a  = 
b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5  m,  a  =  0.01  S/m,  er  = 
15 
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Figure  6.9  Comparison  between  ELNEC  and  Mathcad  of  E-plane 
radiation  patterns  of  a  folded  dipole  with  L  =  A/2  m,  a  = 
b  =  0.0005  m,  d  =  0.00625  m,  h  =  0.5  m,  o  =  0.01  S/m,  er  = 
15 
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VII.  THE  CAGED  PI POLE  ANTENNA 

The  caged  dipole  antenna  is  used  in  communications  for 
its  broadband  antenna  characteristics.  Its  performance  is 
similar  to  that  of  a  thick  cylindrical  dipole.  Since  thick 
cylindrical  dipoles  are  usually  too  clumsy  to  use  except  at 
very  short  wavelengths,  thick  dipoles  are  replaced  with  a 
caged  dipole  which  consists  of  a  number  of  thin  wire  dipoles 
arranged  to  form  a  circular  cross  section.  The  antenna  is 
usually  supported  by  a  central  conductor  and  support  members 
that  provide  rigidity  but  contribute  little  to  the  electric 
properties  of  the  caged  dipole.  The  caged  dipole  (without  the 
central  conductor  and  the  support  members)  is  illustrated  in 
Figure  7.1  [Ref.  12:  pp.  172-174]. 


FIGURE  7.1  Geometry  of  caged  dipole  antenna 
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This  Mathcad  application  models  the  caged  dipole  antenna 
as  a  circular  array  of  center-fed,  equally  excited,  uniformly 
spaced,  thin-wire  dipoles.  This  approach  is  taken  to  account 
for  the  number  of  conductors  that  are  present  in  the  caged 
dipole.  As  with  the  arrays  previously  considered,  mutual 
coupling  between  the  dipoles  is  neglected.  In  this  manner, 
this  application  simulates  the  characteristics  of  a  thick 
cylindrical  dipole.  The  thin-wire  dipoles  are  modeled  after 
finite  length  dipoles  in  free  space  that  have  their  length 
much  greater  than  their  diameter  (L  >  wire  diameter) ,  and  the 
dipoles  are  assumed  to  be  perfect  conductors.  Since  the 
central  conductor  and  the  support  wires  have  little  effect  on 
the  electrical  properties  of  a  caged  dipole,  mutual  coupling 
from  these  components  are  not  taken  into  account  in  this 
application.  As  seen  in  Figure  7.1,  the  inner  ends  of  a  caged 
dipole,  the  feed  points,  maybe  coned  in  order  to  reduce  base 
capacity  [Ref.  12:  pp.  172-174].  This  effect  is  also  not 
addressed. 

This  application  provides  the  characteristics  of  the 
antenna  in  free  space  as  well  as  and  oriented  both  vertically 
and  horizontally  over  the  earth. 


A.  CAGED  DIPOLE  IN  FREE  SPACE 

The  following  inputs  are  required  to  analyze  a  caged 
dipole  in  free  space: 
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N  =  number  of  conductors 

a  =  radius  of  antenna 

L  =  length  of  antenna 

f  =  frequency  of  interest 

I0  =  antenna  feed  current 

Z0  =  input  impedance 

ecd  =  conduction/dielectric  efficiency  of  conductors 

rff  =  distance  for  far-field  parameter  calculations 

aw  =  incoming  wave  electric  field  unit  vector  for 
antenna 

aa  =  unit  polarization  vector  for  antenna 
The  first  three  inputs  are  physical  dimensions  obtained 
primarily   from   photographs,   and   the   *   indicates   input 
parameters  that  are  either  known  or  estimated. 

Assuming  the  caged  dipole  behaves  as  a  circular  array 
with  equally  excited  (In)  ,  uniformly  spaced  elements  in  the  x- 
y  plane  with  a  radius  of  (a)  ,  we  have  the  array  factor 
(AF(8,0))  [Ref.  3:  pp.  274-278]: 

AF(8,4>)  =  £  j^ju-sineco.^-^)]    (dimensionless)      (7.1) 

27  =  1 

where  the  angular  position  in  the  nth  element  in  the  x-y  plane 
is: 

&n  =   2tW  —  j,  n   =  1,  2,  .  .  .  N        (radians)  (7.2) 
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Equation  (7.1)  is  derived  from: 

AF(6,4>)  =  £  Ja^[l:a.ineeo.(*-«.)+«J   {dimensionless)       (7.3) 


n  =  l 


where  the  phase  excitation  of  the  nth  element  is: 

an  =   0  =  -icasin60cos  ($0-$n)         (radians)  (7.4) 

In  (7.3),  and  the  phase  excitation  of  the  nth  element  (an)  is 
zero  as  the  result  of  desiring  the  main  beam  of  the  array  to 
be  directed  perpendicular  to  the  length  of  the  antenna  (60  = 
tt/2)  .   In  (7.1),  k  is  the  wavenumber  (2ir/X)  . 

Assuming  each  conductor  of  a  caged  dipole  behaves  as  a 
thin-wire  dipole,  we  have  the  far-field  electric  field 
intensity  of  a  single  conductor  aligned  with  the  z-axis  as 
[Ref.  3:  p.  120]: 


Bid)  =jr)0±f 


-okrff 


2nxff 


[  — cos(6)j  -  cosl— 
sin(6) 


,      /     v        (7.5) 

( V/m)      v        ' 


In  (7.5),  ?70  =  1207T  and  is  the  intrinsic  impedance  of  free 
space. 

As  previously  discussed  in  Chapter  III,  the  total 
electric  field  of  an  array  constructed  with  identical  elements 
is  the  product  of  the  field  of  a  single  element  and  the  array 
factor  as  a  consequence  of  the  principal  of  pattern 
multiplication   [Ref.  3:  p.  207].      Therefore,   the  total 
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electric  field  of  a  caged  dipole  in  free  space  is  approximated 
as: 

Et(d,$)  =  E(Q)AF(6,  $)    (V/m)  (7.6) 

An  observation  point  is  considered  to  be  in  the  far-field 
if  all  of  the  following  are  satisfied  [Ref.  4:  pp.  24-25]: 

r  ^  1.61    (meters)  (7.7) 

r  z  SsjL2  +    (2a)2         {meters)  (7.8) 


r  ±    2(L2   +    (2a)2)  {rneteis)  (7.9) 

A 


The  minimum  distance  to  the  far-field  is  found  by  taking  the 
maximum  of  (7.7)  through  (7.9).  In  (7.8)  and  (7.9),  the 
maximum  dimension  of  the  antenna  is  the  diagonal  length  of  the 
array. 

Directivity  (D0)  for  a  caged  dipole  is  found  by 
determining  the  radiation  intensity  (U(6,0))  and  radiated 
power  (Prad)  .   Radiation  intensity  is  [Ref.  3:  p.  28]: 

T       2 

£7(0,4))  =  _«_|£  (e,(|>)  I2    (W/  solid  angle)  (7.10) 

The  radiated  power  is  [Ref.  3:  pp.  28]: 

Prad    =     (2"[nU(e,4>)  Sin(0)   d8d<J)  (W)  (7.11) 

JO     Jo 
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Directivity  for  a  caged  dipole  is  [Ref.  3:  pp.  29]: 

4  71  If 
Do  =   —         (dimensionless)  (7.12) 

Prad 

where  Umax  =  U(tt/2,0)  . 

Effective  isotropic  radiated  power  (EIRP)  is  determined 
by  [Ref.  4:  p.  62]: 

EIRP  =   PiadDQ  (W)  (7.13) 

Input  resistance  (Rin)  for  a  caged  dipole  in  free  space 
is  found  from  [Ref.  3:  p.  124]: 

Rin   =    .  JkL\  (Q)  (7.14) 

sin2' 


where  the  radiation  resistance  (Rr)  is: 

2P 


rad 

-  n  I 


Rz   =    — -^    (CI)  (7.15) 

J." 


The  gain  (G)  of  a  caged  dipole  is  obtained  as  the  product 
of  the  antenna  reflection  efficiency  (et)  and  the  directivity, 
where  the  antenna  reflection  efficiency  is  ([Ref.  3:  pp.  43- 
45]  : 

ec  =  tztcd  =   ecd(1~|r'|2)    (dimensionless)  (7.16) 

where 

er  =  l-|r|2        (dimensionless)  (7.17) 
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r  =  -^5 °         (dimensionless)  (7.18) 


The  parameter  er  is  the  mismatch  efficiency,  ecd  is  the 
conduction  and  dielectric  efficiency  which  is  unity  for  a 
lossless  conductor,  and  Z0  is  the  characteristic  impedance  of 
a  transmission  line.  Therefore,  the  gain  of  the  antenna  is: 

G  =  ttD0        (dimensionless)  (7.19) 

G(dB)    =   10log10  (ttD0)  (dB)  (7.20) 

A  caged  dipole  in  free  space  aligned  parallel  to  the  z- 
axis  is  horizontally  polarized  with  respect  to  the  x-y  plane. 
Therefore,  given  the  direction  of  the  incoming  wave's  unit 
polarization  vector  (aw)  ,  the  polarization  loss  factor  (PLF) 
is  [Ref.  3:  pp.  48-53]: 

PLF  =  \ow-  a*a\2        (dimensionless)  (7.21) 

The  maximum  effective  aperture  (Aem)  of  a  caged  dipole  is 
[Ref.  3:  p. 63] : 

Km  =    GU)2  PLF         (m2)  (7.22) 

The  effective  height  (hem)  of  the  caged  dipole  is  [Ref. 
6:  p.  42]: 


Km    =    2 


\ 


RrA™  (m)  (7.23) 
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The  bandwidth  of  a  thick  cylindrical  dipole  antenna  is 
dependent  on  the  ratio  of  the  length  of  the  antenna  to  the 
diameter  of  the  antenna  (L/2a)  .  When  L/2a  «  5000  the 
acceptable  bandwidth  of  a  thick  cylindrical  dipole  is  about  3% 
of  the  center  frequency,  and  when  L/2a  «  260  the  bandwidth  is 
about  30%.  Hence,  if  we  assume  a  linear  relationship  between 
bandwidth  and  L/2a,  and  since  the  caged  dipole 's  performance 
is  similar  to  that  of  a  thick  cylindrical  dipole,  then  the 
bandwidth  of  a  caged  dipole  is  [Ref.  3:  p.  333]: 

BW   =  fhigh   -   flow  (Hz)  (7.24) 

where  the  upper  and  lower  frequencies  are: 

fM*»  =  f„  +  Fl4r\fn         (Hz)  (7.25) 


high         "-c  \  2a]     c 


1^) 


In  (7.25)  and  (7.25),  fc  is  the  center  frequency  and  F(L/2a) 
is  the  linear  relationship  function  between  antenna  length  and 
diameter.  This  relationship  is  obtained  by  solving  y  =  mx  + 
fc  for  m  using  x:  =  5000,  y:  =  3%  and  x2  =  260,  y2  =  30%.  The 
relationship  becomes: 

fI  —  \  =  -5.696-10'5  I  —  \+0.3148  (dimensionless)       (7.27) 


2a  \2a 
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Table  7.1  is  a  comparison  of  the  gain  of  a  thick 
cylindrical  dipole  with  that  of  a  caged  dipole  (N  =  8)  in  free 
space  computed  both  by  ELNEC  and  Mathcad  software  [Ref.  8:  p. 
2-38].  As  previously  mentioned  in  Chapter  VII,  ELNEC  is  an 
easy-to-use  numerical  analysis  program  for  modeling  and 
analyzing  different  antenna  types  [Ref.  11:  pp.  5-7].  As  can 
be  seen,  the  comparison  between  the  measured,  ELNEC,  and 
Mathcad  application  values  are  similar. 

TABLE  7.1  Caged  Dipole  Antenna  Data  Comparison 


ANTENNA 

MEASURED 

ELNEC 

Mathcad 

DIMENSIONS 

Gain  (G) 

Gain  (G) 

Gain  (G) 

L  =  X/2    m 

a  =  0. 025  m 

2.14  dB 

1.87  dB 

1.88  dB 

N  =  8 

Figure  7.2  is  the  radiation  pattern  of  a  caged  dipole  (N 
=8)  in  free  space.  It  is  essentially  eguivalent  to  that  of 
a  center-fed,  thick  cylindrical  dipole  having  the  same 
conductor  length  and  diameter  [Ref.  8:  p.  2-38]: 
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FIGURE  7.2  Radiation  pattern  of  a  caged  dipole  in  free 
space  with  L  =  X/2   m,    N  =  8,  and  a  =  0.125  m 


Figure  7.3  is  a  comparison  of  the  radiation  pattern  of  a 
caged  dipole  (N  =  8)  obtained  using  ELNEC  and  Mathcad.  The 
patterns  are  not  identical,  but  are  similar  and  oriented  in 
the  same  direction  as  expected.  The  difference  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
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coupling  between  elements  while  the  Mathcad  application 
ignores  these  effects. 


-1       -0.8     -0.6     -0.4     -0.2       0        0.2       0.4       0.6       0.8         1 


-  Mathcad  Pattern,   ...  ELNEC  Pattern 


FIGURE  7.3  Comparison  of  radiation  patterns  using  ELNEC 
and  Mathcad  of  a  caged  dipole  in  free  space  with  L  =  A./2 
m,  N  =  8,  and  a  =  0.125  m 
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B.  CAGED  DIPOLE  ORIENTED  VERTICALLY  OVER  EARTH 

The  following  known  or  estimated  inputs  are  additional 
data  required  to  analyze  a  caged  dipole  oriented  vertically- 
over  the  earth.  The  height  of  the  antenna  above  ground  is  a 
physical  dimension  obtained  through  photographs: 

h  =  antenna  height  above  ground 

er  =  relative  dielectric  constant  of  earth 

a  =  conductivity  of  earth 

avw  =  incoming  wave  electric  field  unit  vector  for 
vertical  antenna 

6vp  =  antenna  polarization  direction  for  vertically 
orientated  antenna 

As  with  the  caged  dipole  in  free  space,  the  analysis  of 
a  caged  dipole  positioned  vertically  above  the  ground  does  not 
take  into  account  the  central  conductor  and  support  members 
since  they  produce  negligible  effects  on  the  antenna's 
electrical  properties. 

From  image  theory  and  pattern  multiplication  for  arrays, 
the  total  far-field  electric  field  intensity  of  a  caged  dipole 
oriented  vertically  above  a  flat  earth  is  [Ref.  4:  pp.  229- 
235]  : 
Evt(6,$)  =Ev(B)  AFv(Q,<b)(eikhcosie)  +  Tve-jkhcosie))        (V/m)  (7.28) 

In  (7.28),  Ev(8)  is  the  far-field  electric  field  intensity  of 
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a  finite  length  dipole  in  free  space  which  is  given  by  (7.5), 
and  AFv(0,0)  is  the  array  factor  for  a  circular  array  which  is 
(7.1).  The  total  electric  field  intensity  is  valid  only  above 
the  ground  (0  <  6  <  it/2)  [Ref.  3:  pp.  135-142].  The  distance 
from  the  earth  to  the  center  of  the  antenna  is  (h) ,  and  the 
vertical  reflection  coefficient  (Tv)  of  the  ground  is  [Ref.  4: 
pp.  231-  233]  : 

„  n  e'rcos(Q)  -  Je'-sin2  (0) 

rv(6)  =  — - — v  r      —  (dimensionless)     (7.29) 

e'rcos(0)  +  ^e^-sin2  (6) 


where 


t'     =   e  -  j - —    (dimensionless)  (7.3  0) 

1         r    2n  ft0  v 


In  (7.29)  and  (7.30),  (e'r)  is  the  relative  complex  effective 
dielectric  constant  of  the  ground,   (er)   is  the  relative 
dielectric  constant  of  the  ground,  (a)  is  the  conductivity  of 
the  ground,  and  e0  is  the  permittivity  of  free  space  (l/(367r) 
x  10"9)  .   For  a  perfect  ground  plane,  a  ->   <»  and  rv  =  1 . 

The  radiation  intensity  (Uv(6,</>))  of  a  caged  dipole  over 
earth  is  [Ref.  3:  p.  28]: 

C/v(6,4>)  =£^[|£vc(6,(i>)  |2]    (w/ solid  angle)  (7.31) 

The  radiated  power  (Pvrad)  is  present  only  above  the  earth. 
Thus,  [Ref.  3:  p.  28] : 
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— 

Pvrad  =   f  *f  2Uv(6,<b)  sin(6)  dQdfy        (w)  (7.32) 

Jo  Jo 

Directivity  (Dv0)  and  effective  isotropic  radiated  power 
(EIRPv)  for  a  vertically  oriented  caged  dipole  over  the  earth 
are  given  by  (7.12)  and  (7.13),  respectively,  with  the  results 
of  (7.32)  inserted.  Uvmax  in  (7.12)  is  simply  the  maximum 
value  of  (7.31)  which,  in  Mathcad,  is  found  by  constructing 
(7.31)  as  a  vector  and  then  finding  the  maximum  value  of  the 
vector. 

Input  resistance  (Rvin)  ,  voltage  reflection  coefficient 
(TV)  ,  gain  (Gh)  ,  maximum  effective  aperture  (Avem)  ,  and  maximum 
effective  height  (hvem)  for  a  vertical  caged  dipole  over  the 
earth  are  determined  using  the  eguations  for  those  parameters 
for  a  caged  dipole  in  free  space. 

The  Polarization  Loss  Factor  (PLFv)  for  the  antenna  over 
the  earth  can  be  determined  when  the  direction  of  the  main 
beam  of  the  polarized  antenna  (0vp)  is  estimated.  The  antenna 
polarization  vector  (ava)  is  then  calculated  and,  estimating 
the  unit  polarization  vector  of  the  incoming  wave  (avw)  ,  PLF 
is  found  from  (7.21). 

Table  7.2  is  a  comparison  of  the  gain  of  a  vertical  caged 
dipole  over  a  flat  earth  obtained  with  both  ELNEC  and  the 
Mathcad  application.  As  can  be  seen,  the  gains  are  similar. 
As  previously  mentioned,  the  difference  in  the  gains  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
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coupling  between  the  elements  as  well  as  ground  effects.  The 
Mathcad  application  ignores  mutual  coupling  between  elements 
and  approximates  ground  effects. 

TABLE  7.2  Caged  Dipole  Antenna  Gain  Comparison 


ANTENNA 

ELNEC 

Mathcad 

DIMENSIONS 

Gain  (G) 

Gain  (G) 

L  =  A./2  m 

a  =  0. 125  m 

N  =  8 

4.25  dB 

4.50  dB 

h  =  0.25  m 

a    =    0.01  S/m 

er  =  15 

Since  the  caged  dipole  behaves  in  a  manner  similar  to 
that  of  a  thin  dipole  except  for  bandwidth  characteristics, 
the  radiation  pattern  of  a  vertically  oriented  caged  dipole 
over  ground  will  resemble  that  of  a  thin  vertical  dipole  over 
ground.  Figure  7.4  is  a  plot  of  the  radiation  pattern  of  a 
caged  dipole  over  earth.  There  is  a  very  close  resemblance  to 
the  radiation  pattern  of  a  thin  vertical  dipole  over  ground 
[Ref.  4:  p.  234].  Figure  7.5  is  a  comparison  of  radiation 
patterns  computed  by  ELNEC  and  Mathcad  software.  As  expected, 
the  patterns  are  similar. 
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FIGURE  7.4  E-plane  radiation  pattern  of  a  caged  dipole 
with  L  =  A/2  m,  N  =  8,  a  =  0.125  m,  h  =  0.25  m,  a  =  0.01 
S/m,  and  er  =  15 


OS  - 


0.6 


0.4 


0.2 


FIGURE  7.5  E-plane  radiation  pattern  comparison  between 
ELNEC  and  Mathcad  of  a  caged  dipole  with  L  =  X/2  m,  N  = 
8,  a  =  0.125  m,  h  =  0.25  m,  a   =   0.01  S/m,  and  er  =  15 
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C.  CAGED  DIPOLE  ORIENTED  HORIZONTALLY  OVER  EARTH 

The  following  known  or  estimated  inputs  are  additional 
data  required  to  analyze  a  caged  dipole  oriented  horizontally 
over  the  earth: 

h  =  antenna  height  above  ground 

er  =  relative  dielectric  constant  of  earth 

a   =   conductivity  of  earth 

ahw  =  incoming  wave  electric  field  unit  vector  for 
vertical  antenna 

0hp  =  antenna  polarization  direction  for  vertically 
orientated  antenna 

As  with  the  evaluations  of  a  caged  dipole  both  in  free 
space  and  oriented  vertically  over  the  earth,  the  analysis  of 
a  caged  dipole  positioned  horizontally  above  the  ground  does 
not  take  into  account  the  central  conductor  and  support 
members  since  they  produce  negligible  effects  on  the  antenna's 
electrical  properties. 

Through  the  use  of  image  theory,  the  far-field  electric 
field  intensity  of  a  caged  dipole  oriented  horizontally  above 
a  flat  earth  can  be  obtained.  The  8  component  of  the  electric 
field  of  a  horizontal  dipole  (parallel  to  the  y-axis)  in  free 
space  is  [Ref.  3:  pp.  143-144]: 
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-jkr. 


EQh(e,4>)  =jti0-2 


2rcr 


ff 


(^sin(6)sin(<J)))-cos/— 

v/l  -sin2  (6)  sin2(4>) 


(y/jB) 


(7.33) 


The  0  component  of  the  total  electric  field  over  ground 
is  [Ref.  4:  pp.  229-235]: 

EQht(B,$)  =EQh(e,$)  AFh(Q,<b)(ejkhcos{6)  -  Tve-jkhcos{6))        (V/m)    (7.34) 

The  array  factor  (AFh(0,0))  in  (7.34)  for  a  horizontal  array 
is  determined  by  [Ref.  3:  pp.  274  and  776]: 

AFhid.V    =  £  jne^asinesin*sin*«        (dimensionless)     (7.35) 

n  =  l 

In  (7.33)  through  (7.35),  the  electric  field  (E0h(8,0)),  the 
total  electric  field  (E0ht(0,0)),  and  array  factor  (AFh(0,0)) 
are  valid  only  above  the  ground  (0  <  0  <  tt/2)  [Ref.  3:  pp. 
135-142].  In  (7.34),  (h)  is  the  distance  from  the  earth  to 
the  center  of  the  antenna,  and  the  vertical  reflection 
coefficient  (Tv)  of  the  ground  is  given  by  (7.29). 

The  $  component  of  the  electric  field  of  a  horizontal 
dipole  (parallel  to  the  y-axis)  in  free  space  is  [Ref.  3:  pp. 
135-142] : 


-jkzt 


E4>h(Q,4>)  =jt|0-J 


2:rr 


ff 


cos(  —  sin(0)sin((J>))-cos(  — ) 
v^l-sin2(0)  sin2((J>) 


(V/m) 


(7.36) 


From  image  theory,  the  $  component  of  the  total  electric  field 
over  earth  is  [Ref.  4:  pp.  229-235]: 
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E$ht{e,$)  =E<S>h(d,4>)  AFh(d,<t>)(ejkhcos{e)  +  rhe-jkhcos{6))     (V/m)      (7.37) 

In  (7.37),  (AFh(0,0))  is  given  by  (7.35),  and  (Th)  is  the 
horizontal  reflection  coefficient  of  the  ground  and  is  [Ref. 
4 :  pp.  231-232] : 


t^    cos(6)  -  J?~-   sin2  (6)  .7,   .-,  ,fl, 

Th   =  v  {dimensionless)        (7.38) 

cos  (6)  +  ^e'r-  sin2  (6) 


When  a   -*   »,  rh  =  -1. 

The  radiation  intensity  (Uh(9,0))  of  a  caged  dipole  over 

a  earth  must  take  into  account  both  the  6  and  $  components  of 

the  electric  field.   Therefore,  the  radiation  intensity  is 

[Ref.  3:  p.  28] : 

r  2 
Uh(Q,$)  =^±-[\EQht(e,<\>)  |2  +  |£"<t)hc(e,(t>)  |2]  (w/ solid  angle)      (7.39) 

The  radiated  power  (Phrad)  is  present  only  above  the  earth. 
Thus, 

Phrad   =  [2"[K/ZUh  (6,4))  sin(  6)  d0d(J>    (W)  (7.40) 

Jo  Jo 

The  remaining  parameters  (i.e.,  Dh0,  EIRPh,  Rhin,  etc.) 
for  a  horizontal  caged  dipole  over  the  earth  are  calculated 
using  the  eguations  for  these  parameters  listed  previously. 

Table  7.3  is  a  comparison  of  the  gain  of  a  horizontal 
caged  dipole  over  a  flat  earth  obtained  with  ELNEC  and  the 
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Mathcad  application.  As  can  be  seen,  the  gains  are  similar. 
As  previously  mentioned,  the  difference  in  gains  can  be 
attributed  to  the  fact  that  ELNEC  takes  into  account  mutual 
coupling  between  the  elements  as  well  as  ground  effects.  The 
Mathcad  application  ignores  mutual  coupling  between  elements 
and  approximates  ground  effects. 

TABLE  7.3  Caged  Dipole  Antenna  Gain  Comparison 


ANTENNA 

ELNEC 

Mathcad 

DIMENSIONS 

Gain  (G) 

Gain  (G) 

L  =  A/2  m 

a  =  0. 025  m 

N  =  8 

4.99  dB 

4.73  dB 

h  =  1.0  m 

a  =  0.01  S/m 

er  =  15 

As  previously  mentioned  with  regard  to  the  vertically 
oriented  caged  dipole  over  the  earth,  the  caged  dipole  behaves 
in  a  manner  similar  to  that  of  a  thin  cylindrical  dipole. 
Therefore,  the  radiation  patterns  of  a  horizontally  oriented 
caged  dipole  are  similar  to  those  of  horizontal  thin 
cylindrical  dipole  which  in  turn  resembles  that  of  a  thin 
dipole.   Figures  7.6  and  7.7  illustrate  the  H-plane  and  E- 
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plane  patterns,  respectively,  of  a  caged  dipole  over  the  earth 
[Ref.  4:  pp.  231-237].  The  radiation  patterns  of  the 
horizontally  orientated  caged  dipole  over  a  flat  earth  are 
very  similar  to  those  of  a  horizontal  thin  dipole  over  earth. 
Figures  7.8  and  7.9  compare  the  radiation  patterns  between 
those  computed  with  ELNEC  and  those  computed  by  Mathcad 
software.   As  expected,  the  patterns  are  similar. 


FIGURE  7.6  H-plane  radiation  pattern  of  a  caged  dipole 
with  L  =  A/2  ™/  N  =  8,  a  =  0.125  m,  h  =  1.0  m,  a  =  0.01 
S/m,  and  er  =  15 
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FIGURE  7.7  E-plane  radiation  pattern  of  a  caged  dipole 
with  L  =  A/2  m,  N  =  8,  a  =  0.125  m,  h  =  1.0  m,  o  =  0.01 
S/m,  and  er  =  15 
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FIGURE  7.8  Comparison  between  ELNEC  and  Mathcad  of  H-plane 
radiation  patterns  of  a  caged  dipole  with  L  =  A/2  m,  N  = 
8,  a  =  0.125  m,  h  =  1.0  m,    o   =   0.01  S/m,  and  er  =  15 
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FIGURE  7.9  Comparison  between  ELNEC  and  Mathcad  of  E-plane 
radiation  patterns  of  a  caged  dipole  with  L  =  A/2  m,  N  = 
8,  a  =  0.125  m,  h  =  l.o  m,  a   =   0.01  S/m,  and  er  =  15 
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VIII.  PARABOLIC  REFLECTORS 

Many  communications  and  radar  operations  require  very 
large  values  of  antenna  gain  (G)  that  are  difficult  to  achieve 
with  a  single  device.  One  of  the  most  popular  ways  to 
construct  a  very  high  gain  antenna  is  to  reflect  the 
electromagnetic  energy  of  a  relatively  small  antenna  off  a 
metallic  parabolic  dish  in  the  far-field  of  the  feed  device 
[Ref.  13:  p.  12-2].  If  the  source  antenna  is  located  at  the 
focus  of  the  parabola  and  if  the  radius  (a)  of  the  parabolic 
dish  periphery  is  large  with  respect  to  wavelength  (A) ,  then 
the  reflector  may  be  approximated  as  a  uniform  phase  aperture 
antenna.  Gain  in  excess  of  30  dB  is  common  for  this  type  of 
structure  [Ref.  4:  pp.  423-424].  The  geometry  of  a  parabolic 
reflector  is  illustrated  in  Figure  8.1,  where  the  parameters 
include  the  focal  distance  (f),  the  distance  (p)  from  the 
origin  to  any  point  (P)  on  the  reflector  surface,  and  the 
distance  (r')  from  the  origin  to  the  projection  of  point  (P) 
onto  the  z  =  0  plane.  The  primed  angles  in  Figure  8.1 
correspond  to  the  parameters  associated  with  the  feed  horn. 
The  Mathcad  application  user  should  note  that  the  focal 
distance  should  exceed  the  minimum  distance  to  far-field  for 
the  feed  antenna. 
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FIGURE  8.1  Parabolic  Reflector  Geometry 


The  Mathcad  parabolic  reflector  applications  are  unique 
in  that  the  software  is  not  written  as  a  stand  alone  package. 
Rather,  the  reflector  software  is  run  as  an  addendum  to  the 
programs  that  already  exist  for  various  antennas  [Ref.  2:  pp. 
4-15,  67-114].  Four  types  of  feeds  are  analyzed:  helical, 
planar  equiangular  spiral,  pyramidal  horn,  and  conical  horn. 
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The  following  additional  inputs  are  required  to  analyze 
parabolic  reflectors: 

foe  =  reflector  focal  length 

aref  =  radius  of  the  mouth  of  the  reflector 

E0  =  electric  field  scale  factor 

eb  =  blockage  efficiency 

esp  =  spar  efficiency 

rff2  =  secondary  far  field  observation  distance 

*>y»z  =  polarization  loss  factor  coordinates 

aw  =  incoming  wave  unit  polarization  vector 

tl,  t2  =  radiated  power  increments  for  <p   and  6 

il  =  secondary  field  increments 

N,M  =  summation  increments 
The  first  two  inputs  are  physical  dimensions  obtained 
primarily  through  photographs,  inputs  with  the  are 
parameters  that  are  either  known  or  estimated,  and  the  last 
inputs  are  only  used  to  affect  the  processing  time  of  the 
application. 

In  order  to  calculate  the  far-field  patterns  of  a 
parabolic  reflector,  the  electric  field  at  the  aperture  of  the 
device  must  be  known.  Since  the  electric  field  patterns  of 
interest  are  in  the  far-field,  it  is  assumed  that  the  circular 
aperture  created  by  the  reflector  is  located  at  the  origin  and 
is  parallel  to  the  z  =  0  plane.  The  Cartesian  components  of 
the  aperture  electric  field  are  expressed  as  [Ref.  14:  pp. 
121-123;  Ref.  15:  pp.  29-31]: 
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E.x(r'/*'>  =  [-  E,f(e'(r')  ,0')cos(0') 

(8.1) 

+  E0f  (6'  (r;)  ,0/)sin(0/)  ]ejk("<r,>-2£ocl    (V/m) 


Eay(r',0')  =  [Eef(6/(r/),0/)sin(0/) 

(8.2) 

+  E^(e'(r')  ,0/)cos(0/)  ]ejk("(r,)-2foc)    (V/m) 


In  (8.1)   and  (8.2),   (Eax(r',8'))   and  (Eay(r'  ,6'))   are  the 
aperture  electric  fields  and  k  is  the  wavenumber  (2tt/X)  . 

The  relationships  between  (6'  )  ,  (p) ,  and  (r'  )  needed  to 
evaluate  (8.1)  and  (8.2)  are  expressed  as  [Ref.  4:  p.  426]: 


6'  =  2tan"1 


-^—\         {radians)  (8-3) 

2  foe. 


<r<)  =  ^f°c)^(r'Y         (m)  (8.4) 

A  foe 


The  electric  field  from  the  entire  parabolic  assembly  is 
commonly  referred  to  as  the  secondary  pattern.  Referring  to 
Figure  8.1,  the  electric  field  in  the  -z  half-plane  is 
considered  the  secondary  pattern.  The  secondary  field  of  a 
parabolic  reflector  is  calculated  using  the  electric  field 
integral  equation  solution  and  computing  the  component  vector 
potentials  (Px(6,0)),  (Py(6,0))  as  follows  [Ref.  4:  p.  426]: 

P    (0    <b)    =f'infaTef   E      d'    d)')  \eJkr'sin(e)(cos{^')cos(^)  *  sin  (♦')  sin  (♦))"] 

°       °  (8.5) 

*r'dr'db'  (V/m) 
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Pj<e'*>=J0  j0"r'Xy(r''*') 


e 


jkr'sin(0)(cos(<£'  )cos(<£)  ♦  sin(tf>'  )sin(<i)) 


(8.6) 


*r/dr/d0/  (V/m) 


Double  integrals  in  Mathcad  take  significantly  longer  to 
converge  to  an  answer  than  do  single  integrals  [Ref.  15:  p. 
200].  To  reduce  the  processing  time  of  the  applications, 
(8.5)  and  (8.6)  were  each  separated  into  two  single  integrals. 
Therefore,  for  the  applications,  (8.5)  and  (8.6)  are  each 
computed  as: 

P(0,4>)  =  [2nA(B,4>.4>')dV  (V/m)  (8.7) 

Jo 

where 

A  (6,  (J),  $')    =  fa"f   Ea(r/,^/)   eJ*r/sin(e)[cos«t>')cos(40  ♦  sin«|>')sin«l»]r/dr/    (8.8) 
J  0 

Any  contribution  from  feed  system  back  lobes  is  ignored, 
and  the  component  vector  potentials  calculated  in  (8.7)  and 
(8.8)  are  applied  to  the  subsequent  formulas  to  compute  the 
secondary  field  [Ref.  4:  p.  383]: 

"I  led  ~3krtf2 

£e(8,4>)  =^p [Px(8,40cos(40  +  Py  (8, 40  sin  (40]  (v/m)     (8.9) 

^(6,4>)  =J4f [-Px(6,40sin(40  +  Py(8,40  cos  (40]cos  (6)       (8.10) 

2  n  rff2 
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The  radiation  intensity  (U(6,0)),  radiated  power  (Prad)  , 
and  directivity  (D0)  of  the  secondary  field  pattern  are 
expressed  using  the  intrinsic  impedance  of  free  space  (r)0  = 
1207T)  as  follows  [Ref.  3:  pp.  27-3.0]: 

£7(6,4))  =  {l'f2)    [|ge(e'4>>  l2+  l^(9'4>)  I2]    W solid  angle)       (8.11) 

Prad    =    f2nf*/2U(e,4>)sin(d)ded<k  (W)  (8.12) 


Dq   =  i^ax    (dimensionless)  (8.13) 

Prad 


In  (8.13),  (Umax)  is  the  maximum  value  of  radiation  intensity, 
obtained  by  creating  (Umax)  as  a  matrix  and  finding  the  maximum 
value  of  the  matrix. 

The  computer  processing  time  required  to  evaluate  (8.12) 
is  prohibitively  long.  Thus,  the  parabolic  reflector  Mathcad 
applications  use  the  following  numerical  approximation  to 
evaluate  radiated  power  [Ref.  3:  pp.  37-42]: 

N     (    M  \ 


N  A  2M/JT! 


p*>*-[^ \\-h\L  E^wsinie,) 


\A  =  \ 


(W)  (8-14) 


In  (8.14),  (N  and  M)  are  user  selected  integers.  This 
approximation  for  Prad  reduces  the  processing  time  of  the 
applications  substantially.  As  N  and  M  increase,  the  results 
obtained  with  (8.14)  converge  to  that  obtained  with  (8.12). 
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When  N  =  M  =  3,  the  processing  time  for  only  Prad  of  the 
helical  feed  reflector  option  is  approximately  four  hours  on 
a  33  MHz,  386  machine,  and  when  N  and  M  is  increased  to  5,  the 
time  to  calculate  Prad  for  the  same  application  is  ten  hours. 
For  the  applications,  we  chose  N  =  M  =  10  to  evaluate  Prad  with 
reasonable  accuracy.  With  this  choice,  it  takes  approximately 
one  day  for  the  application  to  run,  but  this  is  much  less  than 
the  thirty-six  hours  it  takes  if  (8.12)  is  used. 

Parabolic  reflectors  are  members  of  the  aperture  antenna 
family.  The  gain  of  an  aperture  antenna  is  the  product  of  its 
aperture  efficiency  (eap)  and  directivity.  Therefore,  the  gain 
of  a  parabolic  reflector  is  written  as  [Ref.  6:  p.  575]: 

G   =  eapD0        (dimensionless)  (8.15) 

Aperture  efficiency  is  the  product  of  several  separate 
terms,  including  spillover  efficiency  (es)  ,  taper  efficiency 
(et) ,  phase  (or  random  surface  error)  efficiency  (ep)  , 
polarization  efficiency  (ex)  ,  blockage  efficiency  (eb)  ,  spar 
efficiency  (esp)  ,  and  ohmic  efficiency  (eohmic)  . 

Spillover  efficiency  (es)  is  the  fraction  of  power 
radiated  by  the  feed  system  that  is  intercepted  and  reflected 
by  the  parabolic  dish.  Spillover  efficiency  is  the  most 
important  efficiency  term  for  any  feed  system  and  is  defined 
as  [Ref.  6:  pp.  583] : 
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f     f  °Uf  (0')  sin  (00  dQ'dtf 

Jo    Jo 


es" 


f    rtueOsinoOde^ 

Jo  Jo 


(dimensionless)         (8.16) 


In  (8.16),  (Uf(6'  ))  is  the  radiation  intensity  function  of  the 
feed  antenna  and  (0O)  is  the  value  of  (8.3)  when  (r'  )  equals 
the  radius  of  the  mouth  of  the  reflector  (aref)  . 

Taper  efficiency  (et)  accounts  for  the  lack  of  amplitude 
uniformity  of  the  feed  pattern  on  the  surface  of  the 
reflector.  Taper  efficiency,  which  is  a  dimensionless 
quantity,  is  expressed  as  [Ref.  3:  pp.  626-627]: 


et  =  2cot: 


l  2 


Jo  Jo  \      piad         v  2 ; 


(8.17) 


/   I — sin  (8')  de'dcb' 


Phase  efficiency  (ep)  is  a  function  of  the  power  loss 
that  occurs  if  the  field  at  the  mouth  of  the  dish  is  not  in 
phase  at  every  point  in  the  aperture.  Phase  efficiency  is 
approximately  given  by  [Ref.  4:  pp.  434-435]: 


4n610-5aref 


ep  -  e 


(dimensionless) 


(8.18) 


The  power  that  is  lost  due  to  cross-polarized  electric 
fields  in  the  aperture  of  the  reflector  determines 
polarization  efficiency   (ex) .   Polarization  efficiency  is 
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difficult  to  precisely  calculate,   but  can  be  reasonably 
approximated  as  0.98  [Ref,  4:  p.  4  35]. 

Blockage  (eb)  and  spar  (esp)  efficiencies  are  determined 
by  the  presence  of  feed  structures,  support  struts,  and 
possible  signal  processing  equipment  in  front  of  the 
reflector's  aperture.  Precise  computation  of  blockage  and 
spar  efficiencies  requires  an  extensive  method  of  moments 
simulation  and  is  not  attempted  by  the  Mathcad  parabolic 
reflector  applications.  However,  if  the  radius  (af)  of  the 
blocking  structure  at  the  focal  point  and  the  number  of  spars 
(N)  that  are  (A/2)  thick  is  known,  blockage  efficiency  and 
spar  efficiency  can  be  estimated  using  Tables  10.1  and  10.2, 
respectively  [Ref.  4:  p. 436]. 

TABLE  8.1  Blockage  Efficiency 


af/aref 

0.05 

0.  10 

0.20 

Eb 

0.99 

0.96 

0.84 

TABLE  8.2  Spar  Efficiency 


aref   =    5/ A 

aref   =    50/A 

aref    =    100/A 

N    = 

3 

0.95 

0.99 

1.00 

N    = 

4 

0.94 

0.99 

1.00 
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Ohmic  efficiency  (eohmic)  for  a  reflector  is  usually  very 
high  but  difficult  to  accurately  calculate.  The  Mathcad 
parabolic  reflector  applications  assume  ohmic  efficiency  is 
0.98. 

With  all  reflector  efficiency  terms  known  or  estimated, 
the  aperture  efficiency  (eap)  ,  effective  aperture  (Ae)  ,  and 
isotropic  radiated  power  (EIRP)  are  given  by  [Ref.  3:  pp.  623- 
630]  : 

taP=  tstttptxtbtsPt  ohmic        (dimensionless)  (8.19) 

Ae  =  tapTt(aref)2        (m2)  (8.20) 

EIRP  =   PiadD0         (W)  (8.21) 

The  magnitude  of  the  feed  current  (I0)  for  the  source 
antenna  cannot  be  estimated  based  on  measured  geometry. 
However,  if  the  feed  current  is  normalized  to  1  amp,  radiation 
resistance  (Rr)  for  the  entire  reflector  structure  is: 

l-Tol 

In  (8.22),  (Prad)  is  the  radiated  power  from  the  aperture  of 
the  reflector,  not  of  the  feed  antenna. 

The  Mathcad  parabolic  reflector  applications  use  the 
radiation  resistance  computed  for  the  parabolic  reflector 
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antenna  system  to  determine  effective  height  (hem)  by  [Ref.  6: 
p.  42]: 


he  =   2 


^r£e     {m)  (8.23) 


The  bandwidth  of  a  parabolic  reflector  is  assumed  to  be 
the  same  as  its  feed  antenna. 

The  polarization  loss  factor  (PLF)   of  the  parabolic 
reflector  for  a  user  selected  point  in  the  secondary  far-field 
and  a  specific  wave  unit  polarization  vector  (aw)  is  found 
using  the  electric  field  components  of  (8.9)  and  (8.10)  as 
follows  [Ref.  3:  51-53]: 

PLF  =  |a,/a*|2    (dimensionless)  (8.24) 

where  the  unit  polarization  vector  of  the  antenna  is: 

aE  +  aE  +  aE 
oa(x,y,  z)   =  x    x        y    y        - — -    (dimensionless)       (8.25) 

y/\E(x,y,z)  | 

and  the  Cartesian  components  of  the  electric  field  are  [Ref. 
3:  p.  776]: 

Ex=   £"Qcos(6)cos((J))  -  £^sin(4>)    (V/m)  (8.26) 

Ey=   £,ecos(0)sin(4>)  +  fi^eos  (<|>)    (V/m)  (8.27) 

Ez=   -£,Qsin(6)    (V/m)  (8.28) 
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All  of  the  secondary  field  parameters  calculated  by  the 
Mathcad  parabolic  reflector  applications  assume  that  the 
observation  point  is  in  the  far-field.  Thus,  all  of  the 
subsequent  formulas  must  hold  for  the  application  calculations 
to  be  valid  [Ref.  4:  pp.  24-25]: 

i   *  1.6X     (m)  (8.29) 

r  z  10aref  (m)  (8.30) 

8(3.        \^ 

r  z      ^  ref'  (dimensionless)  (8.31) 

Table  8.3  is  a  comparison  of  reflector  directivity 
computed  by  the  Mathcad  applications  to  predictions  made  for 
an  identical  antenna  with  unity  aperture  efficiency  [Ref.  3: 
p.  634].  As  can  be  seen,  the  comparisons  vary  by 
approximately  8  dB.  This  can  be  attributed  to  the  fact  that 
Prad  was  approximated  using  (8.14)  instead  of  (8.12),  and  also 
to  the  fact  that  the  numerical  tolerance  variable  (TOL  in 
Mathcad's  math  menu)  for  the  application  was  increased  to  0.5 
verses  0.001  in  order  to  reduce  processing  time.  This 
increase  in  TOL  reduces  both  the  accuracy  of  (8.14)  and  the 
processing  time  of  the  application. 
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TABLE  8.3  Parabolic  Reflector  Data  Comparison 


ANTENNA 
PARAMETER 

PREDICTED 
DATA 

CALCULATED 
DATA 

DIRECTIVITY 

42  dB 

50.2  dB 

Figure  8.2  is  a  comparison  of  measured  data  fed  by  an 
optimum  horn  with  that  calculated  by  the  Mathcad  application 
fed  by  an  equiangular  spiral  for  a  parabolic  reflector  (aref  = 
0.320  m,  f  =  0.147  m)  operating  at  35  GHz  [Ref.  14:  p.  126]. 
As  can  be  seen,  the  comparison  is  good.  The  Mathcad  pattern 
of  Figure  8.2  is  not  symmetric  which  can  be  attributed  to  the 
number  of  increments  used  to  increase  the  speed  of  processing. 

The  reflectors  fed  by  pyramidal  and  conical  horn  antennas 
are  not  included  in  this  report.  The  applications  have  been 
developed  and  are  included  on  the  diskette  containing  the 
software.  The  complexity  of  the  software  results  in 
prohibitively  long  processing  time  for  the  two  applications 
(i.e.,  the  pyramidal  horn  fed  reflector  application  was  still 
calculating  Prad2  after  seven  days  when  run  on  a  50  MHz,  486 
personal  computer  using  Mathcad  4.0). 
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-  Mathcad  Pattern,  ...  Meosured  Pattern 


FIGURE  8.2   Comparison  of  Parabolic  Reflector  E-Plane 
Electric  Field  Patterns 
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IX,  REMARKS  AND  CONCLUSIONS 

This  thesis  and  the  accompanying  software  applications 
are  intended  to  provide  the  Naval  Maritime  Intelligence  Center 
(NAVMARINTCEN)  with  software  for  the  analysis  of  various 
antenna  types  where  information  regarding  the  antenna  is 
restricted  primarily  to  the  antenna  type  and  the  physical 
dimensions  of  the  antenna.  The  goal  of  this  research  was  to 
provide  user  friendly  software  that  an  engineer  familiar  with 
basic  antenna  types  can  easily  use  and  interpret.  Although 
the  software  was  developed  for  that  purpose,  the  complexity  of 
formulas  that  apply  to  many  of  the  antenna  types  of  necessity 
reduced  the  simplicity  of  the  applications. 

The  software  applications  developed  are  written  in 
Mathcad  3.1  and  are  compatible  with  any  IBM  personal  computer 
that  supports  Mathcad  3 . 1  or  Mathcad  4.0  for  Windows.  Many  of 
the  applications  require  extensive  processing  time.  In  order 
to  reduce  processing  time,  which  results  from  the  complexity 
of  various  integrals,  approximate  numerical  techniques  were 
implemented  to  evaluate  the  integrals  and  the  numerical 
tolerance  of  the  more  complicated  applications  was  reduced. 
Even  with  the  processing  time  reduced,  some  programs  require 
up  to  five  days  to  complete  when  run  with  Mathcad  4.0  on  a  50 
MHz,  486  personal  computer.  Applications  that  did  not  exceed 
a  day  to  complete  were  run  on  a  33  MHz,  386  personal  computer 
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with  Mathcad  4.0.  All  applications  were  run  using  Mathcad  4.0 
since  Mathcad  4.0  is  approximately  twice  as  fast  as  Mathcad 
3.1.  As  previously  mentioned,  the  applications  were  written 
in  Mathcad  3.1  due  to  sponsor  requirements.  Mathcad  3.1  is 
compatible  with  Mathcad  4.0,  but  Mathcad  4.0  is  not  compatible 
with  Mathcad  3.1. 

For  this  thesis,  a  number  of  the  antennas  could  not  be 
adequately  researched  through  the  use  of  texts.  In  some 
cases,  such  as  the  caged  dipole,  very  little  information  on 
the  antenna  could  be  found  in  either  texts  or  in  professional 
journals.  For  the  caged  dipole,  the  antenna  is  modeled  as  a 
circular  array  of  center-fed,  equally  excited,  uniformly 
spaced,  thin-wire  dipoles  to  take  into  account  the  number  of 
conductors  that  are  present  in  the  caged  dipole.  This  model 
provides  good  results  when  compared  to  predicted  (the  caged 
dipole's  performance  is  similar  to  that  of  a  thick  cylindrical 
dipole)  and  ELNEC  results. 

It  was  found  that  Mathcad  is  not  the  ideal  tool  to 
analyze  antennas  due  to  its  long  processing  time.  An 
additional  drawback  of  Mathcad  is  that  when  a  function  is  in 
the  beginning  of  an  application  and  the  same  function  is 
placed  later  in  the  application,  Mathcad  does  not  store  the 
values  of  the  function,  but  re-evaluates  the  function.  This 
re-evaluation  increases  processing  time  significantly, 
especially  when  the  function  itself  takes  a  very  long  time 
(hours  vice  minutes)  to  reach  a  numerical  solution. 
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Nevertheless,  the  Mathcad  applications  developed  in  this 
research  allow  the  user  to  analyze  several  antenna  types 
(i.e.,  arrays  with  isotropic  point  sources,  parabolic 
reflectors  with  the  various  feed  antennas)  that  could  not 
easily  be  evaluated  using  some  antenna  numerical  analysis 
programs  such  as  NEC,  ELNEC,  or  MININEC.  For  the  other 
antenna  types  researched,  the  antenna  numerical  analysis 
programs  mentioned  would  be  easier  to  setup,  run,  and  provide 
parameters  closer  to  those  of  a  real  antenna.  The 
disadvantage  of  this  approach  is  that  it  requires  an  ability 
to  evaluate  antennas  using  a  program  such  as  ELNEC,  while  the 
Mathcad  applications  do  not  require  any  advanced  programming 
knowledge  on  the  part  of  the  user. 
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"X> 


LINEAR  ARRAY 
MATHCAD  SOFTWARE  ARRAY_LN . MCD 

This  linear  array  application  is  based  on  N  isotropic  point  sources 
uniformly  spaced  and  aligned  on  the  z-axis.   Mutual  interference  between 
adjacent  sources  are  ignored  when  calculating  antenna  parameters.   Since 
isotropic  point  sources  are  used  as  the  radiating  elements,  polarization 
of  the  antenna  can  not  be  determined.   The  gain  can  only  be  idealized 
when  it  is  assumed  that  the  antenna  input  resistance  is  equivalent  to 
the  radiation  resistance.   Bandwith  calculations  are  not  carried  out  as 
a  result  of  the  narrowband  characteristics  of  arrays,  and  in  large  part 
to  the  dependance  inter-element  spacing  (d)  has  with  the  wavelength  (X) . 
The  following  limits  should  be  adhered  to  when  analyzing  Broadside, 
End-Fire,  and  Phased  (scanning)  antennas: 

0  <  80  <  7i   (Scanning  Antennas) 

=  Tt/2   (Broadside  Antennas) 

0O  =  0  or  7t   (Ordinary  End-Fire  Antennas) 

(Note:   Mathcad  equations  can  not  use  symbolic  subscripts.   Therefore, 
parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 

The  linear  array  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  *  indicate  parameters  that  can  be  computed  but 
their  results  are  idealized  as  the  result  of  using  isotropic  point 
sources  as  the  radiating  elements): 

X      =      Wavelength 

rmin   =   Minimum  Distance  to  the  Far-Field 

P   =   Phase  Difference  Between  Adjacent  Elements 

y   =   Array  Factor  Phase  Shift 

AFn   =   Normalized  Array  Factor 

U(0)   =  Radiation  Intensity 
Pracl   =   Radiated  Power 

DQ   =   Directivity 

EIRP   =   Effective  Isotropic  Radiated  Power 
Rr   =   Radiation  Resistance 

Rj_n   =   Antenna  Input  Resistance  * 

r   =   Voltage  Reflection  Coefficient  * 

et   =   Antenna  Efficiency  * 

G   =   Gain  * 

Aem   =   Maximum  Effective  Aperture  * 

hem   =   Maximum  Effective  Height  * 
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The  following  known  or  estimated  data  must  be  entered: 

N  =  Number  of  Isotropic  Radiating  Elements 

f  =  Frequency  of  Interest 

d  =  Distance  Between  Adjacent  Elements 

0O  =  Direction  of  Main  Lobe 

IQ   =   Antenna  Feed  Current 

Z0   =   Characteristic  Feed  Impedence 


Enter  input  data  here 
N  :=  10  (elements)  d:=0.25  (meters) 


f :  =  3  1 0  (Hz)  9o:=0  (radians) 


Io:=1.0  (A)  Zo:=75.0         (Q) 


Define  constants  and  calculate  wavelength: 

c:=  2.9979  108      (meters/sec)  t]0  :=  120-ji        (Q) 


X'=-  (meters/cycle)  k  ~ —  (m  ■*■) 


X  =0.999  (meters/cycle)  ccd  :=  1  (dimensionless ) 
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Calculate  linear  array  parameters 


Define  angular  offset  9: 


e:=(o-io-4)i^  .  (2  Ti-  10- 4) 

W 


( radians; 


Minimum  Distance    to    Far-Field    r 


mm  j. 


r  :=  1.6-X  (meters) 


r,:=5Nd 


(meters 


2(Nd)' 


(meters ) 


rmin:=max(r)  (meters 


rmin  =  12.509 


(meters! 


Array    Factor  AFn : 


P  :=-kdcos(8o) 


( radians) 


(|»(8):=kdcos(e)+p 


(radians) 


Afn(0)  :=  — 
N 


sin|—  y(0) 


sin|—  >|;(9) 


(dimensionless) 
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Radiation  Intensity  U: 


U(0):=(AFn(6)) 


(W/solid  angle) 


Radiated  Power  P 


rad-L 


Prad:=2  7i 


U(0)  sin(6)de 


(W) 


Prad  =  1 .256 


(W) 


Directivity   D^: 


Umax  •=  1 


(W/solid    angle) 


Do:= 


4  7i  Umax 
Prad 


(dimensionless) 


Do  =  10.007 


(dimensionless) 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :=  Prad  Do 


(W) 


EIRP  =  12.566 


(W) 
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Radiation    Resistance    Rr: 


(Holy 


Rr  =  2.512  (Q) 


Input  Resistance  Rjn: 

Rin:=Rr  (Q) 


Rin  =2.512  (Q) 


Voltage  Rejection  Coeffecient  T: 

_ ._  Rm  -  Zo  .         " 

I  •= (dimensionless) 

Rin+  Zo 


f  =-0935  (dimensionless) 


Reflection  Efficiency  tt : 

et  =  ecd  [  1  —  ( |r| )  J  (dimensionless) 


et  =0.125  (dimensionless) 
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Gain    G: 


G=et  Do 


(dimensionless ) 


GdB-=  10  1og(G)         (dB) 


G  =  1.255 


(dimensionless ) 


GdB=0  986 


(dB) 


Maximum   Effective   Aperture   AP1T1: 
,2 


Aem  := 


G£ 

4-71 


[m< 


Aem  =0.1 


m< 


Maximum   Effective    Height    h^m: 


hem  :=2 


f  Rr  Aem 


T]0 


(m) 


hem  =  0.052 


(m) 
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NORMALIZED  ARRAY  FACTOR  RECTANGULAR  PLOT 

Number  of  radiating  elements: 
N  =  10 

Direction  of  Main  Lobe  (radians): 
0o=  0 


0.75 


JAFn(  9)| 


0.5 


0.25 


0  0.628  1.257  1.885  2.513  3.142  3.77  4.398  5.027  5.655       6.283 

e 
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NORMALIZED  ARRAY  FACTOR  POLAR  PLOT 

Number  of  Elements:  N  =  10 

Direction  of  Main  Beam  (radians)  :     80  =  0 


X(8):=  |AFn(G)|-cos[e  +  - 


Y(G):=  |AFn(0)|  sin  0+  - 


=   0° 


1 

s                                                             J,                    /                 \                                       ,                                       '                   \                   >^                                                            ' 

tt*  \    /        ,             ;             ;         \    /    -. 
K               K   , -"  ',          !           •    "^v.   x                '\ 

*     ,                      *         \     *           1          '          '           'it                     ->    " 

/              x                    >                    x         \           x."           ,              .               >-.(/'                    .                ^XN 

>                   "  \       *        \       '        '         '       A  ' 

y^m 

0 

I                                     1                 ^T^^ffl^r-:     '                '                                                                          ' 

Y(_6) 

•  ,---:      ,'\,';:\\\\y-*       7-..  ;          ,; 

x     /    .'    :    \    \     ;< 

-x  *        ■       ■       ,        »  ,' 

'^                                            '          "  ^          .                           '                           '           ^   '          *                                        N       ' 
V                                       x^_'                          1                          .     _'                 \                                      > 

/      %                                       >                                  *•  •<-  ->                1      .       .   j.  -r    ,                                   K                                      f       \ 

x           s                      ,                        1                 -  7  -                                         >                                 N 

*         '                              '                           '                            .                               »       ^  ' 

'                                                t             --                    >                              1                                                  <-              ^                                                s 

X                                                                    (•                                  ^             C                                           1                                           '              ^                                   1                                                                    N 

-1 

1                                                                                                                                                     ->    ■                                         *                    \                                                                                                                                1 

-1 


xref 


0 
,,X(e) 

9   =    180° 
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PLANAR  ARRAY 
MATHCAD  SOFTWARE  ARRAY_PL . MCD 

This  planar  array  application  is  based  on  M  and  N  isotropic  point 
sources  uniformly  spaced  and  aligned  on  the  x-axis  and  y-axis 
respectively.   Mutual  interference  between  adjacent  elements  are 
neglected  when  calculating  antenna  parameters.   Since  isotropic  point 
sources  are  used  as  the  radiating  elements,  polarization  of  the  antenna 
can  not  be  determined.   The  gain  can  only  be  idealized  when  it  is 
assumed  that  the  antenna  input  resistance  is  equivalent  to  the  radiation 
resistance  and  the  total  antenna  efficiency  is  calculated  for  a  ideal 
lossless  antenna. 

(Note:   Mathcad  equations  can  not  use  symbolic  subscripts.   Therefore, 
parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 

The  planar  array  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  *  indicate  parameters  that  can  be  computed  but 
their  results  are  idealized  as  the  result  of  using  isotropic  point 
sources  as  the  radiating  elements): 


Lx   =   Antenna  Dimension  in  x-direction 


X      =   Wavelength 
Lx   =   Antenna  D 

Ly   =   Antenna  Dimension  in  y-direction 

rmin   =   Minimum  Distance  to  the  Far-Field 

Px   =   Progressive  Phase  Shift  Between  Adjacent  Elements  in  x-direction 

(3y   =   Progressive  Phase  Shift  Between  Adjacent  Elements  in  y-direction 
i|»x   =   Array  Factor  Phase  Shift  in  x-direction 

H»y   =   Array  Factor  Phase  Shift  in  x-direction 

AFn(0,(|i)  =   Normalized  Array  Factor 

U  ( 9,<}> )   =  Radiation  Intensity 
Pracl   =   Radiated  Power 

D0   =   Directivity 

EIRP   =   Effective  Isotropic  Radiated  Power 
Rr   =   Radiation  Resistance 

R^n   =  Antenna  Input  Resistance  * 

T   =   Voltage  Reflection  Coefficient  * 
et   =   Antenna  Efficiency  * 
G  =   Gain  * 


Aem   =   Maximum  Effective  Aperture 
hem   =   Maximum  Effective  Height  * 
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The  following  known  or  estimated  data  must  be  entered: 

M   =   Number  of  Isotropic  Radiating  Elements  in  x-direction 
N   =   Number  of  Isotropic  Radiating  Elements  in  y-direction 
f   =   Frequency  of  Interest 
dx   =   Distance  Between  Adjacent  Elements  in  x-direction 

dy   =   Distance  Between  Adjacent  Elements  in  y-direction 

0O   =   Direction  of  Main  Beam  at  0  =  0o 

$0   =   Direction  of  Main  Beam  at  §   =  $o 

I0   =   Antenna  Feed  Current 

ZQ   =   Characteristic  Feed  Impedence 

j   =   Number  of  Increments  in  Degrees  for  Far-Field  Radiation  Pattern 


Enter  input  data  here 
M  :=  5  (elements)  dx:=0.5  (meters) 


N  =5  (elements)  dy  =0.5  (meters) 


f  :=  3  108  (Hz)  0o:=-  (radians) 


Io  =  1.0  (A)  fo   =-  (radians) 

4 


Zo:=75.0  (fi)  j  :=  180  (degrees) 
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Define  constants  and  calculate  wavelength: 


c:=  2.9979  108      (meters/sec) 


T|o  :=  1 20  n 


(Q) 


l:=* 


(meters/cycle) 


k:= 


2-71 


(m"1) 


X  =0.999 


[meters/cycle; 


ecd  :=  1 


(dimensionless) 


Calculate  linear  array  parameters 


Define  angular  offset  9_andjfe: 


e  :(o-io-4),[-io-4+i?]  ..(2,.  icr4) 


(radians) 


♦  :=(o-icr4),(.i(r4+ii!J..(2n-lor4) 


( radians] 


Minimum   Distance    to    Far-Field    r 


miru 


r  :=  1.6-X  (meters; 


r  :=5-/(M-dx)2t  (N-dy)2 


(meters) 


,_2  /J(Mdx)2^-(Ndy)2- 


(meters) 


rmin  :=  max( r)  (meters) 


rmin  =25.018  (meters) 
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Array  Factor  AFn : 


Px  :=-k-dx-sin(9o)-cos(<jio) 


( radians) 


Py  '=-kdysin(0o)sin(<J>o) 


(radians ) 


x|0<(  G ,  4> )  :  =  kdxsin(0)  cos(<j>)  +  px 


( radians) 


ipy(G,4)  -  =  kdysin(0)  sin(iji)+  Py 


radians) 


AFn(9,<J0  :« 


sin 

(t 

■vx(e,*)J 

M- 

sin 

VX(6,*)\ 

n|—  <fy(M) 


sJbesm 


(dimensionless) 


Radiation  Intensity  U: 


U(G,<jO  .=  AFn(0»-AFn(e,4) 


(W/solid  angle) 


Radiated  Power  Prar): 
f  2-Ti 


Prad  := 


0 


U(0,<|>)  sin(0)d0d<)> 


0 


(W) 


Prad  =0.206 


(W) 
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Directivity  Dn: 

Umax  :=  1  (W/solid  angle) 


^  .  4  n  Umax 

Do  •= (dimensionless) 

Prad 


Do  =61.09  (dimensionless) 


Effective  Isotropic  Radiated  Power  EIRP: 

EIRP  =  Prad  Do  (W) 


EIRP  =  12.566  (W) 


Radiation  Resistance  Rr: 

Rr:  =  2-^  (Q, 


C  |lo|  >- 


Rr  =0.411  (Q) 


Input  Resistance  R-jn: 

Rm:=Rr  (Q) 


Rin  =0.411  (Q) 
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Voltage  Rejection  Coeffecient  f: 
Rin  -  Zo 


r-= 


Rin  +  Zo 


(dimensionless ) 


r  =-0.989 


(dimensionless; 


Reflection  Efficiency  zr 
et:=«xi[l-(|r|)2] 


(dimensionless) 


et  =0.022 


(dimensionless 


Gain  G: 


G:  =  et  Do 


(dimensionless ) 


GdB:=  101og(G)    (dB) 


G  =  1.326 


[dimensionless) 


GdB  =  1.225 


(dB! 


Maximum  Effective  Aperture  A^m: 
,2 


Aem  := 


or 

4  71 


(m^ 


Aem  =0  105 


(m2) 


Maximum  Effective    Height    hpm: 

I       .    .    /Rr  Aem 

hem  •=  2     (m) 


TJO 


hem  =  0.021 


(m) 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:  (x-z)  plane 

Number  of  Elements  along  x-axis:  M=5 
Number  of  Elements  along  y-axis:  N=5 
Direction  of  Main  Beam:  0o  =  0.524 

$o  =  0.785 


( radians) 
(radians) 


X1(0)=  |AFn(8,0)|-cos(e-- 


Yl(e)=|AFn(e,0)|sin|e-- 


=  0° 


-  -   0 

Y1(G) 


-1 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:   (y-z)  plane 


Number  of  Elements  along  x-axis 
Number  of  Elements  along  y-axis 
Direction  of  Main  Beam: 


X2(8)  = 


AFnle,- 


COS|W 


M=5 
N=5 
0o  =  0.524 
<t»o=0  785 

Y2(9)  = 


AFnle,- 


( radians) 
( radians) 

W-  - 


e  =  o( 


3™^ 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:  &  =  7t/4 
Number  of  Elements  along  x-axis: 
Number  of  Elements  along  y-axis: 
Direction  of  Main  Beam: 


M  -5 
N=5 
00  =  0.524 
<(>o  =  0.785 


( radians) 
( radians) 


X3(0)  := 


AFn|9,- 


IU 


Y3(9):  = 


AFn|6,- 


sin  « 


=    0° 


1 

\                                         %                                       —  *~             i              --  —       i                              i 

/ 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:  & 

Number  of  elements  along  x-axis: 
Number  of  Elements  along  y-axis: 
Direction  of  Main  Beam: 


-   *c 


M=5 

N=5 
0o  =  O524 
<fro  =  0.785 


(radians) 
( radians ) 


X4(6)  :=  |AFn(6,*o)|  cos|9-  - 


Y4(e):=|AFn(9,*o)|  sin  6-- 


=    0C 
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CIRCULAR  ARRAY 
MATHCAD  SOFTWARE  ARRAY_CI.MCD 

This  circular  array  application  is  based  on  N  isotropic  point  sources 
uniformly  spaced  on  the  x-y  plane  along  a  circular  ring  of  constant 
radius.   Mutual  interference  between  adjacent  elements  are  neglected 
when  calculating  antenna  parameters.   Since  isotropic  point  sources  are 
used  as  the  radiating  elements,  polarization  of  the  antenna  can  not  be 
determined.   The  gain  can  only  be  idealized  when  it  is  assumed  that  the 
antenna  input  resistance  is  equivalent  to  the  radiation  resistance  and 
the  total  antenna  efficiency  is  calculated  for  a  ideal  lossless  antenna. 


(Note:   Mathcad  equations  can  not  use  symbolic  subscripts.   Therefore, 
parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 


The  circular  array  antenna  Mathcad  application  will  compute  thefollowing 
parameters  (Items  with  *  indicate  parameters  that  can  becomputed  but 
their  results  are  idealized  as  the  result  of  usingisotropic  point 
sources  as  the  radiating  elements): 

X      =      Wavelength 

rmin   =   Minimum  Distance  to  the  Far-Field 

an   =   Phase  Excitation  of  the  nth  Element 
AF(0,4O  =   Normalized  Array  Factor 
U  (0,4>)   =  Radiation  Intensity 
Pracj   =   Radiated  Power 

DQ   =   Directivity 

EIRP   =   Effective  Isotropic  Radiated  Power 
Rr   =   Radiation  Resistance 

Rin   =   Antenna  Input  Resistance  * 

r   =   Voltage  Reflection  Coefficient  * 

et   =   Antenna  Efficiency  * 

G   =   Gain  * 

Aem   =   Maximum  Effective  Aperture  * 


hem   =   Maximum  Effective  Height  * 
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The  following  known  or  estimated  data  must  be  entered: 

N  =  Number  of  Isotropic  Radiating  Elements 

f  =  Frequency  of  Interest 

a  =  Radius  of  circle 

0O  =  Direction  of  Main  Beam  at  8  =  Go 

4>0   =   Direction  of  Main  Beam  at  <j>  =  tyo 

IQ   =   Antenna  Feed  Current 

Z0   =   Characteristic  Feed  Impedence 

j   =   Number  of  Increments  in  Degrees  for  Far-Field  Radiation  Pattern 


Enter  input  data  here 
N  :=  10  (elements)  a  := —  (meters) 

2  71 


f  :=3  108  (Hz)  0o:=O  (radians) 


k>:=1.0  (A)  <|kd  :=  0  (radians) 


Zo:=75.0  (Q)  j:=360  (degrees) 


Define    constants    and   calculate    wavelength: 

c  :=  2.9979  108  (meters/sec)  no  :=  1 20  ti  (Q) 


X:=-  (meters/cycle)  k  := —  (m      , 
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Define  constants  and  calculate  wavelength: 


X  =  0.999 


(meters /cycle) 


ecd  :=  1 


(dimensionless) 


Calculate  linear  array  parameters 


Define  angular  offset  8_andjji: 


e:=(o-lff4)i.lO-4+M  (2  k-  10-4) 


radians ) 


M(o-io-4)i-io-4+i-j  (2,-KT4) 


radians) 


Minimum  Distance  to  Far-Field  r 


minj 


rQ  =1.6-1 


(meters) 


r,  =5(2  a) 


(meters) 


2-(2-a)' 


(meters ) 


rmin  :=  max(r)       (meters; 


rmin=  20.278 


(meters) 


Array  Factor  AF: 

n=l..N  (elements) 


0  =2  7i- 


N 


( radians) 


a  :=  -k  a  sin(0o)  cosf tyo-  <b  \ 


radians) 
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Array  Factor  AF: 


AF(8,*)  := 


N 


j  (kasin(e)cos(iti- 0^+a^ 


(dimensionless) 


Radiation  Intensity  U: 


U(e»:=AF(8,t)-AF(e,t) 


(W/solid   angle) 


Radiated    Power    Prarj: 

2  71 


Prad  := 


•0      • 


U(0,<j>)  sin(0)d9d(t. 


0 


(W) 


Prad  =0.536 


(W) 


Directivity  Dn: 


Umax  :  =  1 


(W/solid  angle) 


Do.= 


4  n-  Umax 
Prad 


(dimensionless) 


Do  =23.427 


(dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :=  Prad  Do 


(W) 


EIRP  =  12.566 


(W) 
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Radiation   Resistance   Rr: 


(Uoiy 


Rr  =  1.073 


Input    Resistance    R-jn: 

Rin:=Rr  (Q) 


Rin  =  1.073  (Q) 


Voltage  Rejection  Coeffecient  T: 

r ._  Rin  -  Zo  ,  .     .     - 

T  •= (dimensionless) 

Rin  +  Zo 


f  =-0.972  (dimensionless) 


Reflection  Efficiency  et  : 

Et  =ccd[l-(|r|)  J  (dimensionless; 


et  =0.056  (dimensionless) 
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Gain    G: 


G-  =  et  Do 


(dimensionless ) 


GdB-=  10  1og(G)    (dB) 


G  =  1.303 


(dimensionless ) 


GdB  =  1.149 


(dB) 


Maximum  Effective  Aperture  Apm: 
,2 


Aem  := 


GX 

4-71 


(m2) 


Aem  =0.104 


(mz) 


Maximum  Effective  Height  hemj 

.  „    /Rr  Aem 

hem  •=  2     (m) 


VP 


hem  =0.034 


(m) 
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TWO-DIMENSIONAL   ANTENNA    PATTERN:     (x-z)    plane 


Number   of    Elements    in   x-y   plane: 
Direction    of   Main    Beam: 

Xl(8):=  |AF(e,0)|  cos  8-  - 


N  =  10 

Go=0 

<(>o=0 


(radians) 
( radians) 


Y1(8):=|aF(8,0)|  sin  8-- 


8   =    0( 


xref    ,X1(9) 
m 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:   (y-z)  plane 


Number  of  Elements  in  x-y  plane 
Direction  of  Main  Beam: 


X2(0)  = 


9,- 


COS|U 

2 


N  =  10 
Go=0 
$o  =  0 

Y2(0) 

=    0° 


( radians ] 
( radians ] 


AFie,- 


«  — 


1 
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V"      ■     / 

\                                                                 ^' 

\                                                                           J^                                                            ' 

*'  \ 

1 

^                     '                    s                                , 
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'              '  %            "»      *  r\ >  \ 

r    /I  t           *                 -    »                    % 

'                                               '                 "    -.                         '           N                    ^     \    \' 

}'/ 1/         '    s             ^        *                    > 

-        /                N                   \     V  1 

-  -/y  /  w     v        ^      .  -  t 

p  jv  '  X x                V  *■ "           '                 * 

y^m 

;  "-AW 

'                                                      '      -  J  x    M 

9  J/ y     sr'                 .,  •»  "*   '                               i                              • 

m 

i                       ■                       •       y^  -»  ^  jtf  ^  ■*  ~\                              *•                       i 

1 vt;'A*:->"»L • t 1 

Y2(6) 

'                 ■         .- "  ~  T   v*"  #/ 

'                    '_-'"'            \V^    <// 1 

\  \\\       w  '                 —  -*  _.                            J                            ' 

i                       '                -  -r                   -X     '//  8 

\  VO     V  -                 i                      /                     ' 

\  \\/s'  \     x           '           *"  -  -'                   * 

1     -  *"  ~           '                            \  '              L^    <T  'JW 

1 1\  V                        '  "•                    ' 

1                 .          "       "     /       '//'I 

\,\  \k                /                     1 

'- 

,\J    i           »                       xi                        ' 

* -> 

'    "    1 

'    -  '                                         < 

'                          s                                  \ 

I                            '                                                   '                                                           " 

•'"^'-i   V 

/    \.>  \ 

-1 

1                                                                             '                                                      - ~    — '  —    —    ~                        ■                                                       ■                                                                              ' 
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xref    ,X2(6) 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:  &  =  n/4 


Number  of  Elements  in  x-y  plane:         N  = 10 
Direction  of  Main  Beam:  0o  =  O 

4>o=0 


( radians ) 
( radians) 


X3(0)  := 


AFie,- 


cos|« 


Y3(6):  = 


6,- 


H 


9   =    0° 


Y3(6) 


-1 
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-  '                                   <                        s                ' 
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/ 

'/9  I  \^« 

*i  -  l-  rv 
]      i      1  i 

'  ~  r  ~ 

\\A        x            ,            *■  -  _,                    ' 
\  v       N          ?  -                > 

1     VI         <•                                                                          XI                                                  ' 

'                                                                X                                                  '        ^ 
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-1 


xrcf    ,X3(9) 
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TWO-DIMENSIONAL  ANTENNA  PATTERN:  &  =  &( 


Number  of  Elements  in  x-y  plane: 
Direction  of  Main  Beam: 


X4(8):=  |AF(9,4o)|  cos|0-  - 


N  =  10 
0o  =  O 
4>o  =  0 


radians ) 
radians ) 


Y4(0):=|AF(0,(tK))|sin|0-- 


=    0< 


xref    ,X4(9) 
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FOLDED  DIPOLE 
MATHCAD  SOFTWARE  FOLDED. MCD 

To  provide  impedance  matching  to  coaxial  input  cables  and  maintain  good 
directional  pattern  characteristics,  a  folded  dipole  is  used  instead  of 
a  single  dipole  element.   The  geometry  of  a  folded  dipole  is  such  that 
the  spacing  between  dipole  elements  (d)  is  much  less  than  the  wavelength 
(X)  and  the  length  of  an  element  (L) . 

This  Mathcad  application  calculates  the  parameters  of  a  folded  dipole 
aligned  vertically  on  the  z-axis  and  a  folded  dipole  positioned 
horizontally  above  the  earth.   For  the  calculations,  the  antenna  is 
assumed  to  behave  as  a  cylindrical  dipole  with  equivalent  radius  (ae) . 

(Note:   Mathcad  equations  can  not  use  symbolic  subscripts.   Therefore, 
parameters  do  not  contain  subscripts  (i.e.,  Rj_n  will  be  written  Rin)  .  ) 

The  folded  dipole  antenna  Mathcad  application  will  compute  the  following 
parameters  (*  indicates  antenna  parameters  that  are  calculated  for  a 
folded  dipole  in  free  space  and  positioned  horizontally  over  the  earth) : 


X      =      Wavelength 

ae   =   Equivalent  Radius 

Rr   =   Radiation  Resistance 

R^   =   Input  Resistance  of  Linear  Dipole 

Xr   =   Reactance  of  Linear  Dipole 

X^   =   Input  Reactance  of  Linear  Dipole 

Z^   =   Input  Impedance  of  Linear  Dipole 

ZQ   =   Characteristic  Impedance  of  Two-wire  Transmission  Line 

Zt   =   Characteristic  Impedance  Transform  of  Folded  Dipole 

Z^n   =   Characteristic  Impedance  of  Folded  Dipole 

ZQ   =   Maximum  Current  of  Folded  Dipole 

rmin   =   Minimum  Distance  to  the  Far-Field 

E(9)   =   Electric  Field  * 

En(9)   =  Normalized  Electric  Field  * 

U(9)   =  Radiation  Intensity  * 

umax  ~      Maximum  Radiation  Intensity  * 

Pracj   =   Radiated  Power  * 

D0   =   Directivity  * 

EIRP   =   Effective  Isotropic  Radiated  Power  * 

Rin   =   Antenna  Input  Resistance  * 

T  =   Voltage  Reflection  Coefficient  * 
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Et   =   Antenna  Efficiency  * 

oa   =   Antenna  Polarization  Vector  * 

PLF   =   Polarization  Loss  Factor  * 

G   =   Gain  * 

Aem   =   Maximum  Effective  Aperture  * 

hem   =   Maximum  Effective  Height  * 

-rp 

Tv   =   Vertical  Reflection  Coefficient  of  Ground  Plane 

r^   =   Horizontal  Reflection  Coefficient  of  Ground  Plane 


The  following  known  or  estimated  data  must  be  entered: 

f  =  Frequency  of  Interest 

a  =  Radius  of  Conductor 

b  =  Radius  of  Conductor 

d  =  Distance  Between  Conductors 

L  =  Length  of  Antenna 

h  =  Height  of  Antenna  Above  Earth 

eccj  =   Conduction/Dielectric  Efficiency 

er   =   Relative  Dielectic  Constant  of  Ground 

o  =   Conductivity  of  Ground 

rff   =   Far  Field  Distance 

i   =   Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 

ow   =   Incoming  Wave  Electric  Field  Unit  Vector  For  Free  Space  Antenna 

°wh   =   INcoming  Wave  Electric  Field  Unit  Vector  For  Horizontal  Antenna 
9ph   =   Antenna  Polarization  Direction  For  Horizontal  Antenna 
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Enter  input  data  here 


aw 


0 


[dimensionless) 


awh  := 


[dimensionless) 


f :  =  3  1 0° 


(Hz) 


Gph  :=0  (radians  ] 


a  :=  0.0005 


(meters ) 


rfif  :=  110  (meters) 


b  :=a 


(meters ) 


l  :=90 


( radians 


d:=  12.5  a 


(meters) 


ecd:=  1 


(dimensionless) 


L:=0.5 


(meters ) 


o:=25 


( siemens/meter ) 


h:=5 


(meters) 


er:=  15 


(dimensionless; 
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Define  constants 


oa  •= 


(dimensionless) 


c  :=  2.9979  10  (meters/sec; 


Tio  :=  1 20- n       <Q) 


l=Z 


(meters/cycle; 


k.=i_" 


(m-1) 


X  =0.999 


(meters/cycle; 


•y  :=  0.57721  (dimensionless) 


eo-=—  10" 9 

36n 


(F/m) 


z-  =  0  .2  i 


( radians ) 


Define  angular  offset  8: 


•  =  0, —  ..  2  Tt 


radians; 


Define  cosine  [Ci(x)]  and  sine  [Si(x)]  integrals 


n  :=  1 ..  50 


Ci(x)  -=0.57721  + 


ln(x)-]T 


(■l)n-'x2'n 
(2-n)(2n)! 


Si(x)-=2 


(-1)-1^2-B-I) 

(2n-  l)(2n-  1 )! 
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FOLDED  DIPOLE  ALIGNED  ON  THE  Z-AXIS  IN  FREE  SPACE 


Equivalent  radius  ap: 


ae  =  exp 


a2ln(a)  +  b2ln(b)+2  abln(d) 


(a+b)' 


(meters) 


ae  =0.002 


(meters! 


Electric  field  E  (0) 


E(9)  =j   -tjo- 


-j    krff 

e  J 
2-n-rff 


/       /kL  A  /k-L\\ 

'cos cos(0)    -COS  —  ■ 

2  I  \2 


\  sin(9) 


(V/m) 


En(e)  :=S8i 


- 


(V/m) 


Minimum  Distance  to  Far-Field  r 


mm  j. 


r0:=1.6-i 


(meters; 


r   :=5  L  (meters) 


V 


2(L)V 


(meters) 


rmin  :=  max( r)  (meters) 


rmin  =  2.5 


(meters) 
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Radiation  Intensity  U(6) 

2-TJO 


rv2 

U(6):  =  —  (|E(e)|)2 


(W/solid  angle) 


Radiated  Power  Prarj: 


Prad  :=  2- 71 


U(8)-sin(8)d9 


0 


(W) 


Prad  =  36.64 


(W) 


Directivity  D( 


Umax:  =  U|- 
[2, 


(W/solid  angle) 


Do:= 


4  71  Umax 
Prad 


(dimensionless 


Do  =  1.641 


(dimensionless; 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :=  Prad  Do 


(W) 


EIRP  =60.132 


(W) 
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Radiation  Resistance  Rr: 


Rr:=^ 
2-n 


y+  ln(kL)-  Ci(k  L)  +  -  sin(kL)  (Si(2kL)-  2  Si(kL)) 
2 

+  -  cos(kL)  h+lnf—  |  - 


(Q) 


Rr  =  73.281 


(0) 


Input    Radiation    Resistance    of    linear    dipole    R,-j : 
Rr 


Rd:= 


sin 


k-L 


(Q) 


Rd=  73.281 


(Q) 


Reactance    of    linear   dipole    Xr: 


Xr:  = 


_no  r 

4n 


2  Si(kL)+  cosfk  L)  (2  Si(kL)-  Si(2 -k-L)) 


1  sin(kL)    2  Ci(kL)-  Ci(2kL)-  Ci 


2  k  &c 


(Q) 


Xr  =  43.142 


(Q) 


Input  reactance  of  linear  dipole  X^ 


Xd:=. 


Xr 


m) 


(Q) 


Xd  =43.142 


(Q> 
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Input  impedance  of  linear  dipole  Z^: 

Zd  =Rd+j  Xd  (Q) 

Zd  =  73.281  +43.1421  (Q) 

Characteristic    impedance   of    two-wire    transmission  line    Zn: 

Zo:  =  ^acoshf — — |  (Q) 

Zo  =  302.312  (Q) 


Characteristic  impedance  transform  for  folded  dipole  Zt : 
Zt:=j  Zotanf— j  (Q) 


Zt  =-2.747  105i  (Q) 


Characteristic  impedance  of  folded  dipole  Z^n: 

4ZtZd 
Zin:  = (Q) 

2Zd+Zt 


Zin  =293.307 +172.465i  (Q) 

Input  Resistance  Rjn: 

Rin   =Re(Zin)  (Q)  Rin=293.307  (Q) 
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Voltage  Rejection  Coeffecient  £: 
Rin  -  Zo 


r- 


Rin  +  Zo 


(dimensionless) 


r  =  -0.015 


(dimensionless) 


Reflection  Efficiency  gt  : 

et  :=ecd  [  1  -  ( |r|  )  J   (dimensionless) 


et  =  l 


(dimensionless) 


Gain  G: 


G:  =  et  Do 


[dimensionless) 


GdB:=  10  1og(G)    (dB) 


G  =  1.641 


(dimensionless ) 


GdB=2.15 


[dB) 


Polarization  loss  factor  PLF: 


IV  2 


PLF 


:=  \|ow  aa\) 


(dimensionless; 


PLF  = 


[dimensionless) 


Maximum  Effective  Aperture  Aprn: 

C  \2 
Aem  :=-^-- PLF     (m2) 

4n 


Aem=0.13 


(m2) 


Maximum   Effective    Height    hPni: 

,        .    «    fRrAem         ,    , 
hem  =  2     (m) 


VP 


hem  =0.318 


[m) 


138 


Far  Field  radiation  pattern  of  vertical  folded  dipole  in  free  space . 


X(0)  :=  |En(0)|-cos|e+- 


Y(8):=|En(0)|  sin  6+- 


=    0° 
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xref    ,X(9) 
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HORIZONTALLY  POSITIONED  FOLDED  DIPOLE  OVER  THE  EARTH: 


Define  angular  offsets  01  and  02: 


61  •=-  — ,-— t  —  ■•—      (radians 
2  2   2i  2 


02  := h z   (radians) 

2   2  i 


Relative  complex  effective  dielectric  constant  err,: 


op  :=  er-  j    - 


2  7t  feo 


(dimensionless) 


ji   component    of    E-field   Eji: 


E<t>(91):=j    no- 


-I    k  iff 

e  J 
2  7i  rff 


cosl — sin(01 )  sin(0)j  -  cosf — 
2  \2 


JT- 


sin(61)2sin(0)2 


(V/m) 


E<t»l    :=j    no- 


,-J    k  rff 

2ti  rff 


cosf sin (02  J  sin(0)|  -  cosf — 


Jl-  sm(0232sin(O)2 


(V/m) 


6   component    of    E-field   E0: 


E9(01)  :  =  j    no. 


,-J 


krff 


2ti  rff 


COS 

It 

sin(91 )  sin 

ill 

— 

:os 

It) 

J 

2 

1-  sin(01) 

sin 

(-5 

(V/m) 


E01    :=j    no- 


-j    krff 

2nrff 


cosf sin  f  02  J  sin]  —  j  1  -  cosf- 


k_L 

2 


'W 


l  I    sinf— 

L2, 


(V/m) 
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Vertical  plane  wave  reflection  coefficient  Yx 


rv(91 )  •-gPcos(ei)-  ^op-  sin(61) 


-  -^e 


■J-'Ji 


op  cos(01 )  4-  V,erp-  sin(91 ) 


[dimensionless) 


rvi 


op  cos 


(92j-Jcrp-s,n(e22): 


erp  cosf 02 J  +  Jap-  sin f 02 J' 


(dimensionless 


Horizontal  plane  wave  reflection  coefficient  T^ : 


nxeD-.i^l) 


-  ^erp-  s 


sin(01) 


4^erp- 


cos(01  )4  Jap-  sin(01 ) 


(dimensionless ) 


cos 


rhi   =■ 

z 


cos 


(92j-jap-sin(82j: 
(02  j  4  Jap  -5111(021 


(dimensionless) 


E-field    (j>   component)    with    ground   plane    reflection    Ejjgp : 
E^gp(01).  =  E^(01)  (e>   khcos(ei)+rh(01)e-J   khcos<ei>)  (V/m) 


E<t»gpl    '■- 


E<|>1 


/  j    k  h  cos(e2^ 


fhl   e 


-j    khcos(92^1 


(V/m) 


EHei).=  |Et8P<9l)l 


max(  E()>gp  1 ) 


(V/m) 
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E-field  (6  component)  with  ground  plane  reflection  EO^p: 

Eegpon^Eeon-P  khCM(ei)-rv(ei)e-j  khcos<91>)  (v/m) 


EGgpl. 


/  j    k  h  cosfo}  -j  ■k-h-cos(92^1 

E81    le  -  Tvl    e 


(V/m) 


EMei):=|Eea,(9l)l 


max(  E0gp  1 ) 


(V/m) 


Radiation  Intensity  Uh(8l): 

uh(9i )  =— ■[(  Ngpcei  )|  )2+  ( lEGgpiei  >|  )2 

2  rp 


(W/solid  angle) 


Radiated  Power  Phrarj: 


Phrad:=2  n 


Uh(81)  sin(91)d61 


(W) 


Phrad  =  70.415 


(W) 


Directivity   Phi-,: 
rfl2 


Ugp_max  := 


2  rp 


(|l*gpTJ)  +(|iWJ) 


(W/solid   angle) 


Uhmax  =  max(Ugpmax) 


(W/solid   angle) 
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Directivity   Dhn: 
4-71-Uhmax 


Dho  =. 


Phrad 


(dimensionless) 


Dho -3.938 


(dimensionless; 


Effective  Isotropic  Radiated  Power  EIRPh: 


EIRPh  :=  Phrad  Dho 


(W) 


EIRPh  =277.3 


(W) 


Input    Resistance    Rh.,n  : 
2  Phrad 


Rhin:=- 


sin 


k-L 


Rhin  =  140  83 


(Q) 


Voltage    Rejection    Coeffecient    TH: 
Rhin-  Zo 


m:= 


Rhin  +  Zo 


(dimensionless) 


TH  =-0.364 


(dimensionless ; 


Reflection  Efficiency  eht : 

eht  :  =  ecd  [  1  -  ( |rH|  )  J     (dimensionless) 


eht  =0.867     (dimensionless) 
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Gain    Gh: 


Gh  :=EhtDho 


( dimensionl ess) 


GhdB  :=101og(Gh)  (dB) 


Gh  =3415 


(dimensionless) 


GhdB  =5.334 


(dB) 


Polarization    loss    factor    PLFv: 


Ehmax=|E9n(9ph)|  (V/m) 


Lobes  -floor I —  +  1 


aah  :  = 


0 

0 

Ehmax 


Lobes 


(dimensionless; 


PLFh:=(l<whoah|) 


(dimensionless) 


PLFh  =  7.736  10 


(dimensionless! 


Maximum  Effective  Aperture  Avpm: 
,2 


Ahem  :=  "LL-PLFr 
4-ti 


(m2) 


Ahem=2  099-10~4(m2) 


Maximum  Effective    Height    hvprn: 


hhem  :=  2 


iRhin  Ahem 


m 


no 


hhem  =0018  (m) 
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H-Plane  Radiation  Pattern  (y-z  plane)  of  Horizontal  Folded  Dipole  Over 
The  Earth: 


Yl(61)  =  |E<t>n(91)|  sinfei  +- 


Xl(91)  =  |E(J.n(01  )|  cosfei  +  - 


yrcf 
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E-Plane  Radiation  Pattern  (x-z  plane)  of  Horizontal  Folded  Dipole  Over 
The  Earth: 


Y2(01)=  |E9n(01)|  sin  61  +- 


X2(61)  =  |E9n(81)|  cos[ei  +  - 


J™*'- 


0.5   - 


Y2(9l) 


xref    ,X2(9l) 
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CAGED  DIPOLE  ANTENNA 
MATHCAD  SOFTWARE  CAGED_DI.MCD 

This  caged  dipole  application  is  modeled  after  a  circular  array  with  the 
radiating  elements  as  dipoles.   The  dipoles  are  modeled  after  conductors 
with  infinitely  small  radii  and  antenna  length  greater  than  X/4.    Mutual 
interference  between  adjacent  elements  are  neglected  when  calculating 
antenna  parameters. 

(Note:   Mathcad  equations  can  not  use  symbolic  subscripts.   Therefore, 
parameters  do  not  contain  subscripts  (i.e.,  AFn  will  be  written  AFn) . ) 


The  caged  dipole  antenna  Mathcad  application  will  compute  the  following 
parameters : 

VERTICAL  CAGED  DIPOLE  IN  FREE  SPACE: 


X      -      Wavelength 

rmin   =   Minimum  Distance  to  the  Far-Field 

AF(0,4>)  =   Normalized  Array  Factor 

E  {6,4»)   -   E-Field  of  Finite  Length  Dipole 

Et(0>)   =   Total  E-Field 

Etn  (0,<|))   =   Total  Normalized  E-Field 

U(0,i(i)   =   Radiation  Intensity 
Pracj   =   Radiated  Power 

D0   =   Directivity 

EIRP   =   Effective  Isotropic  Radiated  Power 

Rr   =   Radiation  Resistance 

^in   =   Input  Resistance 

T   =   Voltage  Reflection  Coefficient 

et   =  Antenna  Effeciency 

G   =   Gain 

PLF   =   Polarization  Loss  Factor 

Aem   =   Maximum  Effective  Aperture 

nem   =   Maximum  Effective  Height 
flow   =   Low  Frequency 
fhigh   =   High  Frequency 
BW   =   Bandwidth 

VERTICAL  CAGED  DIPOLE  OVER  EARTH: 


erp   =   Relative  Complex  Effective  Dielectric  Constant 
AFv(01,$)  =   Normalized  Array  Factor 
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Ev(01,<j>)   =   E-Field  of  Finite  Length  Dipole 

rv(01)   =   Vertical  Plane  Wave  Reflection  Coefficient 

Evt(01,<Ji)   =   Total  E-Field 

Evtn(8l,ij>)   =   Total  Normalized  E-Field 

Uv(01,$)   =   Radiation  Intensity 
Pvracj   =   Radiated  Power 

Dv0   =   Directivity 

EIRPv   =   Effective  Isotropic  Radiated  Power 

Rvr   =   Radiation  Resistance 

Rv^n   =   Input  Resistance 

rv   =   Voltage  Reflection  Coefficient 
ev^   =   Antenna  Efficiency 

Gv   =   Gain 

PLFv   =   Polarization  Loss  Factor 

Avem   =   Maximum  Effective  Aperture 

hvem   =   Maximum  Effective  Height 
HORIZONTAL  CAGED  DIPOLE  OVER  EARTH: 


AFh(01,<|tl)  =   Normalized  Array  Factor 

T^Ol)   =   Horizontal  Plane  Wave  Reflection  Coefficient 

E9h(Gl)   =   6  Component  of  E-Field 

E8ht(ei,$l)   =   G  Component  of  Total  E-Field 

E0htn  (01,  (jil )   =   0  Component  of  Total  Normalized  E-Field 

E$h(01)   =   <ji  Component  of  E-Field 

E$ht (01,^1)   =   4  Component  of  Total  E-Field 

E<|>htn  (01 ,  4>1 )   =   4  Component  of  Total  Normalized  E-Field 

Uh(01,(jil)   =   Radiation  Intensity 
Phrad   =   Radiated  Power 

Dh0   =   Directivity 

EIRPh   =   Effective  Isotropic  Radiated  Power 
Rhr   =   Radiation  Resistance 

Rhj_n   =   Input  Resistance 

TH   =   Voltage  Reflection  Coefficient 
eht   =   Antenna  Efficiency 

Gh   =   Gain 

PLFh   =   Polarization  Loss  Factor 

Ahem   =   Maximum  Effective  Aperture 

hhem   =   Maximum  Effective  Height 
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The  following  known  or  estimated  data  must  be  entered: 


o   = 


ov. 


9vF 
ohv 
6hr 


=   Number  of  Dipoles 
=   Radius  of  Antenna 

Length  of  Antenna 

Antenna  Height  Above  Ground 

Frequency  of  Interest 
=   Antenna  Feed  Current 

=   Input  Impedance 

=   Relative  Dielectric  Constant  of  Ground 

Conductivity  of  Ground 
=   Incoming  Wave  E-Field  Unit  Vector  for  Antenna  (Free  Space) 

=   Antenna  Unit  Polarization  Vector  (Free  Space) 
=   Incoming  Wave  E-Field  Unit  Vector  for  Vertical  Antenna 
=   Direction  of  Incoming  Unit  Vector  (Vertical  Antenna) 
=   Incoming  Wave  E-Field  Unit  Vector  for  Horizontal  Antenna 
=   Direction  of  Incoming  Unit  Vector  (Horizontal  Antenna) 
Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 


ow  := 


Enter  input  data  here: 

I, 


[dimensionless ) 


oa 


[dimensionless ) 


ovw 


(dimensionless ) 


.       _  80-Tt 

9vp  = 

180 


[dimensionless) 


ohw  := 


[dimensionless ) 


0hp:  =  O 


[dimensionless ) 


N  =8 


(elements ) 


Io  =  1.0 


(A) 


a  :=  — 
8 


(meters ) 


Zo:  =  75.0 


:q> 
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L:=0.5 


(meters) 


j  :=90 


(increments) 


h-=0.25 


(meters; 


o  ■=  10' 


(siemens/m) 


f:=3  10° 


(Hz) 


er:=15 


[dimensionless) 


Define  constants  and  calculate  wavelength 


c:=  2.9979- 10 


(meters/sec) 


tio:=  120ti 


(Q) 


l:=< 


(meters/cycle) 


k-22 


(nT1) 


i     -Io 
In  •=  — 

N 


(A) 


rff.=  10J 


(meters) 


eo:=_L  10*9 
36ti 


f arads/m) 


ecd  :=  1 


(dimensionless 


a.  =0.999 


(meters/cycle) 


z=0  2j 


( increments) 


Calculate  caged  dipole  parameters  parameters: 


Define  angular  offset  6_: 


=  0, —  .  2  7i 


( radians ] 


ty  ■=  0, — ..  2  n  (radians; 

J 
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Minimum  Distance  to  Far-Field  r 


rain: 


r0:=1.6X 


(meters ) 


r  :=5-^4 


(2a)    (meters! 


.  2-[L2+(2-a)2]   .   . 

r  .  =  — ' - (meter 


s) 


rmin  :=  max(r)       (meters) 


rmin=  2.795 


(meters ) 


VERTICAL  CAGED  DIPOLE  IN  FREE  SPACE: 


Array  Factor  AF(8): 


n:=  l.N 


(elements ] 


<1>   :=  2  7t  — 

N 


( radians ) 


AF(8,t)  :=Vln  e 


j  -k-a-sin(9)-cos($-  <I>J 


(dimensionless; 


E-Field  of  finite  dipole  E  (9) 


E(G)  :=j    rp 


In  e' 


j   krff 


2-n-rff 


cos cos(9)   -  cos  — 

\  sin(e)  / 


;v/m) 
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Total  E-Field  Et  (9)  &  Normalized  E-Field  Etn_(e) 


Et(8,<t>)  :=E(0)  AF(0,<fr) 


(V/m) 


Etn(9,$):  =  - 


Et(G,<J>) 


ti?° 


(V/m) 


Radiation  Intensity  U  (8) 

rff2 
U(6,(fr).=  — (|Et(9,,|>)|)2 

2-rp 


(W/solid  angle; 


Radiated  Power  P 


radl 


Prad  := 


2  7t 


U(9,<|0  sin(9)d8d$ 


(W) 


Prad  =0.441 


(W) 


Directivity  Dn 


Umax  :=U|-,0|    (W/solid  angle) 


Do 


4  7i- Umax 
Prad 


(dimensionless) 


Do  =  1  543 


;  dimensionless ) 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP  :=  Prad  Do 


(w; 


EIRP  =0.681 


(W) 
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Radiation   Resistance    R, 


Rr:  = 


2  Prad 

(N? 


:Q) 


Rr  =  56.487 


(0) 


Input    Resistance    R1n : 
Rr 


Rin:  = 


(0) 


sin 


kL1 


Rin  -  56.487 


(Q) 


Voltage  Rejection  Coeffecient  T: 
Rm  -  Zo 


r:= 


Rin  +  Zo 


( dimension! ess ] 


r=-o  hi 


(dimensionless ) 


Reflection  Efficiency  gt  : 


et  '=  ecd 


[i-(|r| 


(dimensionless) 


et  =  0.98 


(dimensionless ) 


Gain   G: 


G:  =  et  Do 


(dimensionless) 


GdB:=  10  1og(G)         (dB] 


G  =  1.513 


(dimensionless) 


GdB  =  1.798 


(dB) 


Polarization    Loss    Factor    PLF: 


PLF  =  flow  oa|) 


(dimensionless) 


PLF  =  1 


(dimensionless) 
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Maximum  Effective  Aperture  A^m : 

r  i2 

Aem:=-^PLF     (m2) 

4  71 


Aem  =0  12 


(m2) 


Maximum  Effective  Height  hpm: 


,     „  fRr  Aem 
hem:=2-  (m) 

J   T|0 


hem  =  0.268 


(m) 


Define  Frequency  Relationship  F(x) 


-5 


F(x)  =-5.696  10"  Jx+ 0.31481 


(dimensionless ) 


Frequency   Range    f-|  nvi   and    fnjghJ 


flow:  =  f-F|_]f        (Hz) 
l2- 


flow  =2.056  108  (Hz) 


:f+F|— J    f  (HZ 


thigh  =3  944  10*  (Hz) 


Bandwidth    BW: 


BW  =  fhigh-  flow         (Hz; 


BW  =  1.888- 108  (Hz) 
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E-FIELD  PATTERN:   This  far-field  radiation  pattern  is  for  the  vertical 
caged  dipole  in  free  space: 


Number  of  Dipole  Conductors: 


N  =8 


Radius    of   Antenna: 


a  =0  125 


(meters ) 


X(9).=  |Etn(8,*)|  cosfe--] 


Y(0)  ■=  |Etn(e,4.)|  sinle-  - 


e  =  oc 
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VERTICAL  CAGED  DIPOLE  OVER  EARTH: 


Define  angular  offset  9_1  and  Q2 

( radians ) 


01  Z.1..1  +  J- 


2  2   2-j  2 


*2  =  -z 
z  J 


(  radians ) 


62  = 1- z 


radians 


2   2-j 


Relative  complex  effective  dielectric  constant  c 


£E- 


0"p  -  ET-  J 


2  n  feo 


(dimensionless ) 


Array  Factor  AFv  (8)  : 

El    k  asin(Gl)  cos(<t>-0.) 
Ine  V       * 


(dimensionless) 


APvgp    =  >    In  e 


(dimensionless) 


E-Field   of    finite   dipole    Ev(9l): 


Ev(91)  =j    Tp 


In  e 


j    k  iff 


/       /kL  \  /kL\\ 

cos cos(91)   -  cos  — 

2         y       \  2 


2-n-rff 


1  sin(91)  / 


(V/m) 


Evgp   :  =  j    ip 


In  e 


-j    krff 


2  TirflF 


cos 


±  COS(02)J  -  COS^ 


sin(e2j 


(V/m) 


Vertical  plane  wave  reflection  coefficient  Ev(91): 


rv(01):_ap_cos(01)- Vop-  sin(Gl ) 


i-  «Je 


H-'Jl 


erp  cos(91 )-»-  ^erp-  sm(91 ) 


(dimensionless ) 


156 


rvl   :=■ 

z 


erp-  cos 


(62 }  -  Jap-  s,n(e2j: 


erp  cosf  02j  +  Jap-  sin  f  02  V 


(dimension! ess  \ 


Total    E-Field    Evt(01)     &    Normalized    E-Field    Evtn (01 
Evt(61,t)  =Ev(ei)-AFv(ei,*)-(ej   kh  cos<ei>4  fv(01  KJ   khcos<ei>) 


(V/m) 


Evtgpz  := 


j    khcos(92j  -j    k  h  cos(02^ 

Evgp   Al-vgp    le  s   *  +-  fvl    e 


z      z 


(V/m) 


max(  Evlgp ) 


(V/m) 


Radiation  Intensity  Uv(0j) 

rfT2 
Uv(0l,$):=  — —  ( |Evt(0l  ,4>)| ) 

2  T|0 


(W/solid  angle) 


Uvlz:  =  _(|i^)|)2 

2  T]0 


(W/solid    angle) 


Radiated    Power    Pv 


radl 


Pvrad  := 


2-n 


Uv(8l,t)  sin(0l)d0l  cfy 


(W) 


Pvrad  =  0.483 


(W) 


Directivity   Dvn : 
Uvmax  :=  U v(  9vp ,  0 ) 


(W/solid   angle) 
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Directivity   Dvn 
4-7iUvmax 


Dvo:= 


Pvrad 


(dimension! ess) 


Dvo  =  4.415 


(dimensionless) 


Effective  Isotropic  Radiated  Power  EIRPv: 


EIRPv  :=  Pvrad  Dvo       (W) 


EIRPv  =2.133  (W) 


Radiation    Resistance    Rv, 


Rvr  := 


2  Pvrad 
(Un|)2 


(Q) 


Rvr  =61.857 


(Q) 


Input    Resistance    Rv.jn: 
Rvr 


Rvin  '=• 


(Q) 


sin 


k-L 


Rvin  =61.857 


(Q) 


Voltage  Rejection  Coeffecient  [V: 
Rvin-  Zo 


TV:=- 


Rvin  +-  Zo 


(dimensionless 


TV  =-0.096 


(dimensionless) 


Reflection  Efficiency  ev^  : 

cvt  :  =  ecd[  1  -  (  |rv|  )  J    (dimensionless) 


evt  =0.991 


(dimensionless) 
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Gain    Gv : 


Gv  =evtDvo 


( dimension! ess ) 


GvdB  =  lOlog(Gv)     (dB) 


Gv  =4.374 


(dimensionless ) 


GvdB  =  6.409 


[dB] 


Polarization    Loss    Factor    PLFv: 


Evma\:  =  Evtn(0vp,O)     (V/m) 


lobesv  :=  floor  [  2  —  +  1 


ova  •  = 


Evmax 


lobesv 
0 
0 


(dimensionless ) 


PLFv  =  V|ow  <jva|J         (dimensionl 


ess ) 


PLFv  =0.986 


(dimensionless ) 


Maximum  Effective   Aperture   AOTT1: 
Avem  := PLFv      ( m2  ) 

4-71 


Avem  =0.343 


(mz) 


Maximum  Effective    Height    riv,^: 


,              „    IRvr  Avem 
hvem  :  =  2    f - ( m ) 


hvem  =0  474 


(m) 
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E-FIELD  PATTERN:   This  far-field  radiation  pattern  is  for  the  vertical 
caged  dipole  over  earth: 


Number  of  Dipole  Conductors 


N  =8 


(elements ) 


Radius  of  Antenna 


a  =0.125       (meters ) 


Height  of  Antenna  Above  Earth 


h  =  0.25 


(meters) 


X1(G1)  :=  |Evtn(81,4i)|  cos  61  +  -  Yl(81)  :=  | Evtn( 91 , 4> ) |  sinl 01+-—] 


=    no 


0   =    0 


vef„ 


Yl(Gl) 


xref    ,Xl(9l) 
m 
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HORIZONTAL  CAGED  PI POLE  OVER  EARTH: 

Define  angular  offset  |1  and  J2 : 

(jil  :=0,-..2ti       (radians) 
J 


Array  Factor  AFh(8l,$l); 


n  '=  1 ..  N 


(elements ) 


0   -=2n- 


N 


( radians) 


AFh(01  ,$!)•=  Vln 


j     k  a  (sin(  4>  1  >  sin(Gl)  sin(©1)  4-cos(  9l )  cos(©J) 


(dimensionless 


AFhl 


=.£lnV 


k  a  (sin(i|>2^   sin(92^   sin(©,)  -+cos(e2^   cos(©,)) 


(dimensionless) 


component  of  E-Field  of  finite  dipole  E0h(01) 


E0h(01)  =j    no 


In  e 


j    krff 


2-n-rff 


ikL       ,fln  kV 

cos sin(91)    -  cos  — 

2  \2 


f- 


sin(61) 


E6hl    =j    no 


In  e 


j  krff 


2  7i  iff 


"k-L 

cos  | sin 

2 


M)-«  (t) 


Jl-sin(82j: 


(V/m) 


(V/m) 


Total   E-Field   E0hf  (61,  &1 )  : 

E9ht(ei,*l):=E0h(ei)  AFh(01>l)  -{j   khcos(9l)_  rv(01)e"J   khcos(9l))  (V/m) 


E0htl    := 

z 


E8h 


lI.AIl.lI-^*t-W.Ivli.e-J*t-«) 


(V/m) 
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Normalized   E-Field    Eghtn  (81 ,  $1 ) 
|E9ht(ei,4>l)| 


E6htn(91,*l)  := 


max(E0htl) 


(V/m) 


Horizontal  plane  wave  reflection  coefficient  Eh  (91): 


rh(91): 


•  "Jg?- 


._  cos(61 )-  ^erp-  sin(81 ) 


i-f  fijap- 


cos(61 )  +  Jap-  sin(91 ) 


(dimensionless ) 


Fhl   :=■ 

z 


cos (02]  4  Jap-  sinfG2_) 


(dimensionless) 


jl   component    of    E-Field   of    finite   dipole   Ej>h(9l) 


E4>h(91)  =j    r]0 


In  e" 


2-7i-rff 


,k-L 
1  -  cosl 


(V/m) 


„.,      .             Ine"J  "'"      /.  /k  L' 

h-4»hl     =j    tp 1  -  cos  — 

2-n-rff      \  \  2 


(V/m) 


Total    E-Field   E$ht  (61, $1 )     &    Normalized   E-Field   E$ht-n  (91,  $1)  : 
E^ht(01  ,*!):=  Eih(ei  )-AFh(91  ,*1  )■  (e*   k  h  C0S(9lV  rh(61 )  e'J  khcos(9l>)  (V/m) 


E4htl    := 


E<j>hl    AFhl 

z  z 


(j    khcosfo) 
e  *  4 


rhl    e 

z 


-j    kh  cos(e2^ 


)] 


(V/m) 


E4hwei.ti)»J*«!!Lili 

max(E(tthtl) 


(V/m) 
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Radiation  Intensity  Uh(8l,j>l): 


Uh(91 ,41 )  ■=—•[(  |E9ht(0]  ,<(.!  )|  )24  ( |E<(.ht(01  ,<frl  )|  )2] 


2  T|0 


(W/solid  angle) 


Uhl  :  = 


2  T]0 


*IS2 


(|E0htl|)  +  (|E$htl|) 


(W/solid  angle; 


Radiated  Power  Pv 


rad^ 


■2* 

Phrad  •= 

■ 

0    J 

Uh(01  Jl)  sin(91)d01  d<t»  1 


(W) 


Phrad  =  1.059 


(W) 


Directivity  Dh^ : 


Uhmax-=Uh(6hp,0)   (W/solid  angle; 


n,   _4-ji-Uhmax  , 

Dho  •= (dimensionless) 

Phrad 


Dho  =  5.051 


(dimensionless ; 


Effective  Isotropic  Radiated  Power  EIRPh: 


EIRPh  :=PvradDvo       (W) 


EIRPh  =2.133     (W) 


Radiation  Resistance  Rv 


r_L 


Rhr  •= 


2  Phrad 
(|ln|)2 


(Q) 


Rhr  =  135.58 


(Q) 
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Input    Resistance    Rh-jn: 
Rhr 


Rhin  :=■ 


(Q) 


sin 


k-L 


Rhin  -  135.58  (Q) 


Voltage  Rejection  Coeffecient  fH 
Rhin-  Zo 


TH:= 


Rhin  4-  Zo 


( dimension! ess ] 


TH  =0.288 


(dimensionless ) 


Reflection  Efficiency  cVf  : 

eht  :  =  ecd[  1  -  (  |rH|  )  J   (dimensionless) 


eht  =0.917 


(dimensionless ) 


Gain    Gh 


Gh  =eht-Dho 


(dimensionless ) 


GhdB  =  lOlog(Gh)     (dB) 


Gh  =4.633 


(dimensionless ) 


GhdB  =  6.659 


(dB) 


Polarization    Loss    Factor    PLFh: 

Ehmax:=E<t»htn(ehp,0)     (V/m) 

0 

0 

Ehmax 


cha  :  = 


lobesh 


(dimensionless  \ 


PLFh  =(|ohwoha|)2   <dimensionl 


ess) 


lobesh  :=  floor f 2—  +  1 


PLFh  = 


(dimensionless ) 
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Maximum  Effective  Aperture  Ah,^: 

2 
Ahem  :=-^-L  PLFh   (m2) 

4-71 


Ahem  =35.536 


(m2) 


Maximum  Effective  Height  hh^m: 


„  [Rhr  Ahem 
hhem  :=  2  f - (m) 

J    no 


hhem  =  7  15 


(m) 
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E -PLANE  RADIATION  PATTERN:   This  far-field  radiation  pattern  in  the 
(x-z)  plane  is  for  the  horizontal  caged  dipole  over  earth: 


Number  of  Dipole  Conductors 


N  =8 


(elements; 


Radius  of  Antenna 


a  =0.125       (meters) 


Height  of  Antenna  Above  Earth: 


h=0.25 


(meters ) 


X3(01):= 


E8htn  01, 


cos|81  +- 


Y3(01):= 


EGhtnfei,-) 


■sin|ei4- 


=    0° 


yreL 


Y3(9l) 


xref    ,X3(9l) 
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H-PLANE  RADIATION  PATTERN:   This  far-field  radiation  pattern  in  the 
(y-z)  plane  is  for  the  horizontal  caged  dipole  over  earth: 


Number  of  Dipole  Conductors 


N  =8 


(elements) 


Radius  of  Antenna 


a  =0.125       (meters ) 


Height  of  Antenna  Above  Earth 


h=0.25 


(meters) 


X2(81 )  :=  |E<(>htn(ei  ,0)|  cosfei  t  - 


Y2(01)  :=  |E^htn(01,O)|  sin|01  +- 


=    0° 


yrefn 


Y2(91) 


xref   ,X2(91) 
m 
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THE  HELICAL  ANTENNA  (REFLECTOR  OPTION) 
MATHCAD  SOFTWARE-DISH_HEL.MCD 
When  built  to  the  proper  specifications,  the  helical  antenna  possesses 
many  qualities  which  make  it  suitable  for  a  wide  variety  of 
communications  applications.   If  the  following  conditions  are  satisfied 
the  helix  will  exhibit  a  highly  directional  axial  main  lobe,  low  side 
lobe  level,  negligible  mutual  interference  with  adjacent  antennas,  low 
voltage  standing  wave  ratio  (VSWR) ,  and  resistive  input  impedance  over  a 
wide  frequency  band: 

.8   <   Cjl  <  1.15 

n  >  3 
12   <   a  <   14 

(Note:  X   in  a  subscript  indicates  the  dimension  is  in  wavelengths. 
Mathcad  equations  can  not  use  symbolic  subscripts.   Therefore,  the  symbol 

X   will  immediately  follow  the  parameter  in  equations  (i.e.,  CX)   in  lieu 
of  subscripts.) 

The  helical  antenna  Mathcad  application  will  compute  the  following 
parameters  (Items  with  *  indicate  parameters  that  are  calculated  for 
both  axial  and  peripheral  feed  geometries) : 

C   =  Circumference  of  Helix 

X      -      Wavelength 

a   =   Pitch  Angle 

Do  =   Directivity 

p   =   Relative  Phase  Velocity 

Eg  a   =   Electric  Field  Components 

¥   =   Array  Factor  Phase  Shift 

U   =   Radiation  Intensity 

Prad   =   Radiated  Power 

R   =   Antenna  Input  Resistance* 

r  =  Voltage  Reflection  Coefficient* 

er   =   Reflection  Efficiency* 

hem   =   Maximum  Effective  Height* 

G   =   Gain* 

EIRP   =   Effective  Isotropic  Radiated  Power* 

Aem   =   Maximum  Effective  Aperture* 

AR   =  Axial  Ratio 

PLF   =   Polarization  Loss  Factor 

BW  =  Bandwidth 

^high  =  uPPer  Frequency  Limit 

fj_ow  =  Lower  Frequency  Limit 

Acceptable  Conductor  Diameter 

Ev  ,,  -,   =   Electric  Field  Cartesian  Coordinates 

A/  y ,  z. 
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P' 


=   Unit  Polarization  Vector  Coordinate  Angles 


oa   =   Antenna  Unit  Polarization  Vector 


-mm 


=   Mininum  Distance  to  the  Far-Field 


The  following  data  must  be  input  based  on  known  or  estimated  data 


D  = 

S  = 

L  = 

n  = 

d  = 

f  = 

m  = 

i  = 

lo 

zo 
°w 


Diameter  of  Helix  (Center  to  Center) 

Spacing  Between  Turns  (Center  to  Center) 

Length  Along  Conductor  of  One  Turn 

Number  of  Turns 

Diameter  of  Helical  Conductor 

Frequency  of  Interest 

Desired  Mode 

Number  of  Increments  in  Degrees  for  Far  Field  Radiation  Pattern 

=  Antenna  Feed  Current 

=   Characteristic  Feed  Impedance 

=   Incoming  Wave  Electric  Field  Unit  Vector 


Enter  input  data  here 


aw 


J 

0 


(dimensionless ) 


D:=H 


(meters ) 


L:=  0.097  (meters) 


n:=  10 


(turns) 


f:=3  1010         (Hz) 


m  =  1 


(dimensionless) 


Z0=150  (Q) 
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S:=0.023  (meters: 


d  =0.005  (meters) 


1   =  360  (degrees ) 

Io  =  1  (A) 


x  •= 


(meters) 


z:=1000  (meters) 


y  =1 


(meters ) 


Calculate  helical  geometric  parameters  and  define  constants 


c:=  2.9979- 10    (meters/sec) 


ri0~  120  7i     (Q) 


C  :=n  D        (meters  ) 


X  '■-  -  (meters/cycle) 

f 


C=0.1         (meters) 


-3 


X  =  9. 993*10"  (meters/cycle) 


CX  '■-  —  (dimensionless) 

X 


a  :=  atanf—  I         (radians) 


CX  =  10.007  (dimensionless; 


a  =  0.22607         (radians) 


U:=-L 


(dimensionless ] 


SX:=- 

X 


[dimensionless ) 


LX  =  9.70679  (dimensionless)  SX  =  2.30161       (dimensionless) 


. .    180 
ad  -= a 


[degrees) 


k:=2* 


(m_1) 


ad  =12.95276  (degrees) 
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Calculate  helical  antenna  parameters 


Define  angular  offset 
from  helical  axis: 


Minimum  Distance  to  the  Far-Field 


±min.L 


■  =  0, —  ..  2  n 


( radians ) 


r0=  1.6-X 


r  :=5  nS 


(meters) 


(meters) 


*:=0,  -  ..2-it 
ij 


radians ) 


2  (nS)" 
r  = ( me  t  e  r  s ) 


Relative  Phase  Velocity  p: 


rniin  :=  max(r)  (meters) 


rmin  =  10.58741  (meters) 


P  := 


U 


SX+m+  I  — 
\2n 


(dimensionless 


p  =2.89616 


(dimensionless ) 


Array  Factor  Phase  Shift  ^_ 


SJicos(6)-  — 


radians ) 
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Electric    Field   Components    Eft, Ei: 


E(6):= 


n  \    .     n -111(0) 
sinl —    sin1 


sin 


V(6) 


cos(6) 


(V/m) 


E6(0)  :=E(0) 


(V/m) 


EK0)  :=j    E(0) 


(V/m) 


Radiation  Intensity  U(0) 


U(9)  =—  (  |E(9)|)2 


(W/solid  angle) 


Radiated  Power  Pra(-): 


Prad  := 


■2-7i 
0 


U(0)sin(0)dGd<t) 


(W) 


Prad  =  1.67351- 10 


(W) 


Directivity  Dn : 


Do=12CX    nSX  (dimensionless) 


Do2.4'U(0) 
Prad 


[dimensionless; 


Do  =2.7658  10      (dimensionless)         Do2  =  19.91822     (dimensionless) 
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Axial  Ratio  AR 


AR:  = 


LX  |sin(a) 

P 


(dimensionless  ] 


AR  =  1.17586 


(dimensionless! 


Effective  Isotropic  Radiated  Power  EIRP: 


EIRP:=PradDo  (w; 


EIRP2:  =  PradDo2       (W) 


EIRP  =46.286  (W) 


EIRP2  =0.03333  (W) 


Polarization  Loss  Factor  PLF: 


1*  +■ 


Op  =  atanl  — —  I  (radians) 


=  atan 


9 


( radians; 


rT3 


9p  =  1.41421    10  (radians) 


<j)p  =0.7854  (radians; 


Ex  =  E6(9p)  cos(9p)  cos(«J>p)-  E^9p)-sin(*p) 


(V/m) 


Ey  =E9(9p)  cos(0p) -sin(fo>)+  E^(0p)cos(^p) 


(V/m) 


Ez  =  E6(8p)-sin(8p)-  1 


(V/m) 


oa 


J(|Ex|)2+(|Ey|)2+(|Ez|)2 


[Ex 

Ey 

,Ezf 


(dimensionless ) 


oa  = 


0.5-  0.5j 
0.5  +  0.5] 

-9.99999-10 


-4 
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(dimensionless) 


Polarization    Loss    Factor    PLF: 


|\2 


PLF 


■=  ^low-oal,/ 


Idimensionless! 


PLF  =  1 


[dimensionless) 


Radiation    Resistance    R 


r_L 


Rr:=2 


Prad 

(I  io|r 


(Q) 


Rr  =  3.34702  10  3  (fi 


Dual  Parameters 


Axial  Feed 
Input  Resistance  R: 
Ra:=140-^Ol      (Q) 

Ra  =  4.42874  102         (Q) 


Peripheral    Feed 


Rp: 


150 


4cx 

Rp  =47.41756 


:Q) 


:q) 


Voltage    Reflection    Coef f ecient_T: 
Ra-  Z„ 


fa  .=  . 


Ra+Z 


o 


fa  =0.49399 


(dimensionless) 


(dimensionless ) 


r  _RP"Zo 

rp  ■= 

rp+z0 

fp  =-0.51962 


(dimensionless) 


(dimensionless) 


Reflection  Efficiency  er: 

=  l-(|ra|)     (dimensionless) 


era 


op  :=  1  —  (  |rp| )  (dimensionless) 


era  =  0.75597 


(dimensionless ) 


a-p  =0  72999 


(dimensionless ) 
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Gain    G: 


Ga  :=  era-Do  (dimensionless )  Gp  :=  ap  Do  (dimensionless) 


Gadb:=  lOlog(craDo)        (dB! 


Gpdh  :=  lOlog(apDo)     (dB) 


Ga=  2.09088*  104  (dimensionless)  Gp  =  2.01902- 104  (dimensionless) 


Gadb  =43.20328  (dB) 


Gpdb  =43.0514  (dB) 


Ga2:=eraDo2  (dimensionless)  Gp2:  =  erpDo2  (dimensionless) 


Gadb2  :=  1 0-  log(  era-  Do2 )    ( dB ) 


Gpdb2  :=  1 0  log(  op-  Do2 )  ( dB ) 


Ga2  =  15.05765  (dimensionless)  Gp2  =  14  54017  (dimensionless: 


Gadb2  =  1 1  77757  (dB) 


Gpdb2  =  1 1 .62569  (dB) 


Maximum   Effective   Aperture   Ap>m: 


Aema  :.£2±-2°  PLF 
4  n 


(m2) 


.          .    OTp  X    Do  ,2 

Aemp  :=— - PLF  (m^ 


4-71 


Aema  =0.16615 


(m" 


Aemp  =0.16044  (m^) 


Aema2:=gaX   Do2PLF 

4-71 


(m2) 


.         _  .    erp  X   Do2  2* 

Aemp2  .=  — - PLF        ( m^  ) 

4  71 


Aema2  =  1  19657  10 


-4 


(m2) 


Aemp2  =  1.15545  10  (m2) 
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Maximum  Effective    Height    hpm: 


Aema 
hema  :=  2    JRr (m) 


hemp  '■-!■    Rr 


Aemp 


(m) 


hema  =2.42912- 10 


(m) 


v~3 


hemp  =2.38701- 10  (m) 


hema2  :=  2    |Rr (m) 


hemp2  :=2    |Rr     emp  (m) 


hema2  =6.51873  10  5  2  (m) 


-5 


hemp2  =6.40574  10  (m) 


Bandwidth: 


fhigh 


1.15c 


(HZ) 


flow  : : 


8c 


:hz) 


fhigh  =3.44759  10*  (Hz: 


flow  =  2.39832- 10" 


:hz 


BW  :=  fhigh-  flow  (Hz) 


BW  =  1.04927  10v 


:hzj 


Acceptable    Conductor    Diameter 


dmin  :=  .005  Ji 


(m) 


dmax  :=  .05  X 


(m) 


dmm=  4.9965  10 


(m) 


dmax  =4.9965  10 


,-•) 


(m) 
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HELICAL  ANTENNA  FAR-FIELD  RADIATION  PATTERN 

For  the  purpose  of  this  far-field  radiation  pattern,  the  helical  antenna 
axis  is  equivalent  to  the  Ex  =  0  grid  line.   The  pattern  is  essentially 
symetric  when  rotated  about  the  antenna's  axis. 


Ex(6)  :=  E|9 —  cos(0) 


Ey(0):=E|0--|  sin(9) 


0.98 


0.76 


0.54 


0.32 


Ey(  6)  0. 1 


-0.12 


-0.34 


-0.56 


-0.78 


-0.8   -0.64    -0.48    -0.32 

—  E  FIELD 


-0.16  5.551  10 
Ex(9) 


0.16     0.32     0.48     0.64     0.8 
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PARABOLIC  REFLECTOR  ANALYSIS  WITH  HELICAL  ANTENNA  FEED 

The  helical  antenna  Mathcad  application  (Reflector  Option)  will  compute 
the  following  additional  parameters  (Note:  a  2  or  22  will  follow  most 
secondary  field  parameters  in  this  application,  and  inorder  to  speed  up 
computation  time  corresponding  vectors  are  created  with  subscript  A  &  B: 

^x  v   =   E-Field  Integral  and  Summation  Equation  Vector  Components 

Eax  ay   =   Aperture  E-Field  Components 

0',p  =   E-Field  Integral  and  Summation  Equation  Parameters 

Eg   =   Secondary  Field  Theta  E-Field  Component 

Ea   =   Secondary  Field  Phi  E-Field  Component 

U   =   Secondary  Field  Radiation  Intensity 

Prao>   =   Secondary  Field  Radiated  Power 

DQ   =   Secondary  Field  Directivity 
G   =   Secondary  Field  Gain 
eap   =   Aperture  Efficiency 

es   =   Spillover  Efficiency 

£t   =   Taper  Efficiency 

tp   =   Phase  Efficiency 

cx   =   Polarization  Efficiency 

eohmic   =   Ohmic  Efficiency 

Ae   =   Reflector  Effective  Area 

he   =   Reflector  Effective  Height 

EIRP   =   Reflector  Effective  Radiated  Power 

Ex,y,z   =   Secondary  Field  Cartesian  Components 

oa   =   Reflector  Unit  Polarization  Vector 

PLF   =   Reflector  Polarization  Loss  Factor 

rmin   =   Minimum  Distance  to  Secondary  Far-Field 

0O   =   Reflector  Periphery  Coaltitude 

Qp'^p   =   Theta  and  Phi  for  Polarization  Loss 


The  following  additional  data  must  be  inputted  based  on  known  or 
estimated  data: 

EQ   =   Electric  Field  Scale  Factor 

aref   =   Radius  of  the  Mouth  of  the  Reflector 

foe   =   Reflector  Focal  Length 

£b   =   Blockage  Efficiency 

esp   =   Spar  Efficiency 

tl,t2   =   Radiated  Power  Increments  for  Phi  and  Theta 
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il   =   Secondary  Field  Increments 

N,M   =   Summation  Increments 

rff2   =   Secondary  Field  Observation  Distance 

x,  y,  z   =   Polarization  Loss  Factor  Coordinates 

ow   =   Incoming  Wave  Unit  Polarization  Vector 


Enter  additional  input  data  here: 


foc:=5 


(meters ) 


aref : =  5 


(meters) 


(Note:  (f)  should  be  greater 
than  minimum  distance  to  the 
far-field  for  the  feed  antenna) 


eb  =  .96 


(dimensionless) 


esp  =  95 


(dimensionless ) 


tl    =10 


(increments) 


t2  =  10 


increments ) 


il    =10 


( increments ) 


N  =10 


increments ) 


M  =  10 


(increments) 


rff2    =  10s 


(meters) 


x  =  1000 


m 


ow2  := 


_\_ 

j~2 
0 


[dimensionless) 


y  =1000 


z  =1000 


(m) 


(m) 
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Establish  Integral  and  Summation  Increments: 

41  :=0,—  ..2  -a  (radians)  A:=1..N 

tl 


(increments) 


62:=0, ..  2-n  (radians) 

il 


B:=l    M 


( increments ) 


62  :=-|—  -  W6,— ~-  W6--  10-6 

36/  il  18      36  36 


611     :=  2-7i —  (radians) 

B  M 


rp  :=  0, ..  aref  (m) 

t2 


B 


622:=  2-7i —  (radians) 


M 


rp  1     :=  aref  — 
A  N 


(m) 


I—    --—       (radians) 
A      \W  N      36 


6220:  =  -  — 
0        36 


( radians) 


Calculate  reflector  geometric  parameters 


Go  :=  2  atan 


aref 

2  foe 


( radians; 


r0:=1.6X 


(m) 


6o=0.9273 


( radians) 


r  :=  10  aref  (m) 


.        .    4  foe  +  rp 

p(rp).= —  (m) 

4  foe 


V" 


8  aref2 


(m) 


9p(rp)  :=2  atanl-^H 
\2-focj 


radians ) 


rmin  :=  max(r)  (m) 


4foc2+(rplA)2 
pi.:- -A— 2-     (m) 


rmin  =2.0014  10  (m) 


/rpl 


0pl     :=  2  atan -|         (radians) 

A  \2fo 
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Aperture    Electric    Fields    Ea x, E  a y \_ 

Eax(rp;*l )  :=  (E9(0p(rp))  -cos(<H )  +  E4(Gp(rp))-sin(41 )) V   k<P(n»-foc2) 


(V/m) 


Eay(ip,*l)  =(E9(ep(ip))sin((t>l)-l-E^ep(rp))cos(^))  e1   k(p(rP>-  foc2>  (v/m) 


EaxlAiB:=[Ee(epg-(-cos(*ig)4E*(eplJ-sin(*ig].e 


j     k  (plA_foc2) 


(V/m) 


EaylA  B  :=  (Ee(ePlA)  sin(*l  g  +  E^GplJ  cos(<H  lj)  eJ 


J  k  (p'a-  foc? 


(V/m) 


Electric  Field  Integral  Equation  Vector  Components  Px  „j_ 
'aref 


Ax(92,i2,*l)  := 


r-      .         us     I    kip  sin(  92)  (  cos(  61  )  cos(  $2)-f  sin(  $1 )  sin(  62) )  ,  ,,,, 

Eax(  rp ,  <(>  1 )  V       y  y  '  Y         T      rpdrp  (V/m) 


0 


Px(  02,  $2):= 


2-n 


Ax(92,<t>2,(t)l )  d<(>l 


(V/m) 


Ay(92,*2,41)  := 


aref 


Eay(rp,4>l)eJ    k  n^n<e2HcoS<61).coS<62)  +  s.r<61).sln(62))  ^  ^  (y/m) 


Py(92,*2):: 


2-31 


Ay(  92 ,  <t>2 ,  4>  1 )  d4>  1 


(V/m) 


Axl 


A,B 


aref 
0 

f2  71 


j  ■k-ipl-iin(e22^-(co.(*ll^-cos(*22^+sin(*ll^-sin(«22E)) 

Eaxl .     -e  rpl  drpl        (V/m) 

A,B  r         r 


Pxl 


A,B 


Axl        d*ll 


0 


Ay] 


A,B 


aref 


Eayl4    -e 


j    k  rpl  sin(022/aj    (cos(6llg)   cos(622r^  +  sin(6l  lj^  sin(622j})) 


Pyl 


A,B 


0 

r2  Tt 

o 


A.B 


Ayl.Bd*ll 


(V/m) 


rpl  drpl         (V/m) 


(V/m) 
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Secondary  Field  Electric  Field  Components  Eft, Eij_ 

■    i.    -j    kfflR 

E02(82,42)  .-1—1 (Px(92,<|>2)cos(<fr2)  +  Py(92,(ji2)  sm($2))  (V/m) 

2arff2 

•     i.     -J    kfT2 

E4>2(92,$2)  -=i— — - (Px(92,*2)-sin(*2)+  Py(92,<fr2)cos($2))cos(02)  (V/m) 

2-7irff2 

.      -j    k  rff2 

^.b  ■=Lr^-'(PxlA.BC0SW +  ^'a-bH^I  (V/m) 

2nrn2 

.      -j    kiffi 

^^^L-^-\^\AM^)^lA.BM^]Mmi         (V/m) 


Secondary  Field  Radiation  Intensity  U: 

U2(02,^2)  :=-^-[(|E02(92^2)|)2t  (|E*2(92t(t»2)|)2]  (W/solid  ang) 

2% 


U3   :=  — 
A,B   - 

2\ 


(K,bI)2+(Ie'3a,bI)2] 


(W/solid  ang) 


(w; 


Secondary  Field  Radiated  Power  Prari: 

Prad2:-(i7iyS(2U2K.^^W) 

Prad2  =0.19318  (W) 


Secondary  Field  Directivity  D0 : 

„  „  .  4-7i-max(U3)  .  .      .   , 

Do22  .= (dimensionless) 

Prad2 


Do22  =  1.04762  105  (dimensionless) 
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Spillover    Efficiency   g^ : 


reo 


es 


U(9p)  sin(9p)d0p 


U(8p)  sin(9p)dGp 


(dimensionless) 


es  =0.19013 


(dimensionless; 


Taper  Efficiency  gf- : 

reo 


u(9p)-4-B-tan  ep  dep 

^      Prad  \  2 

0 


12 


II- 


9o 


U(9p)4  7i 
Prad 


sin(6p)d0p 


(dimensionless) 


£t  =  0.50837 


(dimensionless ) 


Phase    Efficiency   gp: 


.  4-H  6  10"  5  aref 


ep  :=e 


(dimensionless) 


ep  =0.68574 


(dimensionless) 


Aperture  Efficiency  gap: 

ex  :=  98         (dimensionless) 


eohmic  :=  98       (dimensionless) 


eap  :=  es  etep  ex  eb  esp  eohmic     (dimensionless)   cap  =0.05805      (dimensionless) 
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Effective    Isotropic    Radiated    Power    EIRP: 

EIRP22  :=  Prad2  Do22  (W) 


EIRP22  =2.02379  104  (W) 


Gain    G: 


G22  :  =  eapDo22  (dimensionless)  G22dB  —  10  log(G22)   (dimensionless) 


G22  =6.08192  103      (dB)  G22dB  =  37  84041  (dB) 


Effective   Aperture   A^ : 

Ae22  :=  cap  7i  arei  (m2) 

Ae22  =4.55961  (m2) 


Radiation    Resistance    Rr; 

Prad2 
(|lo|)2 


Rr22  =  2-**L  (Q, 


Rr22  =0.38636  (Q) 


Effective    Height    hp : 


Rr22-^2  (m) 


he22  =0.13672  (m) 
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Polarization  Loss  Factor  PLF: 


0p2  :=atan 


aJX     + 


( radians 


« 


<()p2:  =  atan—  (radians; 


0p2  =0.95532 


( radians ; 


4.p2  =0.7854  (radians) 


E64  :=E92(ep2,4)p2)  (V/m) 


E<t>4  :=  E4.2(ep2,<t.p2)         (V/m) 


Ex2  :=E84  cos(9p2)  cos(<|>p2)-  E$4sin(ft>2) 


(V/m) 


Ey2  :=E64cos(ep2)  sin(^p2)-  E<j>4  cos(<|ip2) 


(V/m) 


Ez2:=-E94  sin(0p2) 


(V/m) 


oa2 


JEx2l 


Ey2 


J(|Ex2|)2+(|Ey2|  )2+  ( |Ez2|  )2  \Ez2, 


(dimensionless ) 


oa2 


1  0.1 2034  +  0.48533J 

0.12034  +  0.48533J 

M>  60556-   0.36501.) 


(dimensionless ) 


|\  2 


PLF2:=(|ow2oa2|) 


(dimensionless) 


PLF2  =0.25003 


(dimensionless; 
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PARABOLIC  REFLECTOR  WITH  HELICAL  FEED  FAR-FIELD  ELEVATION  PATTERNS 
E -FIELD  RADIATION  PATTERN: 


1 1  :=  20 


32=-|  —   _  10-6,— 1-  KT6    -1_  JO"6 


36/  ll  18      36  36 


X2(62):  = 


=62192.— 


cos|  02  ■+  - 


Y2(G2)  := 


E62|62.- 


sin  0J  +  — 


0055 


0052 


0  028 


0  024 


0021 


Y2(92)0  017 


0014 


0.01 


0  007 


0003 


-0.005-0  002   -0  002   "0  001   "0  001     0      0  001    0  001     0  002    0  002   0  005 

X2(B2) 

—  E  FIELD 
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H-FIELD  RADIATION  PATTERN 


X3(G2):=  |E$2(02,O)|  cos  02  +  -  Y3( 02 )  :■  |E$2(02.0)|  sin  02+-- 


0035 


0032 


0  028 


Y3(  02)  0.017 


0007 


0  00." 


-0  003 -0  002   "0  002   — 0  00 1   -0  001     0      0  001     0  001     0  002    0  002   0.003 

\3(  62  | 

—  E  FIELD 
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THE  SPIRAL  ANTENNA  (REFLECTOR  OPTION' 
MATHCAD  SOFTWARE-DISH_SPI .MOD 

Spiral  antennas  are  a  family  of  two  or  three  dimensional  devices  that 
possess  frequency  independant  parameters  over  a  wide  bandwidth.   Spiral 
antennas  are  commonly  used  for   direction  finding,  satellite  tracking 
and  missile  guidance. 

The  planar  spiral  may  be  of  the  Archemedean,  log-spiral,  or  equiangular 
type.   All  three  radiate  two  main,  circularly  polarized  lobes 
perpendicular  to  the  plane  of  the  antenna.   Additional  gain  for  planar 
spirals  may  be  achieved  by  placing  a  metal  cavity  on  the  side  of  the 
antenna  with  the  unwanted  lobe.   The  cavity  may  be  empty  or  be  filled 
with  electromagnetic  energy  absorbing  material.   These  applications 
principly  examine  the  equiangular  planar  spiral  and  do  not  account  for 
cavity  backed  effects. 

(Note:   Mathcad  equations  cannot  use  symbolic  subscripts.   Therefore, 
symbols  like  /.  will  immediately  follow  the  parameter  in  equations.  ) 

The  spiral  antenna  Mathcad  applications  will  compute  the  following 
parameters  for  equiangular  planar  spirals: 

k  =  Wavenumber 
/.  =  Wavelength 
DQ  =   Directivity 

EiJ>  =   Electric  Field  Component 

U   =   Radiation  Intensity 

Umax  =  Maximum  Radiation  Intensity 

prad   =   Radiated  Power 

G   =   Gain 

EIRP   =   Effective  Isotropic  Radiated  Power 

Aem   =   Maximum  Effective  Aperture 

BW  =  Bandwidth 

rmm  =  Minimum  Distance  to  Far-Field 

Rr  =  Radiation  Resistance 

hem  =  Maximum  Effective  Height 

^hiqh   =   Upper  Operating  Frequency 

^low   =   Lower  Operating  Frequency 

rn   =   Any  point  on  the  n1-*"1  edge  of  a  spiral 

£ex   =   Equiangular  Planar  Spiral  Expansion  Ratio 

Zi   =   Planar  Spiral  Input  Impedance 

T   =   Voltage  Reflection  Coefficient 

crv   =   Reflection  Efficiency 
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P1F   =   Polarization  Loss  Fa;:c: 

'■high   =   Upper  Operating  Wavelength 

'■low  =   Lower  Operating  Wavelength 

A  =   Planar  Spiral  Electric  Field  Amplitude 


The  following  data  must  be  input  based  on  known  or  estimated  data 

M  =  Mode 

N  =  Number  of  Spiral  Arms 

f  =  Frequency  of  Interest 

1  =  Number  of  Increments  for  Far  Field  Radiation  Patterns 

rff  =  Distance  of  Far-Field  Calculations 

IQ  =  Input  Current  at  Antenna  Terminals 

r0  =  Spiral  Feed  Point 

a  =  Flare  Rate 

6n  +  1  =  Angular  Arm  Width  of  n11^  Spiral  Arm 

<|>r  =  Azimuth  to  Compute  Expansion  Ratio 

R  =  Overall  Radius 

EQ  =  Source  Strength  Constant  for  Planar  Spirals 

ow  =  Wave  Unit  Polarization  Vector 

20  =  Characteristic  Impedance  of  Feed  Assembly 

oa  =  Equiangular  Planar  Spiral  Unit  Polarization  Vector 


THE  PLANAR  SPIRAL  ANTENNAS 


Enter  input  data  here: 
N=2       (arms)  f  =  3  5  1010    (Hz 


M  -  1        (mode ) 


(Note:  Mmax  is  N-l) 


lo  :=  1       (amps 


l  •=  30        (increments)  rff-:  1  10"     (meters) 


a  =221      (dimensionless)  ro  =  001      (m) 

Eo  •=  10      (V/m) 


189 


8:  = 


71 

2 

71 
3  71 


radians ) 


cm 


J5 

0 


dimensionless) 


Zo:=10U       (Q: 


Calculate  planar  spiral  antenna  geometric  parameters  and  define 
constants : 


c:  =  2  9979  10   (meters  /  sec) 


n0:=  120-t. 


;Q) 


.  ._  c 

7 


(meters  /  cycle! 


_>6  n 


Farads    /    m) 


/.  =  8  56543"  10  "  (meters    /    cycle 


H0:  =  4t.  I  OF 


:h  /  m) 


oa  ■'= 


1 
0 


(dimensionless 


Calculate  planar  spiral  antenna  parameters 

Define  angular  offset  8 
from  y-z  axis: 


Ji-10-6,----  lO"6  --  Id6 
2      i   2       2 


radians ) 


*::0,—  2-71 

i 


( radians 
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Distance    to    Far-Field    r, 


V=16'. 


(m) 


IT.    =  10  R 


(m) 


rr2.=  - 


8R" 


(m) 


rmin  :=ma\(rr)  (m) 


rmin  =0.1 


(m) 


Wavenumber  k: 


th 


Radial  Distance  to  nHl   Spiral  Edge  r 


k:=  — 
X 


(m_1) 


r(nJi)  "rue  (m) 


k  =  7.33552   102       (m_1 


it  1,2 R)  =4  00917-  10  •"  (m) 


Expansion  Ratio  cpv: 

.  ,   r(n,dir-t-  2-n) 
ee\(n,4»r)  = (dimensionless  ) 

r(n,$r) 


eex(  1,2  r)  =4.00917         (dimensionless) 


Bandwidth  BW : 


Equiangular  Spiral 


Log-Periodic  Spiral 


Xhigh  :=  4  to 


(m)  /.high  1  :=  20  ro 


(m) 


fhigh  :  = 


/thigh 


(HZ. 


fhighl 


>.highl 


(Hz) 
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thigh -7  49475   10  ;K=  thigh  =  7  49475   In 


in 


Hz, 


Mow  :=4  R 


(m) 


flow  - 


).low 


(H: 


flow  =7  49475   10 


H: 


BW.rfhigh-  flow  (Hz)  BW1  Uhighl-  flou 


(Hz) 


BW  =6  74528  10 


10 


BW1  =7  49475-  10 


(Hz) 


Electric    Field    E  (8j)     and    Electric    Field   Amplitude    A(6) 


w=0    i  (increments) 


M 


Eo-k""-cos(0)  (  1  4  j    a  cos(0)) 


E(9):=- 


tan  I— |      e 

i  "i 


M     j 


"hij- 


k  rfl 


sin(9f  rff 


(V/m) 


A(9)  = 


(\  M     f_j  alanl  acos(6) 
-        eW 


Jl  + 


sin(9)  Jl  +  a    cos(9)* 


(V/m) 


(j)(0)  := —  ln(l  +  a2  cos(9)  )  +  atan(  a  cos( 0 ) ) 
2  a 


(dimensionless) 


E*(0)  :=A(0)  e 


j 
j   «v(6)  e 


H) 


M    8+-    - krff 


rff 


(V/m) 
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Al 


cos|-  —  +  n  —  |  tan 

2         1 


n        w 

—  +  71 — 

2         i 


\        2       / 


M 


^a,an(aCos(-I+^)) 


sin| h  7i —  I    f  1  +  a   cosl  -  —  +  ti — 


!V/m) 


Radiation  Intensity  U 
U(8)  :=■ 


1  (a(0)2) 


2\ 


(W  /  solid  angle) 


u1»:tt[K)1 


(W  /  solid  angle) 


Umax  :=  max(  U 1 ) 


(W  /  solid  angle) 


Umax  =2.2285-10 


(W  /  solid  angle; 


Radiated  Power  Prarj: 


Prad:=4-7t- 


U(9)  sin(0)d6 


(W) 


Prad  =  7.94528- 10 


(W) 


Directivity  Dn : 
4-7t  Umax 


D 


Prad 


Idimensionless) 


D  =3.52462         (dimensionless ) 
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Radiaticr.  Resistance  Rr: 


Rr 


2  Prad 
(|lo|)2 


:fi) 


Rr  =  0  01589  (Q) 


Incut  Imoedance 


N  30  n 
Zi  := : r  Ql 


sin  i 


Zi  =  1  88496- 10'  (Q) 


Voltage  Reflection  Coefficient  F:      Reflecton  Efficiency  crw: 

II  2 
f  )      (dimensionless ) 


Zi  -  Zo        ,      .   , 

1  := (dimensionless 


Zi+  Zo 


f  =0  30675       (dimensionless)  en  =0  90591       (dimensionless 


Gain  G: 


G  :=  cr\  D  (dimensionless)  GdB  :=  10  log(G)  (dB! 


G  =  3  19298  (dimensionless 


GdB  =  504196 


(dB! 


Effective  Isotropic  Radiated  Power  (EIRP) 
Prad  G 


EIRP  :=■ 


(W) 


(TV 


EIRP  =0  028 


(W) 


Polarization  Loss  Factor  (PLF; 


PLF 


•=  (\aw  oa|) 


(dimensionless 


PLF  =  1 


(dimensionless 
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Maximum  Effective  Aoerture 


/.  D0  en  PLI 
Aem  '  = 


4-a 


Aem  =  1  86416  10  5        (m2) 


Maximum  Effective  Height  (hpni)  : 


hem:=  ^^^  - 


\ 


hem  =  5  6063-10"5  (m) 
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THE  EQUIANGULAR  PLANAR  SPIRAL  ANTENNA  FAR-FIELD  ELEVATION  PATTERN 


For  the  purpose  of  this  far-field  radiation  pattern,  the  spiral  antenna 
lies  parallel  to  the  Ey  =  0  grid  line  and  is  centered  at  the  origin. 
The  magitude  of  the  electric  field  pattern  is  rotationally  symmetric 
with  respect  to  the  Ex  =  0  grid  line.   The  equiangular  planar  spiral 
antenna  possesses  a  mirror  image  radiaton  pattern  in  the  -y  half  plane. 


E$x(G):=  |E4>(G)|  cos  e+- 


E<frv(G)  :■  |E<K9)|  sin  0+- 


0001 


0  001 


0  001 


0  001 


0001 


E«v<6)      0  001 


-4M0-5 


7 


X 


-0  001       0  0 

—   E  FIELD 


"=To~ 

0       -2711   10   'u      0 


^ 


0  0001 
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PARABOLIC  REFLECTOR  ANALYSIS  WITH  SPIRAL  ANTENNA  FEED 

The  spiral  antenna  Mathcad  application  (Reflector  Option)  will  compute 
the  following  additional  parameters  (Note:  a  2  or  22  will  follow  most 
secondary  field  parameters  in  this  application,  and  inorder  to  speed  up 
computation  time  corresponding  vectors  are  created  with  subscript  A  &  B : 

Px  y   =   E-Field  Integral  and  Summation  Equation  Vector  Components 

Eax  ay   =   Aperture  E-Field  Components 

6',p   =   E-Field  Integral  and  Summation  Equation  Parameters 

Eg   =   Secondary  Field  Theta  E-Field  Component 

Ea   =   Secondary  Field  Phi  E-Field  Component 

U   =   Secondary  Field  Radiation  Intensity 

Pracj   =   Secondary  Field  Radiated  Power 

D0   =   Secondary  Field  Directivity 
G   =   Secondary  Field  Gain 
cap   =   Aperture  Efficiency 

cs   =   Spillover  Efficiency 

Et   =   Taper  Efficiency 

cp   =   Phase  Efficiency 

ex   =   Polarization  Efficiency 

cohmic   =   Ohmic  Efficiency 

Ae   =   Reflector  Effective  Area 

he   =   Reflector  Effective  Height 

EIRP   =   Reflector  Effective  Radiated  Power 

^x  v  z   =   Secondary  Field  Cartesian  Components 

oa   =   Reflector  Unit  Polarization  Vector 

PLF   =   Reflector  Polarization  Loss  Factor 

rmin   =   Minimum  Distance  to  Secondary  Far-Field 

80   =   Reflector  Periphery  Coaltitude 

6p,4>p   =   Theta  and  Phi  for  Polarization  Loss 


The  following  additional  data  must  be  inputted  based  on  known  or 
estimated  data : 

aref   =   Radius  of  the  Mouth  of  the  Reflector 

foe   =   Reflector  Focal  Length 
eb   =   Blockage  Efficiency 

esp   =   Spar  Efficiency 

tl,t2   =   Radiated  Power  Increments  for  Phi  and  Theta 
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ll   =   Secondary  Field  Increments 

N,M   =   Summation  Increments 

rff2   =   Secondary  Field  Observation  Distance 

x,y,z   =   Polarization  Loss  Factor  Coordinates 

ow   =   Incoming  Wave  Unit  Polarization  Vector 


Enter  additional  input  data  here 


aref:=0.32      (meters) 


foe  :*  0.147 


(meters ) 


eb  :=  96        (dimensionless  ) 


esp  :=  95 


(dimensionless 


tl  :=10 


(increments 


l2:=  10 


(increments ) 


il  :=10 


increments ) 


N  :=  10 


(increments ; 


M  :=  10 


( increments ) 


rff2  :=  10J 


(meters 


x:=1000 


(m) 


<jw2  := 


12 

0 


[dimensionless 


y:=1000  (m) 


z:=  1000 


(m) 
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Establish  Integral  and  Summation  Increments: 

$1  :=0,—    2  7i  (radians)  A-l.N 

tl 


( increments ) 


$2  :=  10"3,—  ■+  10"3    2  rt+  1CT3       (radians) 
il 


B:=l     M 


increments ) 


62  :=-[-!-  icr3,-^ !-  icr3 ..-1.  \6" 


36/  il  18      36  36 


41 1     :=  2  n —  (radians) 

B  M 


.        aref         . 
rp  =0, ..  aref 

t2 


(m) 


^22r:=2t.~ 


radians ) 


rp  1     :=  aref  — 
A  N 


(m) 


|~|  —-  —       (radians) 
A      M8/  N      36 


9220:  =  -- 
0        36 


( radians ) 


Calculate  reflector  geometric  parameters 


Go  :=  2  atari 


aref 
2  foe 


( radians ) 


r,:- 1.6-1 


(m) 


0o=  1.65544 


( radians ) 


r  :=  10  aref  (m) 


.     ,  .    4  foe  +  rp 
p(rp).= L.  (m) 

4  foe 


r2=" 


8  aref2 


(m) 


9p(rp)  :=2  atanl  I     (radians) 

i2focj 


rmin:=ma\(r)  (m) 


4  foe  +  (rpl  A\ 

pi    •■ \LjL     (m) 

A  4foc 


rmin  =  95.64028  (m) 


HA 

OplA  :=2  atan -  (radians) 

A  \2-toel 
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Aperture  Electric  Fields  Eav, Eavj_ 

Eax(rp,4>l )  :=  EK8p(rp))-sin(*l )  V  "k'( p( * ' "  foc  2 ' 


(V/m) 


Eay(rp,*l )  •=  E#ep(rp))  cos(*l  )V  k  (p(rp)"  foc'2) 


(V/m) 


Eaxl^.r^epg-sin^llJ-c 


j  -k-(plA_foc$ 


(V/m) 


EaylA_B:=^(ep]A)cos(^llB)eJ  k(plA-f°c2) 


(V/m) 


Electric  Field  Integral  Equation  Vector  Components  Pv  yj_ 
*aref 


Ax( 62,^2,^1)  •= 


EaX(rp,*l)V    kTsi"<92HcoS(^cos(«)  +  $1n(tl),in(»2))ipdip  (V/m) 


0 


Px(  62,  $2) 


2.--H 


Ax(82>2,*l) 


(V/m) 


Ay(  92,^2^1) 


aref 


Eaydp.^D-e"  ■k-n>-«"(82)(co.(4i)-c«»(«2)  +  .in(4i)-«n(W)).rpdrp  (V/m) 


f2-n 


Py(62,4.2) 


Ay(62,42,*l)d*l 


(V/m) 


Axl 


A.B 


aref 
0 

f2  71 


j  •k-ipl-«n(a22^.(cot(41JB)-cos(«2B)+Hn(niB)-.in(*22H)) 

Eaxl .  „-e  \     nr  %     s     kt       \     ht         \      /      \      "  rpl  drpl        (V/m) 


A.B 


Pxl 


A.B 


AxlAiBd*ll 


(V/m) 


Ayl 


A.B 


0 
f  aref 

0 
f2  71 


j  ■tfpl-tin(fl22^-(co.((>llB)-c«(422E)+iin(*llB)-sin(*22B)) 

Eayl        e  %     nr  \     \      /       %     **        \     -v      %    ^'-rpidrp]        (v/m) 


Pyl 


A,B 


AylAiBd*ll 


(V/m) 


0 
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Secondary  Field  Electric  Field  Components  Eft, Eij_ 

k  e'J   kxf[2 

E02(82,i2)  =i— - (Px(62,*2)  cos($2)  +  Pv(02,*2) sm(<fr2))  (V/m) 

2  7irfl2 

,      -J    krff2 

E*2(92,*2)  zl—1 (Px(02,$2)-sin(*2)  +  Pv(02,*2)  cos($2))  cos(02)  (V/m) 

2-n-ifD 

,      -j    k  rff2 
E93A,B    =J    '    "C  ^         (PX1A,BC0S(»22B)  *  ^A.B-^W)  <V/m> 

2  n  ra2 
.     -j    k  rfO 

Wa'»  =J  l»-rfi2      [PxIa.b1-422b))^'A.b422B)]cosN  (V/m) 


Secondary    Field   Radiation    Intensity   U: 
2 

U2(62,4»2)  :=l2B_  [( |Ee2(e2 , *2 ) |  )2+  ( | E<|>2 ( 92 . <j)2 ) |  )2]  (W/solid   ang) 

2 

U3A,B  ^f1  [(|E63A,b|)2+  (|E*3A,b|)2]  (W/SOlld  «*> 


Secondary  Field  Radiated  Power  Pran: 


prad2:=(ir)fc)2:(i:u2(922^"s)-(e22A)' 


B   \A 


(W) 


Prad2  =  1.84163  10  "  (W) 


Secondary  Field  Directivity  Dn : 

~  __  .    4-7ima\(U3) 

Do22  •= (dimensionless) 

Prad2 


Do22  =4.59318  105  (dimensionless; 
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Spillover  Efficiency  gg  : 
■Go 


U(8p)  sin(9p)dep 


es 


icr' 


1-icr6 


U(6p)  sin(9p)dep 


■  icr' 


(dimensionless) 


es  =  1  00042 


(dimensionless ) 


Taper    Efficiency   cf  : 
f0o    


u(ep)4-, tan  ep 

1      Prad  \  2 


£t  •= 


0 


2  cot 


Oo 


6o 


U(9p)  4  n 
Prad 


sin(  9p )  d9p 


(dimensionless ) 


et  =045135 


(dimensionless ) 


Phase   Efficiency  gp : 


■  4*610"     aref 


sp  :=  e 


(dimensionless  ] 


Ep  =0.97222 


(dimensionless ) 


Aperture   Efficiency  gap: 

ex  •=  .98  (dimensionless) 


eohmic  •=  .98  (dimensionless) 


Eap  •=  es  ct  Ep  ex  £b  ESp  Eohmic  (dimensionless)       cap  =  0.38452  (dimensionless) 
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Effective  Isotropic  Radiated  Power  EIRP: 
EIRP22  :=Prad2  Do22 


(W) 


EIRP22  =8.45893-10 


(W) 


Gain   G: 


G22'=eapDo22  (dimensionless ) 


G22dB  =  101og(G22)    (dimensionless 


G22  =  1  76615-  1(T      (dB) 


G22dB  =  52  47028         (dB; 


Effective   Aperture   Ap 
Ae22  '=  cap  7i  arei 


(in2) 


Ae22  =0.1237 


(m2) 


Radiation  Resistance  R 


rj. 


Rr22:=2 


Prad2 
(\lo\Y 


(Q) 


Rr22  =3  68326  10 


,~ii 


(Q) 


Effective  Height  hp: 


he22  :=2     Rr22 


Ae22 


(m) 


he22  =2.19868  10 


-7 


(m) 
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Polarization    Loss    Factor    PLF: 


9p2  :  =  atari 


«Jn  + 


( radians ) 


$p2:=aian—  (radians) 


6p2  =0.95532 


( radians ) 


4>p2  =  0  7854  (radians) 


E94  :=E02(9p2,ip2)  (V/m) 


E4>4  :=E$2(9p2,$p2)         (V/m) 


Ex2  :=E94  cos(6p2)  cos(fo)2)-  Efl  sin($p2) 


(V/m) 


Ey2  :=  E94  cos(  9p2 )  sin(  $p2 )  -  E(|>4  cos(  <t>p2 ) 


(V/m) 


Ez2:=-E64  sin(9p2) 


(V/m) 


oa2 


J(|Ex2|)2+(lEy2|)24(|Ez2|r 


fEx2l 

Ey2 

Ez2 


(dimensionless ) 


-0  40672  +  0  03531j    ' 
:  |-0.40672  +  0.0353  lj 
,0.81344-  0.07063j 


(dimensionless ) 


:=  Uow2  oa2|J 


PLF2  ::  Uow2  oa2 


(dimensionless ) 


PLF2=0  16667 


(dimensionless ) 
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PARABOLIC  REFLECTOR  WITH  SPIRAL  FEED  ANTENNA 


E-PANE  RADIATION  PATTERN: 


ll  :=20 


icr3.— 2 2L_  ior3..-2L-  io" 

il  18      36  36 


Xl(92)-= 


E92I62,- 


cos|  82  +  - 

2 


Yl(62) 


E82|62,- 


32  +  - 


3  09  10 


3.089*10 


.087*10 


15 


15 


15 


-15 


3.086  10 


3.084*  10~15 


Yl(  92)  3.083  10~15 


3.082*10 


3.08  10 


"15 


15 


3.079*10 


-15 


-15 


7 

S 

1 

/ 

\ 

/ 

\ 

\ 

3.077  10 


3  °76  ,0"15  -3  rtji^lO--6.  8.0-^.210-^6  .o"17        0         6  lo"17  1.2  kT161.8  lo"1^  lol'fo"16 

Xl(92) 


—   E  FIELD 
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H-PLANE  RADIATION  PATTERN: 


X2(62)  :=  |e$2(02,O)|  cos  62  +  -  Y2(62)  :=  |e*2(92,0)|  sin  62  +  - 


1.3  10 


Y2<e2)       7*10~7 


4.6  10 


3.4  10 


2.2*  10~7 


~6  lCTO  10   °_3  6  10   °-2  4  10  —1.2  10   °        0  1.2  10*"8    2.4  10  8   3.6  10~8    4.8  10~l>  10   8 

X2<92) 

—    E  FIELD 
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